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Abstract 

In this paper we analyze the image formation problem for undulator radiation 
through an optical system, accounting for the influence of the electron beam emit- 
tance. On the one hand, image formation with Synchrotron Radiation is governed 
by the laws of Statistical Optics. On the other hand, the widely used Gaussian- 
Shell model cannot be applied to describe the coherence properties of X-ray beams 
from third generation Synchrotron Radiation sources. As a result, a more rigorous 
analysis of coherence properties is required. We propose a technique to explicitly 
calculate the cross-spectral density of an undulator source, that we subsequently 
propagate through an optical imaging system. At first we focus on the case of an 
ideal lens with a non-limiting pupil aperture. Our theory, which makes consistent use 
of dimensionless analysis, also allows treatment and physical understanding of many 
asymptotes of the parameter space, together with their applicability region. Partic- 
ular emphasis is given to the asymptotic situation when the horizontal emittance is 
much larger than the radiation wavelength, which is relevant for third generation 
Synchrotron Radiation sources. First principle calculations of undulator radiation 
characteristics (i.e. ten-dimensional integrals) are then reduced to one-dimensional 
convolutions of analytical functions with universal functions specific for undulator 
radiation sources. We also consider the imaging problem for a non-ideal lens in pres- 
ence of aberrations and a limiting pupil aperture, which increases the dimension of 
the convolution from one to three. In particular we give emphasis to cases when the 
intensity at the observation plane can be presented as a convolution of an impulse 
response function and the intensity from an ideal lens. Our results may be used in 
practical cases as well as in benchmarks for numerical methods. 
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1 Introduction 



The majority of experiments based on the use of X-rays are carried out at Syn- 
chrotron Radiation faciUtics, as very high brilhance is achievable by means of 
undulator devices installed in storage rings. Guiding the photons from the exit 
of an undulator to the specimen position requires the development of optical 
beamlines whose main task is to re-image the undulator source to any plane 
of interest. To deal with the image formation problem, one should account for 
the fact that Synchrotron Radiation constitutes a random stochastic process. 
In fact, the shot noise in the electron beam causes fluctuations of the elec- 
tron beam current density. These fluctuations are random both in space and 
time. As a result, the radiation field produced by the electron beam can be 
described in terms of a phasor with random amplitudes and phases and, in all 
generahty, the laws of Statistical Optics must be applied to solve the image 
formation problem. In this paper we study the image formation problem with 
undulator radiation beams based on Statistical Optics. In this framework, 
the basic quantity characterizing Synchrotron Radiation sources is the second 
order correlation function of the fields at two observation points on a given 
transverse plane identified by the coordinate Zg along the optical beamline. 
Once such a plane is fixed, the two points. Pi and are fully characterized 
by their transverse coordinates rj_i and rj_2 respectively. Our presentation will 
be given in the frequency domain. Due to the limited temporal resolution of 
detectors in Synchrotron Radiation experiments, the analysis in frequency do- 
main is much more natural than that in the time domain. As a consequence 
of this choice, and without restrictive assumptions on the system, we study 
the spatial correlation between Fourier transforms of the electric field ^ at a 
fixed frequency a;, that is the cross- spectral density 



In Eq. (1) E is the complex amplitude of the Fourier transform of a given 
Cartesian component of the electric field at the space-frequency point (2:0, fl, a;), 
the asterisk denotes complex conjugation, and brackets < ... > indicate an 
ensemble average over electron bunches. Since wc study an ultra-relativistic 
system, the paraxial approximation can always be enforced so that, here, the 
electric field is understood to obey the paraxial wave equation [1]. The cross- 
spectral density carries all information about the transverse characteristics of 
undulator radiation. A fully general study of the cross-spectral density is not a 
trivial one. Difficulties arise when one tries to include simultaneously the effect 
of intrinsic divergence of the radiation, due to the presence of the undulator, of 

^ Since our analysis deals with the electric field in frequency domain, we will some- 
times refer to the "Fourier transform of the electric field" simply as "the field", 
when this does not generate confusion. 




(1) 



6 



the electron beam size and of the electron beam divergence. In [2] a technique 
was described, based on Statistical Optics, to calculate the cross-spectral den- 
sity from undulator sources in the most general case ^ , at any position after 
the undulator but still without optical elements (i.e. in free space). Although 
self-contained, the present study relies on that work. Expressions from [2] will 
be taken as a starting point to proceed along the optical beamline towards 
the specimen position. 

In general, as we will see, undulator radiation can be thought as originating 
from an equivalent source localized on a transverse plane at a given longitudi- 
nal position. By definition, such source has the following property: it produces 
a field which coincides with that from the undulator at any distance from the 
exit of the undulator. In particular, we localize the equivalent source in the 
center of the undulator, that will be conventionally taken as the beamline 
origin z — 0. Since in this case the equivalent source does not reproduce the 
real electromagnetic field distribution in the center of the undulator, we refer 
to it as virtual source. Further on, throughout this paper we assume that the 
beta functions of the electron beam have their minimal value in the center of 
the undulator. By this, as we will see, the virtual source exhibits particular 
properties which simplify our treatment. Once the concept of virtual source 
is introduced, the problem of describing radiation characteristics at a certain 
observation plane after a given optical element is twofold. First, one has to 
characterize the cross-spectral density at the virtual source and, second, one 
has to propagate the cross-spectral density along the optical beamline to the 
observation plane. 

Let us first consider the problem of characterizing the source. An important 
simplified model which admits an analytical description without loss of essen- 
tial information about the source features is obtained by letting both horizon- 
tal and vertical electron beam emittances be much larger than the radiation 
wavelength {e^^y S> \/(2n)). This is a good assumption for second genera- 
tion Synchrotron Radiation sources. The kind of virtual source obtained for 
^x,y ^ A/ (27r) belongs to the wider class of quasi-homogeneous ones. These are 
characterized by the fact that the cross-spectral density at the virtual source 
plane (i.e. at z — 0) can be written as: 



G'(rj_i,f_L2 



Lu) ^ I {f±i,uj) g{f±2 - r±i,uj) , 



(2) 



where 




(3) 



^ With the point of view of the source parameters. 
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is the field intensity distribution and g{f±2 — is the spectral degree 

of coherence (normalized, by definition, so that g{0,uj) = 1). The definition 
of quasi-homogeneity amounts to a factorization of the cross-spectral density 
as the product of the field intensity distribution and the spectral degree of 
coherence, which contains information about the spatial correlation. A set of 
necessary and sufficient conditions for such factorization to be possible follows: 
(a) the radiation intensity at the virtual source varies very slowly with the 
position across the source on the scale of the field correlation length and (b) 
the spectral degree of coherence depends on the positions across the source 
only through the difference f±2 — 'r±i- 

There are situations when the Statistical Optics description is not the only 
one possible. The asymptotic limit for large electron beam emittances {e^^y S> 
X/{2n)) is one of these. In this limit, the Statistical Optics description of 
the source coincides with the Geometrical Optics (or Hamiltonian) descrip- 
tion of the source, where a photon-beam phase space is defined and can be 
described in terms of rays specified by position-angle coordinates. In the Ge- 
ometrical Optics approach, based on the uncertainty principle, only the ra- 
diation originating by a photon-beam phase space area of order [A/(27r)]^ is 
spatially coherent. When the emittance is much larger than the wavelength, 
27rea;^j^/A 3> 1, the divergence of the electron beam is much larger than the 
diffraction angle of undulator radiation, and the transverse size of the elec- 
tron beam is much larger than the diffraction size of undulator radiation. As 
a result one can completely neglect diffraction effects, and the Geometrical 
Optics approach can always be applied. Since Geometrical Optics describes a 
limiting situation of Statistical Optics there must be a relation between the 
fundamental Geometrical Optics quantity, the phase space distribution, and 
the fundamental Statistical Optics quantity, the cross-spectral density. It can 
be shown [3] that the radiant intensity of the field generated in free space by a 
quasi- homogeneous source in the direction of a unit vector s can be expressed 
as 



r{s,uj) being the two-dimensional spatial Fourier transform of the degree of 
transverse coherence g{f±2 — f±i,u!): 



The expression for the phase space distribution is given by the product of the 
intensity distribution of the source and the radiant intensity 




(4) 




(5) 



* (s, r±) = I {f±,Lj) X{s, uj) (xl (fj_, u;) V (s, a;) , 



(6) 
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where the variables {f±,s} characterize a ray in phase space. A comparison 
between Eq. (6) and Eq. (2) shows that cross-spectral density and phase space 
distribution contain the same information in the limiting case e^^y S> A/(27r). 

In other words, if a source has a large angular divergence (compared with the 
diffraction angle) and a large transverse size (compared with the diffraction 
size), one can completely neglect diffraction effects and treat the problem of 
the characterization of the source by means of Geometrical Optics ^ . 

Despite the previous discussion, the possibilities of using Geometrical Optics 
to describe undulator sources are quite limited in many realistic situations. 
Applications of Synchrotron Radiation make use of a very wide range of wave- 
lengths which span over four order of magnitude, from O.lA to lO^A. For third 
generation light sources, either planned or in operation, the horizontal electron 
beam emittance Sx — cTxCTx' is of order of 1 3 nm. The vertical emittance is 
given by ey = CyCTy' = x^x, X being the so called coupling factor. Typical values 
of X for third generation light sources are of order x ~ 0.01, corresponding 
to vertical emittances of order 0.1 O.sA. These values are always near or 
within the diffraction limit for wavelength ranges up to the hard X-rays in the 
vertical direction, and Geometrical Optics descriptions fail. In particular, in 
the VUV wavelength range, both vertical and horizontal emittances are much 
smaller than the radiation wavelength {e^^y <^ X/{2tt)). One recovers, then, 
the perfectly coherent situation when the source is diffraction limited in both 
horizontal and vertical directions. This is another situation when the Statisti- 
cal Optics description is not the only one possible. In this case, deterministic 
Wave Optics may be used as well. As the wavelength becomes shorter, in the 
soft X-ray range, one obtains <^ A/(27r), but ^ A/(27r). At wavelengths 
of about lA the vertical emittance reaches the same order of magnitude of 
the wavelength ey ~ A/(27r), while €x ^ A/(27r). Finally, in the hard X-ray 
region, at a wavelength of about O.lA, both emittances are much larger than 



^ Condition 2TT€x^y/X S> 1 is sufficient, but not necessary. In general, we cannot say 
that Geometrical Optics is never applicable for electron beam emittances smaller 
than the radiation wavelength A. We will treat this subject in a more extensive 
fashion in Section 7. There we will see that there are situations when Geometrical 
Optics can be applied to describe the (virtual) source even when the electron beam 
emittance is smaller than A. We will find that a sufficient (less restrictive, but 
still not necessary) condition for the applicability of Geometrical Optics to the 
description of a given undulator source is that such source can be characterized in 
terms of a quasi-homogeneous virtual source. Moreover, as it will also be discussed 
in Section 7, our comparison of the emittance with the radiation wavelength is done 
under the assumption that the electron beam beta function is comparable with the 
radiation formation length at wavelength A. This is often, but not always, the case 
for undulator sources, since the radiation formation length is the undulator length, 
which is at least a few meters. However, it is not the case for bending magnet 
radiation. 
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the wavelength {e^^y ^ A/(27r)), and Geometrical Optics can be used along- 
side Statistical Optics. It follows that, for third generation light sources, only 
the limiting cases for wavelengths around 100 nm and 0.1 A can be treated, re- 
spectively, my means of Wave Optics or Geometrical Optics. The intermediate 
situation can be treated in a rigorous way only with the help of Statistical Op- 
tics, which includes both Wave Optics and Geometrical Optics as asymptotic 
cases. 

Strictly related to the problem of source characterization, but separate from 
it, is the issue of propagating the photon beam through the optical beamline 

to the observation plane. In the case of quasi-homogeneous virtual sources, 
if diffraction effects from the optical elements can be neglected. Geometrical 
Optics can be taken advantage of. The virtual source can be described in 
terms of phase space distribution, and interactions with optical media can be 
conveniently modelled in terms of symplectic transformations, very much like 
electron beams in storage rings optics. Several computer codes (e.g. SHADOW 
[4]), usually referred to as ray-tracing codes, have been developed and are 
standard tools used to carry out Geometrical Optics-based calculations. How- 
ever, this approach is not always possible as the virtual source may not be 
quasi-homogeneous or diffraction effects may not be neglected in the optical 
beamline. A rigorous analysis of the object-image coherence relationship is of 
fundamental importance in the context of several coherence-based techniques 
like fluctuation correlation dynamics, phase imaging, coherent X-ray diffrac- 
tion and X-ray holography, whose development has been fostered by the high 
flux of coherent X-rays provided by state-of-the-art third generation facihties. 
It should be noted that, in the case of partially coherent wavefronts, even 
the calculation of the intensity distribution at the specimen position should 
involve Statistical Optics techniques. In fact, to obtain the intensity at the 
specimen position as some optical element is present one first needs to track 
the cross-spectral density through the beamline i.e. one needs to study the 
evolution of the partially coherent wavefront. 

Computer codes have been written [5] in order to deal with beamline design 
in the case of partially coherent radiation. These are devoted to the solu- 
tion of the image formation problem starting from first principles. Results 
may in fact be obtained using numerical techniques alone, starting from the 
Lienard-Wiechert expressions for the electromagnetic field and applying the 
definition of the field correlation function without any analytical manipula- 
tion. Yet, a first-principle calculation of the field correlation function between 
two generic points or, in particular, calculation of the intensity at a single 
point involves very complicated and time-expensive numerical evaluations. To 
be specific, one needs to perform two integrations along the undulator device 
and four integrations over the electron-beam phase space distribution to solve 
the problem in free space. Then, modelling the optical beamline as a single 
convergent lens, other four integrations are needed to characterize coherence 
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properties on the image plane, for a total of ten integrations. The development 
of a universal code for any experimental setup is then likely to be problematic. 
A more conservative approach may suggest the use of computer codes based 
on some analytical manipulation of first principle equations suited for spe- 
cific experimental setups. Prom this viewpoint our most general expressions 
may be used as reliable basis for the development of numerical methods. Yet, 
computer codes can calculate properties for a given set of parameters, but 
can hardly improve physical understanding, which is particularly important 
in the stage of planning experiments. Our theory will allow treatment and 
physical understanding of many asymptotes of the parameter space and their 
applicability region with the help of a consistent use of dimensional analysis. 
In the most general asymptotic cases treated here, this will allow to reduce 
first principle calculations (i.e. ten- dimensional integrals) to one-dimensional 
convolutions of analytical functions with universal functions specific for the 
undulator source case, and still to retain a certain degree of generality. It is 
also worth to underline that our asymptotic results may also be used as a 
benchmark for numerical methods. 

One of the main difficulties in applying a Statistical Optics approach to Syn- 
chrotron Radiation sources stems from the fact that Statistical Optics has 
principally developed in connection with problems involving thermal light. 
Solutions to all these problems share approximations that allow major sim- 
plifications, but are specific of thermal sources only. For instance, thermal 
sources can be modelled as perfectly incoherent, and the cross-spectral den- 
sity assumes the form 

G(r±i,fj_2,a;) oc / (fxi,a;)(5 (f_L2 - f_Li) , (7) 



where / is the source intensity distribution and 5 is the two-dimensional Dirac 
(5-function. However, there is a close connection between the state of coher- 
ence of the source and the angular distribution of the radiant intensity (see 
Eq. (4)). The physical interpretation of Eq. (7) is that the source is correlated 
over the minimal possible distance (which is of order of the wavelength). This 
has the consequence that the radiant intensity is distributed over a solid angle 
of order 27i. This is correct for thermal sources, but is in contradiction with 
the fact that any Synchrotron Radiation source is confined within a narrow 
cone in the forward direction. The high directionality of Synchrotron Radia- 
tion rules out the use of Eq. (7) as a model for Synchrotron Radiation sources. 
However, such high directionality is not in contrast with the poor coherence 
which characterizes the quasi-homogeneous limit. Quasi-homogeneous sources 
are only locally coherent over a distance of many wavelengths but, by defini- 
tion of quasi-homogeneity, the linear dimension of the source is much larger 
than the correlation distance. Even though a quasi-homogeneous source can be 
described with Geometrical Optics techniques, a coherence distance of many 
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wavelengths rules out the use of Eq. (7) as a model for Synchrotron Radia- 
tion sources. A more precise knowledge of the cross-spectral density (that is 
equivalent to the knowledge of the correct phase space density) is necessary to 
solve the image formation problem. For instance, suppose that a light source 
is placed at arbitrary distance in front of a lens. If the source is perfectly 
incoherent (thermal light case) the area of the light incident on the lens is al- 
ways the area of the lens. In the case of Synchrotron Radiation source though, 
such area may be smaller than the lens. Even in the limit for a large beam 
emittance (compared with the radiation wavelength), information about the 
small angular distribution must be printed in the wavefront at the exit of the 
undulator leading once more to the same conclusion: Eq. (7) cannot be used 
in order to model Synchrotron Radiation sources. In [2] we treated, among 
other cases, the asymptote for a large electron beam size and divergence. The 
expression for the cross-spectral density of the source in free space simplifies 
and a particular quasi-homogeneous model can be given. In the same work, 
we specified also the region of applicability of such model, and we showed 
that it cannot be applied outside the limit for a large electron beam size and 
divergence. 

In relation with these remarks it should be mentioned that an attempt to fol- 
low the path proposed in this paper is described in [6]. To our knowledge, [6] 
constitutes the first remarkable attempt to use Statistical Optics techniques in 
order to characterize the evolution of partially coherent X-ray beams through 
optical systems. In that work, as well as in [7], the beamline optics from the 
undulator to the specimen can be modelled as a critical illumination system 
[8], the beamline behaving as the condenser. After this. Statistical Optics 
techniques are consistently used to calculate coherent properties on the im- 
age plane. However, the authors of [6] reduced the general ten-dimensional 
integrals to four-dimensional integrals by postulating that the cross-spectral 
density distribution at the exit of the undulator can be written as Eq. (7), i.e. 
a perfectly incoherent source is assumed at the exit of the undulator. As we 
have just seen though, this assumption is always inconsistent in the case of 
Synchrotron Radiation, even in the Geometrical Optics limit and, a fortiori, 
in the case treated by the authors (the undulator beamline 12 at ALS), where 
€y ~ 0.1A/(27r) and ~ 3A/(27r), which is highly spatially coherent. 

We organize our work as follows. Besides this Introduction, in Section 2 we 
describe the optical system under study and some concepts from Statistical 
Optics that will be widely used in the following Sections. In Section 3 we 
review some general expressions pertaining undulator radiation from a single 
particle. In particular, following [2] we present an analytical expression for the 
Fourier transform of the electric field generated by a single electron with offset 
and defiection which is valid at any distance from the exit of the undulator. We 
also present the analytical solution of the imaging problem for a deterministic 
model of undulator radiation (absence of electron beam emittance). In Section 
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4 we give a derivation of the cross-spectral density for undulator radiation 
based sources. Subsequently we analyze the evolution of the cross-spectral 
density function through the optical system with particular attention to the 
focal and to the image plane. The following two Sections 5 and 6 describe 
quasi-homogeneous sources, respectively Gaussian and non-Gaussian, in the 
ideal case when the lens is aberration-free and the pupil aperture is non- 
limiting. A digression is then taken in Section 7, where we analyze in detail 
the relation between Geometrical Optics and quasi-homogeneous sources. Such 
Section may therefore be skipped in a first reading, without interrupting the 
main logical stream of our work. The next Section 8 describes the effects of 
a finite aperture size on the radiation characteristics from quasi-homogeneous 
sources at the image plane, and is followed by Section 9 that assumes a quasi- 
homogeneous source as well and deals with the consequences of lens aberration 
on the intensity at the image plane. In Section 10 we introduce a particular 
setup, a pinhole camera, capable of producing images of the sources in the 
absence of lenses. The study of this particular setup is of particular interest 
because it allows the reader to recognize mathematical analogies between cases 
otherwise physically different and serves as a junction between the previous 
Section 9 and Section 11 treating physical cases when one obtains, surprisingly, 
the image of the source on the focal plane. In the following Section 12 we 
extend the treatment for the focal plane to any plane of interest. In Section 
13 we discuss the depth of focus, including the case of a large non-limiting 
aperture and the effects of aperture size. In the next Section 14 solutions 
for the image formation problem in non-homogeneous cases relevant for third 
generation Synchrotron Radiation sources are given. Before conclusions, in 
Section 15, we discuss the accuracy of quasi- homogenous source asymptotes. 
Finally, in Section 16, we come to conclusions. 



2 Elements and definitions of image formation theory 

2. 1 Wave propagation in free space 

Let us indicate with E±{z, f±, t) any fixed polarization component (along the 
direction x or y) of the electric field at time t calculated on a plane at position 
z down the beamline at a certain transverse location f±. E±{z,f±,t) obeys, 
in free space, the homogeneous wave equation: 

^^V^^. - ^ = . (8) 
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Let us now introduce the Fourier transform E±{z,'r±,u!) of the electric field 

E±{z,r±,t): 

oo 

E±{cu) = J dtE^{t)e'''^ , (9) 



so that 



1 °° 

(i) - ^ / dujE^{uj)e-"^ . (10) 



As already remarked in footnote 1, we will sometimes refer to E_[_ as "the 
field" , understanding that we are working in the frequency domain. 

Let us consider the field propagation problem. To this purpose, we first intro- 
duce the complex envelope of the field: 

E = E^ex.^[-iujz/c] . (11) 



It is always possible to give such a definition. However, its utility is restricted 
to the case when is a slowly varying function of z with respect to the 
radiation wavelength A. When the paraxial approximation is applicable (i.e. 
always, for Synchrotron Radiation sources), this condition is fulfilled [1]. 

In paraxial approximation and in free space, the following parabolic equation 
holds for the complex envelope E of the Fourier transform of the electric field 
along a fixed polarization component: 



The derivatives in the Laplacian operator V^^ are taken with respect to the 
transverse coordinates. One has to solve Eq. (12) with a given initial condition 
at which is a Cauchy problem. Indicating with To the transverse coordinate 
of an observation point on a plane at longitudinal position Zo we have 



E{Zo, fo) 



lUO 



2'!rc{zo — z) 



j dr' E{z, r') exp 



21 



Ta — r 



2c{zo - z) 



(13) 



where the integral is performed over the transverse plane. 
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Next to the propagation equation for the field in free space, Eq. (13), we can 
discuss a propagation equation for the spatial Fourier transform of the field, 
which can also be derived from Eq. (12) and will be useful in the following 
parts of this work. We will indicate the spatial Fourier transform of E{z,r') 
with F(z,'ir)^: 



{z, — j dr'E{z, r') exp 



ir' ■ u 



(14) 



Eq. (12) can be rewritten in terms of F as 

|y du¥ {z,u) exp [—ir • ■u| | = . 



„ 9 2iuj d 
c dz , 



(15) 



Eq. (15) requires that 



-\u\ + 



2 , 2iuj d 



(16) 



Solution of Eq. (16) can be presented as 



F {z, u) = F (0, u) exp 



IC\U\ z 



2a; 



(17) 



It should be noted that the definition of F(0, u) is a matter of initial conditions. 
In many practical cases, including the totality of the situation treated in this 
paper, F (or E) may have no direct physical meaning at 2; = 0. For instance, in 



^ For the sake of completeness we explicitly write the definitions of the two- 
dimensional Fourier transform and inverse transform of a function g{f) in agreement 
with the notations used in this paper. The Fourier transform and inverse transform 
pair reads: 



I 



dr g{r) exp 



ir ■ k 



; 9{r) = ^ I <^k g{k) exp 



—ir ■ k 



the integration being understood over the entire plane. If g is circular symmetric 
we can introduce the Fourier-Bessel transform and inverse transform pair: 

g{k) = 2Tr J dr rg{r)Jo{kr) ; g{r) = J dk kg{k)Jo{kr) , 



r and k indicating the modulus of the vectors r and k respectively, and Jo being 
the zero-th order Bessel function of the first kind. 
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all cases considered in this paper, ;i; = is in the center of the undulator, well 
within the radiation formation length. However, F(0, u) can be considered as 
the spatial Fourier transform of the field produced by a virtual source. Such 
a source is defined by the fact that, supposedly placed at z — 0, it would 
produce, at any distance from the undulator, the same field as the real source 
does. The result in Eq. (17) is very general. On the one hand, the spatial 
Fourier transform of the electric field exhibits an almost trivial behavior in 
z, since |F(^)p = const. On the other hand, the behavior of the electric field 
itself is not trivial at all (see Section 2.4 and Fig.s 3-6). These properties follow 
directly from the propagation equation for the field and its Fourier transform. 

Let us now discuss the physical meaning of Eq. (17). The spatial Fourier 
transform of the field, F{z, u), may be interpreted as a superposition of plane 
waves (the so-called angular spectrum) . Once the frequency uj is fixed, the wave 
number k — cu/c is fixed as well, and a given value of the transverse component 

— * 

of the wave vector k± — u corresponds to a given angle of propagation of a 
plane wave. Different propagation directions correspond to different distances 
travelled to get to a certain observation point. Therefore, they also correspond 
to different phase shifts, which depend on the position along the z axis (see, 
for example, reference [9] Section 3.7). Free space basically acts as a Fourier 
transformation. This means that the field in the far zone is, phase factor and 
proportionality factor aside, the spatial Fourier transform of the field at any 
position z. To show this fact, we first recall that if we know the field at a given 
position {z,r') we may use Eq. (13) to calculate the field at another position 
{zo,K) ■ Let us now consider the limit Zg — > oo, with finite ratio fo/ ^o- In 
this case, the exponential function in Eq. (13) can be expanded giving 



E{Zo, fo) 



I'KCZr, 



dr' E{z, r') exp 



i6i;|fo| ioj{ro ■ '"O iujz\fo 



^ 12" 



2cZn 



CZn 



+ 



2czl 



(18) 



— * 

Letting 9 — ro/^o we have 



2'KCZ, 



cxp 



2c 



{Zo + Z) 



(19) 



With the help of Eq. (17), Eq. (19) may be presented as 



E{Zo,fo) 



lUJ 



27rcZo 



exp 



2c 



^0 

F|0,-- 



(20) 



Eq. (20) shows what we wanted to demonstrate: free space basically acts as a 
Fourier transformation. 
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Fig. 1. Single lens imaging system with an undulator source as object. 
2.2 Image formation with coherent light 



As has been remarked in [6], any beamhne optics used to re-image undulator 
radiation to an observation plane of interest can be modelled as a critical 
illumination system [8]. Therefore, the setup considered in this paper can be 
sketched as in Fig. 1. It consists of an undulator of length centered at 
z = 0, a. convergent lens positioned ai z = Zi, characterized by a focusing 
strength /, and a plane of observation at position 2; = 2;2. In principle, zi, zi 
and / are unrelated parameters. However our main case of interest is a critical 
illumination system. Therefore, a given source plane at coordinate z = Zg - 
the object - is imaged at a particular observation position z = Z2, that defines 
the image plane. Using Ray Optics we can calculate the distance Z2 along the 
axis behind the lens where the image is formed. This gives the well-known 
lens-maker equation 



1 

7 



+ 



Z2 — Zl 



(21) 



The size of the image is magnified by a factor |M| and a real image is inverted 
with respect to the object because 



Z2 — Zl 
Zl — Zg 



<0, 



(22) 



as expected from Ray Optics. We also define a scale factor m, that is the 
inverse of the magnification power |M| of the lens: 



m 



1 

W\ 



Zl — Zs 
Z2 - Zl 



(23) 
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Once the lens position is fixed ai z = zi, the position Zg of the source is a 
matter of choice. Such a choice fixes the position Zi = z^ of the image plane 
in agreement with Eq. (21). For instance, with [6] and [7], one may set the 
critical illumination system to image the radiation at the undulator exit, in 
which case the choice Zg — Ly,/2 is made. However, one is not obliged to do so. 
In particular, in this paper, we will make the choice Zs = 0. Then, the critical 
illumination system images the center of the undulator, in the sense that the 
image plane Zi = z^ obeys Eq. (21) with Zg = 0. 

From a mathematical viewpoint, specifying the source is equivalent to fixing 
the initial conditions for Maxwell equations in terms of a field distribution on 
a certain transverse plane. When a certain field distribution (in the frequency 
domain) is fixed on a plane at position z, Maxwell equations automatically set 
the way radiation propagates in free space, and the source is univocally defined. 
The choice made in [6] and [7] corresponds to the choice of a real source. 
The denomination "real" is justified by the fact that the initial condition for 
Maxwell equations amounts to the specification of a field distribution which is 
actually present, and in principle measurable, at the exit of the undulator. On 
the contrary, our choice Zg — Q corresponds to the position down the z-axis 
in the middle of the undulator, well within the radiation formation length. 
Although it makes sense to talk about the distribution of the electromagnetic 
field in the middle of the undulator, it does not make any sense to identify 
such distribution with the initial condition for Maxwell equations, i.e. with 
the source. However, it makes sense to define a virtual source as in Section 
2.1. There is a particular reason for the choice z = as the position for the 
virtual source. In the case of a coherent undulator source (that is being treated 
in the present Section), the wavefront of the radiation at a virtual source 
located at 2; = is plane. In other words, the far zone field from an undulator 
has spherical wavefronts centered in the middle of the undulator. This fact 
alone makes the center of the undulator a privileged point with respect to 
others. Moreover, in the more general case of partially coherent radiation we 
will assume (as it is often verified in practice) that the beta functions of the 
electron beam have their minima (in both horizontal and vertical directions) 
in the center of the undulator. Subject to this assumption, as we will see, the 
center of the undulator constitutes a privileged point of interest in this case 
as well. As a result of the previous discussion we set 

Zg^Q . (24) 

Let us assume that the position zi and the focal length / of the lens are 
set. Throughout this paper we will be particularly interested in the radiation 
characteristics at two privileged positions down the beamline: 

• the image plane, at position Z2 identified by Eq. (21). 

• the focal plane, at position Z2 identified by the equation Z2 — Zi — f. 
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As it will be seen in the next Sections, both image and focal planes have 
special properties that can be expressed in terms of Fourier Optics. These 
properties are valid for any wavefront. We will first take advantage of them in 
Section 3, where we will deal with wavefronts generated by an electron beam 
with zero emittance and further on in Section 4, where emittance effects will 
be discussed in the realm of Statistical Optics. In the present Section, after 
having set the configuration under study, we will limit ourselves to describe 
these properties and to define basic quantities to be used in the Statistical 
Optics formulation of the image formation problem. 

Consider the problem of mapping the plane immediately in front of the lens 
onto the plane at longitudinal position Z2 > Zi behind the lens (see Fig. 1). If 
one knows E at position z = zi immediately in front of the lens, one can also 
obtain the expression for E immediately behind the lens multiplying by the 
transmission function: 



T (r'^ = P {r'^ exp 



—llxJ 
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2c/ 



(25) 



For simplicity of notation we consider here identical focal distances in the 
horizontal and in the vertical direction, i.e. f = fx = fy Results for more 
complicated optical systems (e.g. a combination of cylindrical mirrors) are 
found substituting consistently Eq. (25) with its straightforward generaliza- 
tion. We assume with [8] that the complex pupil function P is zero outside the 
lens aperture. Its phase accounts for aberrations and its modulus may vary 
along the lens to describe apodizations: the simplest possible study case is for 
\P\ = 1 and arg (P) = within the lens aperture. Accounting for Eq. (25), the 
propagation equation for any field E{zi,r') immediately in front of the lens to 
the point {z2, r2) on the observation plane behind the lens can be written as: 



E{z2, fa) 



2'KC{Z2 - Zi) 

iuj ( 



exp 



iiO\r2\ 



X exp 



2c(^2 - ^1) 

1 ^ 



c yi{z2-zi) 2/^ 



dr' \ E{zx,r')P{r') 
iijj{f2 ■ r') 

C {Z2 - ^1) 



exp 



.(26) 



Since Synchrotron Radiation is highly collimated, it is practically relevant 
to discuss the case when the area of the spot of the incident radiation is 
small compared with the area of the lens. It is also simpler and more natural 
to start with this situation. Then, effects from a finite pupil dimension can 
be neglected. Note that this is not the case for thermal sources: since these 
are emitting into a solid angle of I-r (they are perfectly incoherent) the finite 
pupil dimensions cannot be ignored and result in the so-called vignetting effect 
[8]. Considering a perfect lens with no aberrations too, Eq. (26) assumes a 
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particularly simple form at the focal and on the image plane. Initially we will 
consider situations when the pupil presence can be neglected. Later on we 
discuss how to include the effects due to the presence of the pupil. 



2.2.1 Large non-limiting aperture 

We will now specialize the result in Eq. (26) in the asymptote for a large 
non-limiting aperture and in the case of the focal and of the image plane. Let 
us denote with {zf, r/) a point on the focal plane, and with {zi, r^) a point on 
the image plane. Prom Eq. (26), on the focal plane we have 



tio\rf\ 



j d? 



E{zi,r') exp 



— * 

iu!{ff ■ r') 



(27) 



With the help of Eq. (14) we can write Eq. (27) as 



27rcf 



exp 



F Ui, 



cf . 



Substitution of Eq. (17) in Eq. (28) gives 







2' 


exp 


iu}Zi\rf 


2' 


[ 2c/ J 


[ '2cP J 



F 0, 



cf . 



(29) 



For the image plane, remembering that 

1 1 1 

= — + 



/ Zi Zi — Zi 



(30) 



and that 



di — Zi — Zi — 



m 



(31) 



we obtain, from Eq. (26) 



E{zi,fi 



271CZ] 



exp 



iujm 


Ti 


2" 


/*' { 


2czi 









exp 


2czi 


jexp 



ia;m(rj ■ r') 



CZi 



(32) 
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On the image plane, according to Eq. (32), we have to calculate the Fourier 
transform of the product of two factors: exp [—iuj\r'\'^/ {2czi)\, representing the 
phase of a spherical wave in paraxial approximation, and E{zi,r'). A direct 
calculation shows that the Fourier transform of the phase factor is 



/ 



dr' exp 



ioj\r 



/|2 



2czi 



exp 



iojv[i{ri ■ r') 

CZi 



—Aizi exp 



2 1^12 



2czi 



(33) 



Since the Fourier transform of a product is equal to the convolution of the 
Fourier transforms of each separate factor, from Eq. (32) one obtains 



E{zi,fi) 



m 



47r2 



exp 



ia;m|nP 



2czi 



X j duF (zi, u) exp 



ICZi I TClUTi 



2u 



u 



CZi 



(34) 



Substitution of Eq. (17) in Eq. (34) gives 



E{zi,fi) 



m 

4^ 



exp 



iuim 


n 


2' 


exp 






2" 


2czi 




2cz\ 



X j du F(0, u) exp [imfi ■ u] , 



(35) 



that is 



E{zi,fi)^mexp 



iujm 




2' 


exp 






2" 


2cz\ 


2cz\ 



E (0, — mrj) 



(36) 



The phase factor in the Fourier transform of the electric field, given in Eq. 
(17), cancels the quadratic phase factor in \u\^ in Eq. (34). Therefore, the 
convolution integral in Eq. (34) transforms to a Fourier integral. As a result, 
in the image plane we always obtain (aside for a scaling and a net phase factor) 
the inverted field distribution in the virtual source plane. More in general, at 
any observation plane located aX z = z^ behind the lens and the focal plane, 
one observes (aside, again, for a scaling and a net phase factor) the inverted 
field distribution on an object plane located z — Zs, where Zg satisfies the 
lens condition Eq. (21). 

Eq. (29) and Eq. (36) are reflections of well-known theorems of Fourier Optics. 
Neglecting the effects from a finite pupil dimension and assuming a perfect 
lens with no aberrations, the focal plane has the following property [9]: 
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• For any position of the object in front of the lens, the field distribution 
(amplitude and phase) on the focal plane differs from the spatial Fourier 
transform of the field distribution on the object plane by a scale factor 
—u)/{cf) and a net phase factor. 

At the image plane, instead, the following property applies [9] : 

• For any position of the object in front of the lens, the field distribution 
(amplitude and phase) on the image plane differs from the field distribution 
on the object plane by a scale factor — m and a net phase factor ^ . 

These properties can be interpreted in terms of intensity distributions. The 
first tells that the intensity profile on the focal plane has the same shape of 
that on a distant plane and is obtained taking, essentially, the square modulus 
of the FoTiricr transform of the field on the object plane. The second tells that 
the intensity profile of the object is inverted and magnified by the lens on 
the image plane. For the image plane we just obtained, for perfectly coherent 
light, the same result which is obtained in Geometrical Optics in the perfectly 
incoherent limit. 



2.2.2 Effect of aperture size 

Pupil effects are taken into account, from a general standpoint, in Eq. (26). 
Eq. (26) takes a specific form in the focal and in the image plane. One has 



E{zi.rf) 



2^ 



exp 



iuj\rf\ 



j dr' E{zi,r')P{r') exp 



iu!{ff ■ r') 



(37) 



and 



~ ^ turn 
E{Zi,ri) = - exp 



27rczi 
X exp 



icum 




2" 


jiP { 


2czi 



dr' < E[zi, r')P{r') exp 



iui\r 



7|2 



2czi 



iujm.[ri ■ r 

CZi 



(38) 



Eq. (37) and Eq. (38) are formal extensions of Eq. (27) and Eq. (32), respec- 
tively. Use of the convolution theorem on Eq. (37) and Eq. (38) allows to write 
analogous extensions of Eq. (29) and Eq. (36). To this purpose we define 



^ The prefactor ioa{zi) is a consequence of the conservation of the total energy 
associated with the propagating field. 
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V{u) = j dr'P{r') exp 



—ir' ■ u 



(39) 



Then, indicating with Ep the field in the presence of the pupil, on the focal 
plane we have 



2c/ 



/ 



X duV 



u ■ exp 



icf\u\ 
2uj 



cfu 



UJ 



(40) 



and on the image plane 



X 



duV 



2czi 

( ujvafi 



\ CZi 



u ■ exp 



iczi\u\^ 
2a;m 



E\zi, 



CZiU 

com , 



• (41) 



One may apply the following mnemonic rule to include the effects of the pupil 
in Eq. (29) or Eq. (36). First, divide Eq. (29) or Eq. (36) by the first phase 
factor, corresponding to the phase factor outside the integral sign in Eq. (26). 
Second, convolve with V. Third, put the phase factor back. 

It should be noted that, in the limit for large apertures, V can be substituted 
by a (5-Dirac function in both Eq. (40) and Eq. (41). In this case, from Eq. (40) 
we recover Ep{zf,ff) = E(zf,ff), given in Eq. (29). Form Eq. (41) instead, 
we have Ep{zi,ri) — E{zi,'ri), given in Eq. (36). 

Unless particular conditions are met, the phase factors under the integral signs 
in Eq. (40) and Eq. (41) compensate only partially the phase of E, which can 
be found in Eq. (29) and Eq. (36) respectively. This fact complicates the 
evaluation of the convolution integrals. Let us restrict our attention to the 
image plane. We can treat analytically the case when the phase factors under 
integral in Eq. (41) completely compensate the phase factor in Eq. (36), i.e. 
when we can neglect the second phase factor in Eq. (36). This happens when 
we are in the far field limit. The far field limit of Eq. (36) can be obtained 
by substitution of Eq. (20) in Eq. (32). After the inverse Fourier transform in 
Eq. (32) is calculated one obtains 



E{zi,fi) =mexp 



2cz\ 



E (0, -mfi) , 



(42) 



that can also be obtained directly from Eq. (36) neglecting the second phase 
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factor on the right hand side. This is possible when 



2 



< 1 



(43) 



2czi 



for any point fi on the image pattern. 

Let us indicate with the characteristic size of the image. Condition (43) can 
be interpreted as the following requirement for zi: 



Since we are interested in the parametric dependence only, a factor 2 has been 
neglected in Eq. (44). From Eq. (36) follows that mo"j is the characteristic size 
of the virtual source. As such it is independent of the position of the lens and 
of the magnification factor \M\ = as well. Condition (44) is often met 
in practice and means that the radiation spot size on the lens, czi/ (o^mcxi), is 
much larger than the characteristic size of the virtual source, mcTj. Then, the 
lens is placed in the far zone with respect to the virtual undulator source by 
definition of far zone \zi/ (27r(To) ^ do, (^o being the source size. We conclude 
that the condition for the lens to be placed in the far zone is equivalent to 
the condition that Eq. (36) can be reduced to Eq. (42). This result will be of 
importance in what follows. However, this kind of reasoning is only valid on 
the image plane. On the focal plane the second phase factor in Eq. (29) follows 
directly from the phase factor in Eq. (17), that is related with the propagation 
of the angular spectrum: in this case one concludes that plane waves with 
different directions of propagation lead to an increasing phase difference as 
the distance zi increases. Therefore, in the focal plane, some simplification 
may be obtained in the near field only, as zi is small enough that the phase 
difference between different plane wave components is negligible. We will not 
investigate this situation further. Coing back to the image plane one can see 
that a term of the expansion of the phase factor under the integral sign in 
Eq. (34) cancels the phase factor in Eq. (17). It follows that the second phase 
factor in Eq. (36) is not related with the propagation of the angular spectrum. 
Therefore, in the far field region, when condition (44) holds, we have that 
the pupil effects can be accounted for by means of a simpler convolution. In 
the following parts of this paper we will restrict to this particular case when 
treating pupil effects. It is interesting to remark that pupil effects due to finite 
pupil dimension are especially important in the far region. In this limit the 
radiation spot size on the lens is much larger than the size of the radiation 
spot at the virtual source, and is often larger that the size of the pupil. In the 
near field instead, the radiation spot size on the lens is of order of the radiation 
spot size at the virtual source. Therefore, in this hmit, one can neglect effects 



Zi > 



(44) 



c 
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from any finite pupil aperture larger than the radiation spot size at the virtual 
source. 

In conclusion, explicit substitution of Eq. (42) in Eq. (41) yields the following 
far field limit expression on the image plane: 

2.3 Propagation of partially coherent light in free space 

Let us now consider the stochastic nature of the Synchrotron Radiation field 
in general terms. Synchrotron Radiation is a Gaussian stochastic process. As 
has been discussed in detail in reference [2], the electromagnetic signal at 
any position down the beamline is completely characterized, from a statistical 
viewpoint, by the knowledge of the second order field correlation function in 
space-frequency domain 

r^{Zo,fol,ro2,U},u') ^ (^E±{Zo,fol,^)E*j_{Zo,fo2,^^')) ■ (46) 

In this paper, the averaging brackets (...) will always indicate an ensemble 
average over bunches. As it will be better explained in the following Section 
3, we will restrict ourselves to the treatment of radiation from planar undu- 
lators in resonance with the fundamental harmonic. Therefore we can neglect 
vertically polarized radiation components and consider E as a scalar quantity. 
The shot noise in the electron beam is responsible for random fiuctuations of 
the beam density, both in space and time. As a result, the temporal Fourier 
transform of the Synchrotron Radiation pulse at a fixed frequency and a fixed 
point in space is a sum of a great many independent contributions: 

E±{Zo,fo,Uj) = J2^s±iffk,lk,Zo,fo,(^)GXp{lUjtk) , (47) 

fe=l 

— * 

where N is the number of electrons in the bunch. Here ffk, h and tk are random 
variables describing random angular direction, position and arrival time of an 
electron at the reference position z^ = 0. As has been demonstrated in [2], 
under the assumption - generally verified for X-ray beams and third generation 
hght sources - that the radiation wavelengths of interest is much shorter than 
the bunch length we can write Eq. (46) as 

ru,{zo,foi,fo2,uj,uj')^NF^{uj - Uj') 



Ep{zi,ri)=exp 



zum\ri 
2czi 
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x(Es±{f),l,Zo,foi,i^)E*±{'nJ,Zo,ro2,i^')) , (48) 

\ / rfX 

where F{lj) is the Fourier transform of the bunch longitudinal profile function 
Ft{tk), that is 



(exp {iixitk))t 



(49) 



Note that the ensemble average on the right hand side of Eq. (48) is done over 
the product of the electric field produced by the same electron. In other words 
each electron is correlated only with itself. 

If the dependence of Es± on u and uj' is slow enough, so that Es± does not vary 
appreciably on the characteristic scale of F^ we can substitute E*j_{'r], l,Zo,fo2, 

— * 

with E*j_(f],l, Zo,fo2,i^) in Eq. (48) thus obtaining: 



^oj{zo, roi, ro2, UJ, u') = NF^{uj - uj')G{zo, foi, fo2, uj) 
where 



(50) 



Gizo, foi, fo2, (^) = ( Es±{fi, I, Zo, foi,uj)E*j^{ff, I, Zo, ro2, uj) ) ■ (51) 
\ / ffX 

As has been shown in [2] this assumption is by no means a restrictive one. 

From now on we will be concerned with the calculation of the correlation func- 
tion G{zo,roi,^o2,uj), while the correlation in frequency (which may be com- 
plicated by other factors describing, for instance, the presence of a monochro- 
mator) can be dealt with separately. 

On the one hand, the cross-spectral density as is defined in Eq. (51) includes 
the product of fields which obey the free space propagation relation Eq. (13). 
On the other hand, the averaging over random variables commutes with all 
operations involved in the calculation of the field propagation. More explicitly, 
introducing the notation E{zo) = 0[E{z)] as a shortcut for Eq. (13) one can 
write 



G{zo) = {Es^iz,)E:^{z,] 



O 



E.^Az] 



O* 



E*(z] 



O ■ O* 



e,^{z)e:^{z] 



O ■ O* [G{z) 



(52) 



Note that O may represent, more in general, any linear operator. 



26 



As a result, one can obtain a law for the propagation of the cross-spectral 
density in free space in analogy with Eq. (13) from position z to position Zq- 



G{Zo,fol,fo2) 



UJ 



A7r^c^{zo - zY 



X exp 



2c{zo - z) 



dr'idr'2 G{z,r'i,r'2) 



roi - r'l - fo2 - r'2 



(53) 



where the integral is performed in four dimensions. We may now proceed in 
parallel with Section 2.1. In analogy with Eq. (14) let us first define 



G {z, Ml, U2) = j dr'idr'2 G (2;, r'l, r'2) exp i (ui ■ r'l - U2 ■ r'2) . (54) 

Eq. (54) is a second-order correlation function between spatial Fourier trans- 
forms of the field. In the following, for simplicity, and with some abuse of lan- 
guage, we will denote Q as the "Fourier transform of G" . The spatial Fourier 
transform of E depends on the position along the beamline through a phase 
factor only. Moreover, as already said, the operation of ensemble average com- 
mutes with the operation of Fourier transform. It follows that also G depends 
on the position along the beamline through a phase factor only. To be specific, 
the analogous of Eq. (17) is given by 



g {z, ui, U2) = g (0, ui, U2) exp 



ic (.^ ,2 



F2 



(55) 



Continuing in analogy with Section 2.1 we find Eq. (20), which relates the 

far field expression for E{zo^ To) (in the limit Zo — > 00 and for a finite ratio 
fo/zo) to the spatial Fourier transform F. We can take advantage of Eq. (20) 
to obtain, with the help of Eq. (51) and Eq. (54), a useful relation between 
the cross-spectral density in the far field and the Fourier transform of G at the 
virtual-source position. In the limit Zg — > 00 and for finite ratios 9i — foi/^o, 

— * 

^2 = T02I Zo we have 



G(-o,r„i,r.2) = ^^exp[^(^^i - Q2 )\^\^^-—-—) " 

(56) 

This expression will be very useful later on. In [2] we obtained an explicit 
expression for the cross-spectral density of the undulator source in free space 
at any distance from the undulator. In particular we can calculate the cross- 
spectral density in the far field, which assumes a simplified form. Consequently, 
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the use of Eq. (56) allows to calculate the Fourier transform of the cross- 
spectral density at the virtual-source position. As a result we can characterize 
the virtual source and operate with it. Our starting point, here as in [2] is 
the electron beam in the undulator device. In contrast with this, previous 
literature dealing with application of Statistical Optics to undulator radiation 
assumes a priori the validity of a postulated expression for the cross-spectral 
density. 



2.4 Image formation with partially coherent light 



In analogy with Section 2.2 we will now consider the problem of propagating 
the cross-spectral density immediately in front of the lens through the optical 
system up to the screen at z = Z2 behind the lens. Suppose that we know 
G{zi) immediately in front of the lens. The cross spectral density immediately 
behind the lens, Gi{z-i), is related to G{zi) by 



Gi {zi, r'l, r'2) = G (zi, r'l, r'2) T (rl) T* (rl) , 



(57) 



where the transmission function T is defined by Eq. (25). One can obtain a 
law for the propagation of the cross-spectral density in free space using Eq. 
(53), in agreement with [8], from position zi immediately behind the lens to 
the image plane at position Z2, that is 
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Substituting Eq. (57) in Eq. (58) and remembering Eq. (25) one finds 
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(59) 



Manipulation of the argument in the exponential function under integral allows 
the more suggestive representation 
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(60) 



that is analogous to Eq. (26). The quantity in brackets (...) is basically Fourier- 
transformed. As mentioned before, similarities between the way G and E 
evolve through the beamline have to be ascribed to the fact that the average 
over random variables commutes with all other operations in the calculation 
of the cross-spectral density. As a result, the reason why Eq. (60) is basically 
a Fourier transformation is due to the particular way E evolves. 

Similarly as has been explained above for the fields, the image and the fo- 
cal plane are privileged planes for which the cross-spectral density assumes 
particularly simple forms, that can be found in terms of Fourier Optics. 



2.4-1 Large non-limiting aperture 

We will now proceed in analogy with Section 2.2.1. On the focal plane, simi- 
larly to Eq. (28) we have 
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of cf J ' 
while using Eq. (55) we obtain, similarly to Eq. (29), 
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(62) 



In analogy with Eq. (32), for the image plane we can write 



zum /.^ ,2 1^ |2^ 
2^ (l^^^l - l^^'l . 
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X j dr'idr'2 G{zi,r'i,r'2)exp 



X exp 
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(63) 



Using the convolution theorem in analogy with Eq. (34) we obtain 
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Finally, taking advantage of Eq. (55) we have 
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X J du dv Q (0, -u, {?) exp [im {fu ■ u — f2i ■ v)] , 
that can be rewritten as the analogous of Eq. (36) : 
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(66) 



^.^.^ Effect of aperture size 

Effects of aperture size can be included in strict analogy with Section 2.2.2. 
Similarly to Section 2.2.2 the following mnemonic rule can be applied to in- 
clude the effects of the pupil in Eq. (62) or Eq. (66). First, divide Eq. (62) 
or Eq. (66) by the first phase factor, corresponding to the phase factor out- 
side the integral sign in Eq. (60). Second, convolve twice with V and V*, V 
having already been defined in Eq. (39). Third, put the phase factor back. 
We will denote with Gp the cross-spectral density including the effects due to 
the presence of the pupil. On the focal plane, in analogy with Eq. (40), one 
obtains 
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Gp{zf,fif,f2f) = exp 
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while on the image plane, in analogy with Eq. (41), one has 
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(68) 



In the limit for large apertures, V and T* can be substituted by 5-Dirac 
functions in both Eq. (67) and Eq. (68). In this case, from Eq. (67) we recover 
Gp{zf,fif,f2f) — G{zf,fif,f2f), that is given by Eq. (62). Prom Eq. (68) 
instead, we have Gp{zi,'rii,f2i) — G{zi, 7^11,7^21), that is given by Eq. (66). 

Similarly to the case analyzed in Section 2.2, unless particular conditions are 
met, the phase factor under the integrals in Eq. (67) and Eq. (68) compen- 
sate only partially the phase in Eq. (62) and Eq. (66). In the following we 
will restrict our attention to the image plane. In this case, complete phase 
compensation is achieved when the lens is in the far field. 

The far field limit of Eq. (66) can be obtained by substitution of Eq. (56) in 
Eq. (63). After calculating the inverse transformation of Eq. (54) one obtains: 



G{Zi,fii,f2i)=TXl^exp 



imu ( 2 1^ 
,2^ (l^-l - 
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(69) 



Eq. (69) can also be obtained directly from Eq. (66) neglecting the second 
phase factor on the right hand side, that is possible when 



m^a; 

1CZ\ 



- \r2i 



^)«1 



(70) 



for any pair of points (rii,r2i) on the image pattern. When the coherence 
length on the image plane is much smaller than the characteristic size of the 
image, condition (70) constitutes, similarly to condition (43) before (that holds 
in the case of coherent light), the requirement to be satisfied for the lens to 
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be in the far zone. When condition (70) is satisfied, Eq. (69) holds instead of 
Eq. (66). As explained in Section 2.2, in this paper we will study pupil effects 
only under this assumption. 

Explicit substitution of Eq. (69) in Eq. (68) yields, the following far field limit 
expression: 



Gp{zi,fii,f2i) = e:iq) 



imUJ / 2 1^ |2\ 

.2^ (l^-l - l^-l ). 
^^"^''^^ v]g(o. 



du dv V 
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— u 



\ CZi J \ UJ UJ 

that is the analogous of Eq. (45) in Section 2.2.2. 



(71) 



3 Image formation with a perfectly coherent undulator source 



In order to give an explicit expression for the cross-spectral density at posi- 
tion Z2 we need to know an explicit expression for the cross-spectral density 
G{zi,r'i, r'2) immediately in front of the lens. This was the subject of our pre- 
vious work [2] . In this and in the following Sections we will present the image 
formation problem and its solution with the help of that reference. In particu- 
lar, in the present Section 3 we will begin with the simpler deterministic case 
of zero electron beam emittance. One may think of a filament electron beam 
or equivalently, as concerns the cross-spectral density, of a single electron. We 
will start neglecting the presence of the pupil. At the end we will generalize 
our results to account for it. In the following parts of this paper we will then 
generalize the results obtained in the present Section 3 to include emittance 
effects, thus taking full advantage of the Statistical Optics formulation. 

Our starting point is an expression, derived in reference [2], for the complex 
envelope Es± of the Fourier transform of the electric field produced by a single 
electron moving through a planar undulator at any distance from the exit of 
the undulator. That expression accounts for a given offset and deffection angle 
of the particle trajectory with respect to the undulator axis. Referring to Fig. 
2 we found 
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J Zn — Z' 
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X exp < t 



2c 
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(72) 
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Fig. 2. Illustration of the undulator geometry and of the observation plane after the 
undulator. 

Here K is the undulator parameter, L^, as before, is the undulator length, 
(— e) is the electron charge and 7 is the relativistic Lorentz factor. Moreover 



Jn being the Bessel function of the first kind of order n. Also, 
"° = A. (1 + if 72) 



is the fundamental frequency of the undulator, A^^ being the undulator period. 
Finally, 

C = (75) 



is the detuning parameter, which accounts for small deviations in frequency 
from resonance. 

Eq. (72) is valid for frequencies about the fundamental harmonic ujo- This 
means that we are considering a large number of undulator periods A^^ 1 
and that we are looking at frequencies near the fundamental at angles within 
the main lobe of the directivity diagram of the radiation. In this situation 
one can neglect the vertical ^/-polarization component of the field with an 
accuracy {AT^Nyj)~^. This constitutes a great simphfication of the problem. At 
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any position of the observer, we may consider the temporal Fourier transform 
of the electric field as a complex scalar quantity corresponding to the surviving 
x-polarization component of the original vector quantity. Moreover it should 
be noted that, in deriving Eq. (72), we assumed that no influence of focusing 

— * 

is present inside the undulator. ff and I are to be understood as deflection 
angles and offset of the electron at the position z — 0. 

Let us introduce normalized units ^ 



4^ ^ ^ 

KujeAjj ' 



lujL. 



w 



UJ — UJo 



C = L^C = 2nN,^ 
I UJ 

r±o = r±o ' 



LyjC 



UJ 



, LyjC 

i = — . (76) 

As shown in Appendix B of [2], after some algebraic manipulation, Eq. (72) 
can be rewritten in normalized units as 



1/2 

r dz' 
J Zo — z' 

-1/2 



exp < 



+ Cz' + 




2(^0 - z') 



,(77) 



where 



Q = ^ (78) 

Zo 

^ The relation between Egi, and Egi, in Eq. (76) differs from the analogous one in 
Eq. (25) of reference [2] for a factor = Zo/Ly,. The reason for this discrepancy 
is related to the different subjects treated. In [2] we considered only the free space 
case, while in this paper we extend our considerations to an optical element, thus 
introducing another privileged longitudinal position (the lens position) other that 
the observation plane. The definition of Eg±_ in [2] is no more a convenient one 
here and would lead to artificial complications in the following parts of this paper. 
Therefore it has been slightly modified as in Eq. (76). This leads to slight changes 
(related with the factor Zo) in some of the following equations when compared to 
the analogous quantities in [2]. 
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represents the observation angle and is given by 



^rr= \ 9-- 



2 



(79) 



Eq. (77) is of the form 



Es^\C,Zo,e - J -f] \ ^ exp {i^u)S 



C,Zo, I 6* - |- - 17 



(80) 



It is possible to show that the expression for the function S{-) reduces to a 
sinc(-) function as 3> 1 . In this limiting case, the expression for the electric 
field from a single particle, Eq. (77), is simplified to 
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Eq. (81) can be integrated analytically giving 

1 {C 
Es± = exp (i^c/)— sine 77 + -r , 



(81) 



(82) 



where 



Zn 



(83) 



and where the sine function has been defined as 

sin(a;) 



sinc(x) = 



X 



(84) 



For simplicity, in this paper we will restrict our attention to the case (7 = 0. 
In the particular case (7 = 0, the function S can be represented in terms of 
the exponential integral function Ei(-) as 



5(0,i„,C') =exp(-zi„CV2) 



Ei 



-1 + 24. 



2a2 



-Ei 



.1 + 24. 



(85) 
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Fig. 3. Comparison between f(Q = sinc(C^/4) (solid line), the real (dashed line) 
and the imaginary (dash-dotted line) parts of f{Q = ZoS{0, Zo, C^) at Zp = 1. 



It is interesting to study the behavior of the S function as the distance from the 
undulator center increases. This gives an idea of how good the asymptotic 
approximation of the S function for 3> 1 (that is a sine function) is. A com- 
parison between sinc(C^/4) and the real and imaginary parts of ZoS{0, Zg, C^) 
ioT Zo — 1, Zo — 2, Zo — 5 and io = 10 is given respectively in Fig. 3, Fig. 4, 
Fig. 5 and Fig. 6. 

When the electron beam has zero emittance we are dealing with a perfectly 
coherent wavefront. The evolution of the radiation wavefront through our op- 
tical system can be obtained with the help of Eq. (80). In the following we will 

— * 

study such evolution assuming C = 0, / = and fj = 0. These assumptions 
mean that the radiation frequency is perfectly tuned to the fundamental fre- 
quency of the undulator and that the electron beam is moving on the z axis. 
In this case Eq. (82) describes, in the far field region, a spherical wave with 
the source in the center of the undulator. This remark allows one to consider 
the undulator center as a privileged point. In other words, the phase factor in 
Eq. (82) represents, in paraxial approximation, the phase difference (charac- 
terizing a spherical wave) between the point {xg, jjo, Zo) and the point (0, 0, Zo) 
on the observation plane. 
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Fig. 4. Comparison between /((") = sinc((^^/4) (solid line), the real (dashed line) 
and the imaginary (dash-dotted line) parts of /(C) = ZoS{0, Zo, C^) at Zo = 2. 




Fig. 5. Comparison between /(C) = sinc(C^/4) (solid line), the real (dashed line) 
and the imaginary (dash-dotted line) parts of /(C) = ZoS{0, Zo, C^) at Zp = 5. 
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Fig. 6. Comparison between f(Q = sinc(C^/4) (solid line), the real (dashed line) 
and the imaginary (dash-dotted line) parts of f{Q = ZoS{0, Zq, C^) at Zp = 10. 

3.1 Large non-limiting aperture 



We will first study the case when the pupil function can be neglected. Let us 
introduce a normalized version of the spatial Fourier transform of the field, 
analogous to Eq. (14), that is 



(^z,iij — J dr'E ^z, r'^ exp 



(86) 



The spatial Fourier transform F i^z, uj can be calculated directly from Eq. (77) 
(compare also with Eq. (184) of reference [2]) and gives 
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(87) 



Eq. (28) and Eq. (34) can be respectively rewritten in normalized units as 
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(89) 



where / = f /L^. Substitution of Eq. (87) in Eq. (88) and in Eq. (89) yields 
results respectively for the focal and the image plane. In the focal plane we 
have 
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(90) 



Note that the relative distribution of intensity in the focus reproduces the 
angular distribution of intensity in the far field, that is 



(if) ■ 

A plot of the universal function sinc^(a2/4) is given in Fig. 7. Note that all ex- 
pressions pertaining undulator radiation from a single electron are azimuthal 
symmetric but do not admit factorization in the product of functions sepa- 
rately depending on the x and the y coordinates. As we have shown in [2] , the 
absence of factorization leads to an influence of the presence of the horizontal 
omittance on the coherent properties of undulator radiation in the vertical 
direction. As a result, even in the case of zero vertical emittance one cannot 
have perfect coherence in the vertical direction. 

For the image plane we obtain: 
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As already seen in Eq. (190) of [2] (Appendix C), the Fourier transform in Eq. 
(92) can be calculated in terms of a Fourier- Bessel transform: 
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Fig. 7. Universal function sinc^(a^/4) used to calculate (according to Eq. (91)) the 
focal intensity of a single electron at the fundamental harmonic at perfect resonance. 
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Fig. 8. Universal function ^'(q;) used to calculate (according to Eq. (96)) the in- 
tensity distribution on the image plane from a single electron at the fundamental 
harmonic at perfect resonance. 
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Eq. (92) can now be written as 



E±{zi,fi) = -imexp 



im 
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It is convenient to introduce the following universal function normalized to 
unity: 



^(a) = i- [tt - 2Si (a^) 



(95) 



The relative intensity on the image plane is related to the universal function 
^ through the scaling factor m: 



/(|fl|)=^'(m|fi|) . 
A plot of the universal function ^(a) is given in Fig. 8. 



(96) 



Comparison of Eq. (94) with Eq. (36) allows one to conclude that the equiva- 
lent source for a single electron moving on the z axis is a source characterized 
by a plane wavefront and an intensity distribution related to the universal 
function ^, that can be written as 



E^{Q,r) = -i [7r-2Si(| 



(97) 



By means of an inverse Fourier transformation it follows that 



F (0, ii) ^ -2m si 



smc 




(98) 



in agreement with Eq. (87). Comparison of Eq. (97) with a normalized version 
of Eq. (20) show that the phase of the field of the virtual source is shifted of a 
quantity — 7r/2 with respect to the spherical wave in the far zone. Such phase 
shift is the analogous of the Guoy phase shift in laser physics. A single electron 
produces a laser-like radiation beam that has a (virtual) waste much larger 
than the radiation wavelength located in the center of the undulator. 



The intensity distribution from the most elementary undulator source, i.e. the 
radiation from a single electron (or, equivalently, from an electron beam with 
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zero cmittance) was just described analytically. Such analytical description, 
Eq. (96), can immediately be applied in situations of practical relevance. In 
[10, 11] a characterization of the vertical emittance in Spring-8 is reported. It 
is based on the measurement of the X-Ray beam coherence length in the far 
zone. The experiment was performed at the beamline BL29XU. Based on the 
assumption of validity of the van Cittert-Zernikc theorem, it was found that 
the rms electron beam size at the undulator center (corresponding to the min- 
imal value of the beta function) was Sy ~ 4.5 /xm, and that the coupling factor 
between horizontal and vertical emittance was down to the value x ~ 0.12%, 
which corresponds to an extremely small vertical emittance Cy — 3.6 pm-rad. 
A resolution limit of this method was also discussed, based on numerical cal- 
culations of the radiation size from a single electron Sp ~ 1.6 ^m. at Ep = 14.41 
keV for the 4.5 m long undulator used in the experiment. The resolution limit 
of the measurement of Sy was estimated to be about 1 /im. 



Based on Eq. (96), we can determine the virtual source size of undulator 
radiation from a single electron. Let us consider the case when the single 
electron is emitting photons at the fundamental harmonic with energy E = 
14.41 keV. The angular frequency of light oscillations is given, in this case, 
by a; = 2.2 • 10^^ Hz. For an undulator length = 4.5 m, the normalization 
factor for the transverse size introduced in Eq. (76), (L^c/o;)-'^/^, is about 8 //m. 
From Fig. 8 obtain the dimensionless Half Width Half Maximum (HWHM) 
radiation size from a single electron (i.e. the HWHM width of the intensity 
distribution at the virtual source, located at the center of the undulator). This 
HWHM dimensionless value is about 0.7. It follows that the HWHM value of 
the radiation spot size from a single electron is about 0.7 • {cLyj/uY^'^ ~ 
6 /im. Therefore, the rms value Sp ~ 1.6 /im in [10, 11], which was calculated 
numerically [12], is an underestimation of the correct value. 



Note that the HWHM radiation spot size from a single electron is larger 
than the rms electron beam size Sy ~ 4.5 /xm found by means of coherence 
measurements. One concludes that the uncertainty due to finite resolution is 
larger than the measured electron beam size. This suggests that the method 
used in [10, 11] may be inconsistent. Such inconsistency may be traced to the 
fact that authors of [10, 11] assume the validity of the van Cittcrt-Zcrnike 
theorem in the vertical direction. If one assumes their result of a vertical 
emittance Sy ~ 0.3A/(27r), it follows a posteriori that the van Cittert-Zernike 
theorem could not have been applied in first instance (in this experiment the 
value of the beta function was f3 ~ L^^). Hence the inconsistency of the method 
follows. Analysis of experimental results should have been based, instead, on 
the study of transverse coherence for non-homogeneous undulator sources in 
free space made in [2]. 
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3.2 Effect of aperture size 



The effects due to the presence of the pupil can be included in the treatment 
by means of a normalized version of Eq. (40) and Eq. (41) on the focal plane 
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and on the image plane 
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where 



V{u) = j dr'P{?) exp 
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(101) 



Eq. (100) is valid independently of the position of the lens. However, as ex- 
plained in Section 2.2 and in Section 2.4 we will limit ourselves to the case 
when the lens is in the far zone. From Eq. (96) we see that the characteristic 
size of the source is of order unity, because \t' is a universal function. As a 
result, the far field zone for a single particle is defined by the condition 
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By substitution of Eq. (20) in Eq. (26), followed by use of the lens equation 
Eq. (30) and normalization, one obtains an expression for the field valid in the 
case the lens is placed in the far zone. Such expression is equivalent, in that 
limit, to Eq. (100). If the pupil function is not set to unity, we have 
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(103) 



Eq. (103) can also be obtained by substitution of Eq. (89) in Eq. (100) followed 
by application of the convolution theorem. 

In the far field case Eq. (100) (or Eq. (103)) can be more easily used to 
calculate the effects of the pupil on the intensity. A natural example to study 
is the case of a lens with azimuthal symmetry and no aberrations. Consider a 
pupil of radius a. After introduction of the normalized pupil radius: 



a — tti 
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we set the pupil function: 



P(f) 



1 if |f| < d 
otherwise . 



(105) 



Using the Fourier- Bessel transform one obtains 
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Substitution of Eq. (94) in Eq. (100) and use of the far zone assumption leads 
to 
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Eq. (107) is, essentially, the convolution product of the Fourier transform of 
two known functions with circular symmetry. Therefore, recalling Eq. (93) and 
Eq. (105), Eq. (107) can be written in terms of the following Fourier-Bessel 
transform: 
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Fig. 9. Comparison between the relative intensity for a single electron at the image 
plane Ip, Eq. (109), the universal function Eq. (96), and the (Airy) diffraction 
pattern from a circular hole, Eq. (110), as a function of mfj. Here a/zi = 1. 

Eq. (108) can also be directly obtained using the Fourier-Bessel integration 
formula and Eq. (103). Eq. (108) corresponds to a relative intensity 



/p(N) = j 



4zi 



J du uJq {Tii\fi\v}j sine f — 



(109) 



In the limit a/ zi ^ 1 Eq. (108) gives back Eq. (94) without the second phase 
factor, and Eq. (109) gives back Eq. (96) as it should be. 



When a/zi <^ 1 the sinc(-) drops out of the integral in Eq. (109) giving the 
diffraction pattern from a circular hole: 
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Further, analysis of Fig. 9 actually shows that Eq. (110) retains its vahdity 
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Fig. 10. Comparison between the relative intensity for a single electron at the image 
plane Ip, Eq. (109), the universal function Eq. (96), and the (Airy) diffraction 
pattern from a circular hole, Eq. (110), as a function of mfj. Here a/zi =5. 

also for a/zi < 1. 

Comparisons between the relative intensity Ip, Eq. (109), the universal func- 
tion ^, Eq. (96), and the diffraction pattern from a circular hole, Eq. (110), 
are plotted as a function of mfj in Fig. (9) for a/zi — 1, in Fig. (10) for 
a/zi — 5 and in Fig. (11) for a/zi — 10. 



4 Image formation with partially coherent undulator source 

4-1 Coherence properties of undulator source in the presence of electron 
beam emittance 

In the last Section wc dealt with the image formation problem in the case of a 
filament beam, i.e. when the electron beam emittance is zero. In this Section 
we will generalize the previous results to the case when the electron beam 
has finite emittance. In this situation, methods from Statistical Optics must 
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Fig. 11. Comparison between the relative intensity for a single electron at the image 
plane Ip, Eq. (109), the universal function Eq. (96), and the (Airy) diffraction 
pattern from a circular hole, Eq. (110), as a function of mfj. Here a/zi = 10. 

be applied in order to solve the image formation problem. As discussed be- 
fore, the cross-spectral density of an undulator source must first be calculated 
at the lens position and subsequently propagated through the lens and the 
forthcoming optical beamline up to the experimental plane. 

In [2] we proposed a method, based on Eq. (77), to calculate the cross-spectral 
density from undulator sources at any position in free space after the undula- 
tor. Let us follow [2] and use Eq. (77) to calculate the cross-spectral density. 
The cross-spectral density G is given, in dimensional units and as a function 
of dimensional variables, by Eq. (51). Since the field in Eq. (77) is given in 

normalized units and as a function of normalized variables Zo, 9x,y and C*, 
it is convenient to introduce a version of G defined by means of the field in 
normalized units 



G{zoA.02,c) = lEsAcrzo.e^ 
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Transformation of G in Eq. (51) to G (and viccversa) can be performed shifting 
from dimensional to normalized variables and multiplying G by an inessential 
factor 



G 



J- 



G . 



\KueAjj ^ 

Substituting Eq. (80) in Eq. (Ill) we obtain 



(112) 
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(113) 



77,; 



Expanding the exponent in the exponential factor in the right hand side of 
Eq. (113), one can see that terms in l\ y cancel out. Terms in Ql. y contribute 

for a common factor, and only linear terms in l^^y remain inside the ensemble 
average sign. Substitution of the ensemble average with integration over the 
beam distribution function leads to 



G'(io,^i,^2,C') = exp 



^^^^ I ^' ^ 



X exp 



x5* 



%{d2 -ei)-i 



s 







2- 




'i-U 




C, Zo, 













2n 



C, Zo, 



02- ^-V 



V 



(114) 



Here integrals dfj and in dl are to be intended as integrals over the entire plane 

— * j\ 

spanned by the r) and I vectors. Eq. (114) is very general and can be used as 
a starting point for computer simulations. 

We assume that the distribution in the horizontal and vertical planes are not 
correlated, so that F-^ = F. . If the transverse phase-space is specified 
at the virtual-source position Zg — corresponding to the minimal values of 
the beta functions, we can write F~ t — Ff,^Ft and F~ t = Ff,,Ft with 

'/x s^x Hy ■s''y '■y 
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27riV„ 



From Eq. (76) and Eq. (78) it is possible to see that 



2 ^Lyj 



N =a' — 



(115) 

(116) 
(117) 



where a^^y and a^'^y' arc the rms transverse bunch dimension and angular 
spread. Parameters N^^y will be indicated as the beam diffraction parameters 
and arc analogous to Fresnel numbers. They correspond to the normalized 
square of the electron beam sizes. Dx,y represent the normalized square of 
the electron beam divergences instead. It is also convenient to introduce the 
square of the apparent angular size of the electron beam at the observer point 
position Zo, that is 



x,y 



'-x,y 



(118) 



Substitution of relations (115) in Eq. (114) yields^ at perfect resonance [C — 
0): 



exp 



^2 ^o/2 



ATT^ D^DyN^Ny 



dfjx exp 



2D, 



In Eq. (119), for notational simplicity we substituted the proper notation 

G{zo, 9i, O2, C) with the simplified dependence G{zo, 9i, 62). This is justified because 

we will be treating the case (7 = only. Consistently, also S[zo, iO — 1/ Zq — 57)^] is 

to be understood as a shortcut notation for S[C,Zo, {6 — l/zo — "fff] calculated at 
(7 = 0. 
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exp 



-oo 
oo 



2D„ 



j di^ 



exp 



2N^ 



X / (i/j^ exp 



/2 



2iV„ 



exp 



xS 

Let us now introduce ^ 



01- --fj 
z„ 









2- 


s* 









, 9x 



9x1 + 9x2 



■ (119) 



(120) 



and 



9yl - 9y2 ^ _ 9yl + 9y2 



(121) 



With this variables redefinition we obtain 
G — G{zo, 9x, 9y, A9x, A9y) 



and, explicitly, 



exp 



i2zo (9xA9x + 9yA9y 



ATT-'^DxDyNxNy 



dr]x exp 



2D, 



X 



/ dfjy 

— OO 
OO 



exp 



X / aiy exp 

-oo 



2Dy 

27V„ 



oo 

/ 



dlx exp 



2A^. 



xS* 



■2i (A9xix + A9yiy) 
Zo, {dx + A9x-j-- Vx^ + (^9y + A9y -j-Vy^ 
Zo, (9x - A9x -J-Vx"^ + (^^y - ^^y - J ~ '^^) 



21 



(122) 



■ (123) 



^ Note that the definition of A9x and A9y differ for a factor 2 and a sign with 
respect to notations in [8, 9]. 



50 



A double change of variables fjx,y — ^ 'fjx,y + 0x,y followed by the substitution 
ix,y/ Zo — ^ ^x,y — 'f}x,y and by analytical calculation of the integrals in dfix,y 
leads to 



An'^^jAxbxAyby 



xexp 
xexp 

X 



i2exZoMx 



i26yZoA6y 



exp 



exp 



,exp 



91 + AAJ^AelDx + AiAx9xZoA9x 
2{Ax + bx) 

2{Ay + by) 

4>l + 24>x [h + 2iAxZoMx) 
2{Ax + bx) 



X exp 
xS 



<f>l + 2(l)y {9y + 2iAyZoA9y 
+ by) 

zo, {4>x - A9xf + (4 - A9yf 
Zo, {4>x + A9xf + {4>y + A9yf 



(124) 



Eq. (124) is a valid expression for the cross-spectral density in free space 
after the undulator device (i.e. for Zg > 1/2) and can be used together with 
equations from (60) to (66). 

Let us now introduce the dimensionless version of Eq. (54) with the help of 

^ =? -, 

r — Zo9 and Af — ZoA9 ^ : 



^ A short digression about Eq. (125) is due here. As the reader may have noticed, Q 
coincides with the Fourier transform, done with respect to r', of the Wigner distri- 
bution r,u) = J dAr' G{z, r, Af') exp[2z('u- Af')]. The knowledge of the Wigner 
distribution is mathematically equivalent to the knowledge of ^ or G. A formalism 
based on the Wigner distribution may be thus developed, which is mathematically 
equivalent to the one developed here. In the case of quasi-homogeneous sources, 
the Wigner distribution amounts to Eq. (6), that is the phase space distribution. 
Interpretations of such a function as a sort of generalized phase space distribution in 
more generic cases for non-homogeneous sources have been proposed. However, there 
is no practical advantage in considering such an approach in our case. Moreover, 
the Wigner distribution is a quantity that cannot be directly measured. Therefore, 
we prefer to use a formalism based on the cross-spectral density which is a physi- 
cally measurable quantity. The cross-spectral density may be directly measured by 
means of Young's double pinhole interferometer, whereas the Wigner function is a 
mathematical transformation of the cross-spectral density. 
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g (5, u, Auj ^ j df' dAr' G [z, f', Ar') exp \2i (u ■ Ar' + Au ■ r')] .(125) 
Its inverse is given by 



G (^z, f, Af 



(27r) 



du' dAu' 



xQ {z, u', Au'^ exp —2i (u' ■ Af + Au' ■ f) 



(126) 



Similarly as before, we consider coordinates 6 — r/zo and A^ = Af/zo in the 

limit Zo — ^ 00 but for finite ratios 9 and A^. The dimensionless version of 
Eq. (56) then reads: 



Gizo,e,Ae 



2izr,9 ■ A9 



g{0,-9, -A9 



(127) 



This result will be widely used in what follows. Moreover an analogous of Eq. 
(55) is: 



g (^z, 9, A9j = g (0, 9, A9j exp 



-i2z9 ■ A9 



(128) 



This result means that, aside for a phase factor, the spatial Fourier transform 
of the cross-correlation function, g, does not depend on z, as it follows from 
the analogous property of the Fourier transform of the electric field discussed 
in Section 2. 

Before proceeding, let us introduce the spectral degree of coherence g, which 
can be presented as a function of f and Af: 



(r, Af ) 



G (f, Af ) 



1/2 



4±(r + Af)| \ (\Es± (r- a 



1/2 



(129) 



With reference to Fig. 12, the modulus of the spectral degree of coherence, \g\, 
mathematically describes the fringe visibility of the interference pattern from 
a Young's double-pinhole interferometric measure. The phase of the spectral 
degree of coherence is related, instead, to the position of the fringes. The 
cross-spectral density gives amplitude and position of the fringes. In general, 
the process may not be quasi- homogeneous. In this case, the result of Young's 
experiment varies with f . In this case, the relation between the visibility V of 
the fringes and g{r, Af) reads 
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Fig. 12. Measurement of the cross-spectral density of an undulator source, (a) 
Young's double-pinhole interferometer demonstrating the coherence properties of 
undulator radiation. Radiation must be spectrally filtered by a monochromator or 
detector (not shown in figure), (b) In the quasi-homogeneous case the fringe visi- 
bility V of the resultant interference pattern is equal to the absolute value of the 
spectral degree of coherence: V = j^l. 



V^2 



Es± {r + Ai 



1/2 



Egj \r — Af 



1/2 



Es± (r + Ar) 
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In the quasi-homogeneous limit V 



Es± (r- Af) 
(Af) 
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(r, Af ) 
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4.2 Large non-limiting aperture 



As explained before we start neglecting, at first, the effects from a finite pupil 
dimension, assuming a perfect lens with no aberrations. The imaging problem 
for an ideal lens is solved once we find the cross-spectral density of the equiv- 
alent virtual source for the undulator source. On the focal plane we can write 
Eq. (62) in normalized units as 



Tf, Aff) = -—J- exp 



2i - 
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X exp 



j2 •' •' 



6^ 0, 



Tf Af J 



/ / 



while on the image plane, Eq. (66) in normalized units reads 



(131) 



G{zi,ri, Afi) =mexp 
X exp 



2im-, 

— Ti ■ AVi 

— Ti • Ar 

Z\ 



G (O, -mfj, -mAfj) 



(132) 



In all cases considered in this paper the position i = is well within the radi- 
ation formation length of the undulator. Therefore the cross-spectral density 
G, calculated at i = has no direct physical meaning, and must be considered 
as a quantity characterizing the virtual source only. From the definitions of 
virtual source and cross-spectral density follows that the virtual source pro- 
duces not only the same field but also the same cross-spectral density of the 
real undulator source, at any distance from the exit of the undulator. 



In the present study case of a radiation spot size smaller than the area of the 
lens and of a lens with no aberrations, Eq. (124) and Eq. (128), together with 
Eq. (131) and Eq. (132), solve the problem of characterizing the cross-spectral 
density on the focal plane (with the help of Eq. (131)) and on the image plane 
(with the help of Eq. (132)). The situation of a radiation spot size smaller 
than the area of the lens is practically achievable for Synchrotron Radiation 
due to its high directionality. In this case, vignetting effects are not present. 
However, even in this case, in order to use Eq. (131) and Eq. (132) one must 
further assume that aberrations can be neglected. 



J^.S Effect oj aperture size 



Accounting for the presence of the pupil, in analogy with Eq. (60) one has the 
following normalized expression for the cross-spectral density on any observa- 
tion plane at position zi along the beamline behind the lens: 



(7(^2, r,Af) 



1 

47r2(z2 - h) 



exp 



Hr ■ Ar 

Z2 - Zi 



X j dr' dA? I^g(^Zi, r', Ar' ]P{r' + Ar' ] P* [r' - Ar' 
1 



X exp 



f ^2-^1. 



r' ■ Ar' 
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Let us consider, more specifically, the focal and the image plane. Results can 
be obtained directly from Eq. (133). Alternatively, in analogy with Eq. (67) 
and Eq. (68), one can use our previous results, Eq. (131) and Eq. (132), divide 
them by the first phase factor, convolve them twice with V and V* and, finally, 
put the phase factor back. A normalized version of the cross-spectral density 
Gp including pupil effects at the focal is then found and can be written as 



Gp{zf,rf,Aff) = exp 

X du dAu exp —2ifu ■ Au G{zf, fu, fAu) 



xV 



— u — Ait 



— u + Au 

(134) 



while on the image plane one obtains 



Gp{zi,ri,Afi) = exp 



2im 



fi ■ Af,: 



du dAu exp 



xdlziAuAA^ir 

\ m m / 



u ■ Au 



m 



^ (ri + Af -u - Au 



IZl 



xV* 



— (ri — Af j) — u + Au 

IZi^ ' 



(135) 



As said before, we will treat particular situations in the image plane when the 
lens is in the far field. Using coordinates and Af j the analogous of condition 
(70) reads 



-^ri • Ari < 1 

Z\ 



(136) 



for any pair of points on the image pattern. 

Explicit substitution of Eq. (132) in Eq. (135) yields the following far field 
limit expression, which accounts for condition (136): 



Gp(ij,fj, Afj) = mexp 
xV 



2im 



r,- ■ Af ,• 



Z\ 

— [fi + Af j \ —u — Au 

z\^ ' 



J du dAu G (o, —Ziu, —ziAuj 
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— (ri - Afi) - u + Au 



(137) 



Eq. (137) is in analogy with Eq. (71) and Eq. (45). 



5 Imaging of quasi-homogeneous Gaussian undulator sources by 
a lens with Icirge non-Umiting aperture 



In this Section wc spcciaHzc our discussion to the particular case of quasi- 
homogeneous Gaussian undulator sources, assuming a lens with large non- 
limiting aperture and no aberrations. A Statistical Optics treatment is not 
the only one possible in this particular study case. A Geometrical Optics 
approach can also be applied, practically consisting in ray-tracing techniques. 
In this Section we will consider the image formation problem from a Statistical 
Optics viewpoint. In Section 7 an analysis in terms of Geometrical Optics will 
be given, and agreement between these two methods will be demonstrated. 

From this point on, we will systematically ignore unimportant pre-factors ap- 
pearing in the expressions for the cross-spectral density. Moreover we will 
assume iVj, ^ 1 and ^ 1 which is a reasonable approximation for third 
generation light sources in the X-ray region. We will show that this assump- 
tion leads to a major simplification: namely, horizontal and vertical coordi- 
nates turn out to be factorized in the expression for the cross-spectral density 
in free space. As a consequence, Eq. (131) and (132) can also be factorized in 
the product of a factor depending on the horizontal coordinates and a factor 
depending on the vertical coordinates. These separate factors will be obtained 
from Eq. (131) and (132) substituting all vector quantities with scalar quan- 
tities (horizontal or vertical components). 



5. 1 Evolution of the cross-spectral density function in free space 



Eq. (124) is a vahd expression for the cross-spectral density in free space at 
perfect resonance, calculated under the only assumptions that the system is 
ultra-rclativistic (and, therefore, the paraxial approximation can be applied) 
and that the insertion device is characterized by a large number of undulator 
periods. In this case the resonance approximation is enforced. Eq. (124) is 
quite generic and, with respect to first-principle calculations, it involves the 
computation of a two-dimensional integral, whereas the most generic calcula- 
tions would require a total of six integrations, two over the undulator length 
and four over the electron beam transverse phase space (assuming that the 
cross-correlation terms between different electrons is neglected). From a com- 
putational viewpoint, the advantage of reducing the number of integration is 
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obvious and it can be appreciated even more after the cross-spectral density is 
propagated through an optical system with limiting apertures, which naturally 
increases the dimensions of the integration to be performed. 

When TVa; 3> 1 and Dx ^ I the cross-spectral in free space, Eq. (124), can be 
written as the product of factors separately depending on the x and on the y 
coordinate, as has been shown in [2]. In fact, analyzing the exponential factor 
outside the integral sign in Eq. (124) it is possible to see that the maximum 
value of A9l is of order (A^ + Dx)/{AxDxzl) ^ 1, where we remember — 
Nl/z^. As a result, AO^ can be neglected inside the S functions in Eq. (124). 
Moreover, since Dx ^ 1 one can also neglect the exponential factor in 0^ -|- 
2(px0x inside the integral. This leads to 



G = exp i29xZoA9x 
X exp i29yZoA9y 

oo 

X / a0x exp 



exp 
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j d^x 

-oo 

oo 
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+ 2(j)y (9y + 2iAyZoA9y) 

2(4 + 4) 



xS* 



ZO, + (4 - ^Oyf S Z„ + {4)y + A9y) 



(138) 



Following the same reasoning in [2] , we can also neglect the phase factor in ^x 
under the integral in d(^x in Eq- (138). As a result, when A^^- ^ 1 and Dx ^ 1 
horizontal and vertical coordinates are factorized and we obtain the following 
equation for G: 



G{zo, 9x, 9y, A9x, A9y) — Gx{zo, 9x, A9x) Gy{zo, 9y, A9y) , 



(139) 



where 
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(140) 



and 



Gy = exp 



i29yZoA9y 



exp 
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2{Ay + Dy) 

Zo, 4>l + {4>y - Aeyf]s\zo, 4>l + {4>y + ^Oyf 



(141) 



To begin our investigation of quasi-homogeneous sources we will consider the 
limit N ^ 1 and ^ 1, when the photon-beam phase space in a certain 
(horizontal or both horizontal and vertical) direction is an exact replica of the 
electron-beam phase space. Calculations can be performed in one dimension, 
suppressing indexes x and y. In the case of Second Generation light sources 
the results that wc arc going to derive constitute a realistic description of the 
radiation characteristics in both horizontal and vertical directions. Then, Eq. 
(141) coincides with Eq. (140). 

Let us present Eq. (140), i.e. the asymptotic expression for the cross-spectral 
density in the limit iV ^ 1 and Z) ^ 1, in terms of coordinates f — z9 and 
Af = zA6. We have 



G{z, r, Af) =exp 



where 
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X exp 



2{A + D)z\ 
AfAf 





r Af 


exp 
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-2i- 



z{A + b) 



exp 



-2 



AD{Aff 

{A + b) 



(142) 



(143) 



Note that here, depending on the situation, r may assume the meaning of 
either variable x or y. 

In the far field limit, when A ^ D one obtains the following limiting expression 
of Eq. (142): 



G{z,r, Af) =exp 



2i- 



rAf 



exp 



2Dz^ 



exp 



-2A(Af)2 



(144) 



With the help of Eq. (144) and using Eq. (127) one can find the expression 
for Q and G at the virtual source position in the center of the undulator. In 
the case under study (A^ ^ 1 and D ^ 1) the virtual source is a Gaussian 
quasi-homogeneous source. Aside for unessential multiplication constants we 
have 
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g (o, e, Ae) = exp 



exp 



2b 



(145) 



Therefore, using Eq. (126) we also obtain 



G(0, f, Ar) = exp 



^2 



2Ar 



exp 



2D(Af)^ . (146) 



From Eq. (146) we conclude that the intensity distribution of the virtual source 
is a replica of the electron beam density distribution at the position of minimal 
beta function of the undulator (i.e. at the undulator center). Moreover, in this 
particular study case, if the position of the minimal beta function does not 
coincide with the undulator center, the virtual source corresponding to the 
description in Eq. (146) is simply translated, and is always located at the 
position where the beta function of the electron beam is minimal. 

It should be noted that the far field limit A <^ D corresponds with the ap- 
plicability region of the van Cittert-Zernike theorem. In virtue of the van 
Cittert-Zernike theorem the modulus of the spectral degree of coherence in 
the far field, i.e. cxp[— 2/l(Af)^] from Eq. (144), forms a Fourier pair with the 
intensity distribution of the virtual source, i.e. exp[— f^/2A^] from Eq. (146). 

In particular one concludes that the rms width of the virtual source is 
as it can be seen directly from Eq. (146). In our study case for TV » 1 and 
D ^ 1, such a relation between the rms width of the spectral degree of co- 
herence in the far field and the rms dimension of the virtual source is also a 
relation between the rms width of the cross-spectral density function in the 
far field and the rms dimension of the electron beam at the plane of minimal 
beta function in the center of the undulator. In dimensional units one can 

— * 

write the value Uc of the rms width of the spectral degree of coherence g{Af) 
in the far field as 



a being, as usual, the rms dimension of the electron beam. These few last re- 
marks help to clarify what is the size of the source in the van Cittert-Zernike 
theorem, that is far from being a trivial question. For instance, assume that the 
van Cittert-Zernike theorem can be applied. Then, the rms electron beam size 
can be recovered from the measurement of the transverse coherence length. In 
this regard, in [10] Section V, one may find a statement according to which 
the rms electron beam size "is only the average value along the undulator" be- 
cause "the beta function has a large variation along the undulator". However, 
as we have seen before, the concept of virtual source does not require a small 
variation of the beta function. In the most general case, any variation of the 



59 



beta function does not affect the virtual source size and, in our case of quasi- 
homogeneous Gaussian source, the virtual source size is also the transverse size 
of the electron beam at the position where the beta function is minimal. An- 
other example dealing with the same issue is given in reference [13]. This paper 
(as well as reference [10]) reports experimental results. However, authors of 
[13] observe a disagreement between the electron beam rms size reconstructed 
from the van Cittert-Zernike theorem and beam diagnostics result of about a 
factor 2. They ascribe this variation to the variation of the electron beam size 
along the undulator. In footnote [25] of reference [13], one may read: "The 
precise shape and width of the x-ray intensity distribution in the source plane 
are directly connected to the properties of the electron beam. It would not be 
surprising if the limited depth of focus of the parabolically shaped electron 
beta function in the undulator translates into a virtually enlarged x-ray source 
size." . At first glance it looks like if the Synchrotron Radiation source has a 
finite longitudinal dimension. However, based on the previous discussion we 
conclude that the virtual source size is equal to the electron beam size at the 
point where the beta functions have their minimum and that it is not affected 
by variations of the beta function along the undulator. As a result, one should 
not observe any virtually enlarged X-ray source size because of this reason. 



5.2 Evolution of the cross-spectral density function behind the lens 



If we neglect the effect of the pupil in the particular case under examination 
{N ^ 1 and D 3> 1), it is possible to find an analytical expression for the 
cross-spectral density for any observation plane, and not only for the focal 
or the image plane. This is due to the fact that, in this particular case, the 
virtual source is gaussian. In the most general case instead, one has to make 
use of Eq. (133). 

As usual we will neglect, at first, the effect of the pupil function and group all 
the phase terms in AfV in Eq. (133) with the help of the definition 

Q = ^ - i + , (148) 



where 



A^y,. (149) 

Zl 



With this in mind, after substitution of Eq. (142) calculated aX. z = Zi in Eq. 
(133), one obtains 
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' [2AD + 2{A + DYQH^{z2-ziY 

and to a modulus of the spectral degree of transverse coherence 



(151) 



\g{z2,r,Ar)\=exp 



2Ad{A + D)zf{Af) 



AD + {A + DYQHf {z2-z,y 



;i52) 



Letting Q = —A/[zi{A + D)] the reader can specialize the results to the case 
of the image plane. For Q — 1/zi — A/[zi{A + D)] one gets the results for 
the focal plane. Also, the intensity and the modulus of the spectral degree of 
coherence on the image plane can be obtained from those on the focal plane 
exchanging A with D. This symmetry can be explained in terms of Fourier 
transforms. Phase factors aside, the cross-spectral density on the image plane 
is equal to the cross-spectral density on the object plane. The cross-spectral 
density on the focal plane instead, is equal (phase factors aside) to the Fourier 
transform of the cross-spectral density on the object plane. 

As we have seen, in the case for a Gaussian electron beam with iV 3> 1, 
D ^ 1 and for a perfect lens with non-limiting aperture and no aberrations, 
the Gaussian approximation for the cross-spectral density at the virtual source 
in Eq. (146) can be used, and the cross-spectral density in free space at any 
position z can be calculated with the help of Eq. (142). Then, Eq. (133) can 
be simplified to recover both the intensity and the modulus of the spectral 
degree of coherence (Eq. (151) and Eq. (152) respectively), which are Gaussian 
functions for any value of Z2 and Zi. Even for quasi- homogeneous sources 
though, there are a number of examples when it is difficult to obtain analytical 
results from Eq. (133) for any value of Z2. Nevertheless it is possible to calculate 
the cross-spectral density at the image plane and at the focal plane (for any 
value of Zi) with the help of Eq. (131) and Eq. (132). This can be done relying 
on the calculation of the cross-spectral density at the virtual-source position 
(and its Fourier transform), which allows further use of Eq. (131) and Eq. 
(132). We will first use Eq. (131) and Eq. (132) to deal with the case that we 
just discussed when N ^ 1 and D S> 1. This is not a simple repetition of 
already known results, because the particular way of reasoning used for the 
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focal and the image plane, through Eq. (131) and Eq. (132), will be widely 
used in the following parts of this paper too. The Statistical Optics method 
conjugated to Fourier Optics results allows us to predict, by manipulations 
of Eq. (142), the cross-spectral density (and, therefore, the intensity and the 
absolute value of the spectral degree of coherence) on the focal and on the 
image plane. In order to use Eq. (131) and Eq. (132) we must take advantage 
of the expressions for Q and G at the virtual source position, Eq. (145) and 
Eq. (146) respectively. 

With the help of Eq. (131) and Eq. (145), on the focal plane we obtain 



G{zf,ff,Aff)^exp 



^2(/-^i) - 
r fAr f 

j2 J J 



exp 



2NArj 

P 



exp 



2DP_ 

(153) 



For any value of Zi we have a relative intensity on the focal plane given by 



^2 ^ 



i — exp 



2b p\ ' 



(154) 



while the modulus of the spectral degree of coherence (again, for any position 
zi of the lens ^° ) is 



|5(%,r/, Af/)|=exp 

These results are intuitively sound. In Section 2 we explained that we expect 
to find, on the focal plane, the spatial Fourier transform of the wavefront on 
the object plane (except for a phase and a proportionality factor). Therefore, 
it is intuitive that the intensity on the focal plane must depend on the electron 
beam divergence only and that the modulus of the spectral degree of coherence 
must depend on the electron beam size only. In fact, the exchange of roles of 
and D passing from the virtual source plane to the the focal plane is related 
to the operation of Fourier transform. Also note that the Fourier transform of 
the field depends on zi through a phase factor only, and free space basically 
acts as a Fourier transform itself (see Section 2): what we find on the focal 
plane in terms of intensity and modulus of the spectral degree of coherence we 

Note that the modulus of the spectral degree of coherence in Eq. (155) is inde- 
pendent of fj. However, the spectral degree of coherence depends on fj through a 
phase factor. This situation corresponds, according to a definition given by us in [2], 
to a weakly quasi-homogeneous wavefront. This remark is vahd for many expressions 
of the modulus of the cross-spectral density given in this work. 



2NAfj 



(155) 
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must also find in the far field after propagation in free space. The reader may 
check that, after substitution ff/f — > 6 and Arj/f — > A^, Eq. (154) and 
Eq. (155) can be found from Eq. (54) and Eq. (61) in [2] in the limit Zo ^ 1, 
which describe propagation in free space, as it should be. 

A similar simplified reasoning can be applied for the cross-spectral density on 
the image plane. With the help of Eq. (132) and Eq. (146), describing the 
cross-spectral density of the virtual source we obtain: 



G{zi,ri,Ari)^ex-p 



2im(m -|- l)riAfj 



X exp 



2 





2 — 2 " 


exp 




2N _ 



2Dm\Afiy 



The relative intensity on the image plane is given by 
/ = exp 



2N 



(156) 



(157) 



while the modulus of the spectral degree of coherence is 



\g {zi,ri,Afi)\=exp -2m^Z)(Af 



(158) 



These results are very natural. By definition of image plane, when we image an 
object with an ideal lens with a large non-limiting pupil aperture, we obtain 
a magnified version of the object (virtual, in this case). 

We remarked before that the intensity and the modulus of the spectral degree 
of coherence on the image plane can be obtained from that on the focal plane 
exchanging A with D. This symmetry though is not evident from the expres- 
sions in Eq. (154), Eq. (157), Eq. (155) and Eq. (158): to display it one has 
to express these equations in terms of Z2 and zi. 



6 Imaging of queisi- homogeneous non-Gaussian undulator sources 
by a lens with leirge non-Umiting aperture 



In the previous Section 5 we treated the case for TV 1 and » 1. In 

the present Section 6 we will deal with other quasi-homogeneous cases, always 
assuming A^^ ^ 1 and ^ 1. The quasi- homogeneous situations that remain 
to be treated under this assumption are for either Ny ^ 1 or Dy ^ 1. In fact, 
the situation for both Ny^l and Dy^ lis automatically included in Section 
5. Moreover, in all quasi- homogeneous cases, the cross-spectral density in the 
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horizontal direction obeys Eq. (150). Therefore, we will focus our attention on 
the cross-spectral density in the vertical direction only. 



6.1 Source with non- Gaussian angular distribution in the vertical direction 



Let us use of our Statistical Optics approach to solve a somewhat complicated 
image formation problem. After assuming separability of the horizontal and 
vertical directions [N^ ^ 1, -Da- 1) we suppose that the electron beam has 
a vertical transverse size much larger than the diffraction size, TVj^ ^ 1, and 
a finite divergence Dy > 0. As usual, we will first neglect the influence of the 
pupil function. The difference with respect to the case treated in the previous 
Section 5 is that Eq. (133) cannot be explicitly calculated for any value of 
Zi and Z2. However, as said before, the Statistical Optics method conjugated 
to Fourier Optics results allows us to predict, for any value of ^i, the cross- 
spectral density on the focal {z2 = Zf) and on the image {z2 — Zi) plane by 
means of Eq. (131) and Eq. (132). In order to use these equations we must 
first calculate G and Q at the virtual source at position Zo = 0. This can be 
done taking the limit z'^ ^ N/D of Eq. (141), i.e. calculating the far zone 
hmit of Eq. (141), and using Eq. (127). First, under the assumption Ny 1 
we can neglect A^^^ in the S functions in Eq. (141), thus obtaining the vertical 
cross-spectral density function in the far field limit 



G = exp i29yZoA.0y 



exp 



-2NyA9l 



exp 



{4>y + 9y) 



2D„ 



(159) 



where 



Wy) = ^ y ^'^^ [{^l + ^l) 1^. 



(160) 



is a universal function related to undulator radiation. A plot of Is is given in 
Fig. 13. Eq. (159), substituted into Eq. (127), gives the Fourier transform of 
the cross-spectral density at Zo — 0, i.e. at the virtual-source position: 



Q (0, M, A-u) — exp 



-2NyAu^ 



^ d(i>y 



exp 



2A, 



h{4>y) 



(161) 
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Fig. 13. The universal function Is, used to calculate the focal intensity of a 
quasi- homogeneous source at ^ I, Dx ^ ^ and iV^ S> 1. 

Inverse transforming Eq. (161) according to the definition in Eq. (126), we 
obtain the cross spectral density at the virtual source position 



G'(0,y,Ay) =exp 



2K, 



exp 



-2DyAf j{Ay) , 



(162) 



where function '~f{Ay) is an inverse Fourier transform, normalized to unity, of 
Is, defined as from Eq. (93) in [2]: 



1 °° 

7(A^) = ^ / [^(-2A^)0j/s(0j/) 



(163) 



It has been shown in [2] that 7 can be expressed in terms of the sine integral 
function Si(-) and of the cosine integral function Ci(-). One has 



7(Ay) 



TT 



TT 



+ 



2Ay2Ci (2Af) - sin (2Af) - Si (2Af)\ . (164) 
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This means that 7 is a real function. Moreover, in Eq. (162), and y are 
separated and, since A^j, ^ 1, the typical correlation length is much smaller 
than the radiation spot, independently of the value of Dy. This shows that 
Eq. (162) models a quasi-homogeneous source. From Eq. (131) and Eq. (161) 
we obtain the cross-spectral density on the focal plane 



G'(^/,y/, Ay/)=exp 



2i 



P 

00 

X / aOy exp 



(/ - ^1) y/Ay/ 



exp 



-00 



{vf/f + k 
2A, 



P 
Wy) ■ 



(165) 



The relative intensity on the focal plane is therefore given by 



, exp 



{vf/f + k 



2L»„ 



isik) 



00 



X < / o(^j,exp 

v.— 00 



2D„ 



isik 



(166) 



while the modulus of the spectral degree of coherence reads 



b(%,y/, Aj//)|=exp 



-2- 



P 



(167) 



For the image plane, Eq. (132) and Eq. (162) give the following cross-spectral 
density: 



G(ii,yj, Ayj) =exp 



im(m + l)yiAyi 



X exp 



2zi 
' 2K, 



exp 



-2Dym'Ayf ^{mAyi) , 



(168) 



corresponding to a relative intensity on the image plane 



I{zi,yi) = exp 



27V", 



y -I 



(169) 



Eq. (169) is the (magnified) image of the electron beam in the object plane 
z — 0. The modulus of the spectral degree of coherence is: 
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Fig. 14. Absolute value of the universal function 7, used to calculate, on the image 
plane, the spectral degree of coherence of a quasi-homogeneous undulator source 
when iV^ » 1, Z)^ » 1 and Ny » 1. 



g(zi, yi, Ayi) \ = exp -2Dym Ay^ \-f{mAyi) \ . 



(170) 



A plot of |7(Ay)| is given in Fig. 14. 



6.2 Source with non-Gaussian intensity distribution in the vertical direction 



Let us now consider the case when Dy ^ 1 and Ny assumes arbitrary values. 
In this case, Eq. (115) and Eq. (125) of [2] allow reconstruction of the cross- 
spectral density in the far zone, that is 



G — exp 



2ieyZoAey 



exp 



-2NyA9l 



exp 



2D. 



(171) 



where the function P{Ay) is defined in Eq. (113) of [2] and reads: 
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oo oo 
/5(A^,) = — J d4>y J dk 



xsmc 



smc 



(172) 



Eq. (171) may be obtained directly from Eq. (141) in the limit zl ^ Ny/Dy, 
i.e. in the far zone. Note that if Ny/ Dy <C 1 [Ny ^ 1 is our main case of interest 
since we have already treated the case when both iV^^ S> 1 and Dy ^ 1) the 
far zone begins already at the exit of the undulator, when ii ~ 1 (see also [2] 
for details). Let us introduce, in analogy with Eq. (163) the following (inverse) 
Fourier transform of the function 13: 



^{y) = -jQ j d(j)yexp [i{-2y)(j)y\ (3{(f)y) . 



(173) 



Here JC is the normalization factor 



oo 

/C = d4)yl3{4)y) ~ 2.200 , (174) 

— oo 

and has been calculated numerically. 

Both /3{A9y) and B{y) admit representations in terms of a one-dimensional 
integral (note that the representation for l3{A9y) has been already introduced 
in [2]). In order to see this, let us first consider the function: 



oc oo 



/(A^;,A^;) = — J dcpy J dcp. 



-oo — oo 



xsmc 



xsmc 



(4 + A9U2f + {4>y + A9y2f 



(175) 



The function / is circularly symmetric. This can be seen switching to polar 
coordinates: 



4>x^f^ cos(i7</,) 
0j/ = f<^sin(?7<^) 



(176) 



and 
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Then, Eq. (175) can be rewritten as 



(177) 



oo 2tt 
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dr^ J dfj^ 



rl + fo- 2r^f0 cos (?7^ - rje) 



rl + re + 2r^r0 cos {fj^ - fje) 



(178) 



which does not depend on rje, as can be seen switching to the integration 
variable rj' — rj^j, — rje- The following relation follows: 



(179) 



for any (A^^, A^^) such that 



(180) 



The function P can be seen as a restriction of the function /. The reason why 
/ has been introduced is that it allows the use the autocorrelation theorem to 
obtain the following relation: 



/ dAO', J dA^;exp[z(a.A^; + a,A^;)]/(A^;,A^;) 



-oo — oo 

oo oo 



27r2 



j d^x j d4)yexj£)[i{aj:^j: + ay^y)]smc 



(181) 



The integral in the right hand side of Eq. (181) has been already calculated, in 
practice, in Eq. (93), and it can be expressed in terms of the universal function 
^ defined in Eq. (95). It follows that 



otj oo 

j dAe', j dA^;exp[^(a,A^^; + a,A^^;)]/(A^;,A^;) 

-oo — oo 

TT - 2Si(a^ + aD] ' = 277^* {^^^^+^ ■ 



(182) 
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We can now inverse transform Eq. (182) using the Fourier-Bessel formula, thus 
obtaining 



(183) 



Letting A9'^ = and using Eq. (179) and Eq. (180) we obtain the following 
representation for jS: 



P{My) J da aJ^ (2aA^J *(«) . 



Applying the definition of B in Eq. (174) we obtain 



(184) 



CXJ CXJ 

B{y) da a^{a) j dAOyBxp [i{-2y)A9y\ Jo (2aA9y) . (185) 



The Fourier integral in A9y can be performed analytically (see [8], Appendix 
A. 3.), thus giving the following representation of B: 



da 



rect [y/{2a) 
1 - {y/a) 



^(«) , 



(186) 



where the function rect(a;) is defined, following [8], to be unity for \x\ ^ 1/2 
and zero otherwise. A plot of the universal B is given in Fig. 15. 

Using Eq. (127) and Eq. (171) we obtain the Fourier transform of the cross- 
spectral density at the virtual-source position 



Q (0, M, Au) = exp 



-2NyAu' 



exp 



u 



2D„ 



KM 



(187) 



Inverse transforming Eq. (187) we can write the cross-spectral density at the 
virtual-source position 



CXJ 

G{0,y,Ay) = exp [-2DyAf] j drjexp 



iv + yY 

2K. 



(188) 
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Fig. 15. The universal function B, used to calculate intensity on the image plane of 
a quasi- homogeneous undulator source when ^ I, Dx ^ ^ and Dy^^l. 

Then, from Eq. (131) and Eq. (187) we obtain the cross-spectral density on 
the focal plane 



G{zf,yf,^yf)^eW 



2i 



If 
X exp 



^ (/ - ^i) Vf^Vf 



exp 



(189) 



The relative intensity on the focal plane is therefore given by 



H^f^Vf)^ exp 



2DyP 



(190) 



while the modulus of the spectral degree of coherence reads 
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Fig. 16. Absolute value of the universal funetion f3, used to calculate coherence on 
the focal plane of a quasi- homogeneous undulator source when 3> 1, Dx S> 1 
and Dy » 1. 



If(%>y/>^y/)l=exp 



'if. 



(191) 



A plot of \P{A9y) \ is given in Fig. 16. For the image plane, Eq. (132) and Eq. 
(188) give the following cross-spectral density: 



im{m + VjyAm 



X exp 



2zi 
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corresponding to a relative intensity on the image plane 



{fj + myi 
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The modulus of the spectral degree of coherence is 



(193) 



g{zi, yi, Ayi) \ = exp -2bym^Ay1 



(194) 



Finally, in the particular case for Ny <^ 1, Eq. (193) reduces to 



I{zi,yi) = B{myi) . 



(195) 



7 Analysis of the image formation mechanism for quasi-homogeneous 
undulator sources in terms of Geometrical Optics 

In the Introduction we have stressed that the image formation problem is 

twofold: one should be able to provide a characterization of the virtual source 
as well as to track the cross-spectral density of the source trough the optical 
beamline. 

Let us first analyze the problem of source characterization. In Section 1 we 
have seen that, in the asymptotic limit for a large electron beam emittance 
^x,y ^ '^/(27r), Geometrical Optics may be used equally well as Statistical 
Optics to fulfill this task. Here we will discuss more in detail the relation 
between the Statistical Optics approach and the Geometrical Optics approach 
with particular attention to the applicability region of the latter. 

Let us start with a remark, which applies not only to undulator radiation 
sources but also to sources of other kind (e.g. bending magnets). In the In- 
troduction, in order to decide wether Geometrical Optics or Wave Optics 
is applicable, we compared the electron beam emittance with the radiation 
wavelength. This is acceptable in many cases when undulator radiation is in- 
volved but not, for instance, when bending magnet radiation is considered. 
In all generality one should separately compare the photon beam size and di- 
vergence with the radiation diffraction size and diffraction angle, which are 
quantities pertaining the single electron radiation. Let us fix a given direction 
X or y. The square of the diffraction angle is defined by (o"^)^ ~ A/ (27rL/), Lf 
being the formation length of the radiation at wavelength A as defined in [1]. 
The diffraction size of the source is given by a^, ~ f^'d^f- In calculating the 
photon beam size and divergence one should always include diffraction effects. 
As a result, if and (c')^ indicate the square of the electron beam size and 
divergence, the corresponding square of the photon beam size and divergence 
will be respectively of order max[cr^, crj] and ma.x[{a')^, (cr^)^]- These quantities 
can be rewritten in terms of the electron beam emittance as max[e/3, crj] and 
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max[e//3, (c^)^], (3 being the minimal beta function value, defining the virtual 
source position for the radiator (undulator, bending magnet, or other). Di- 
viding these two quantities respectively by cr^ and (cr^)^ give natural values, 
normalized to unity, for the photon beam size max[27re/?/(LjA), 1] and diver- 
gence max[27reL//(/3A), 1]. When the product between these two quantities is 
much larger than unity one can use a Geometrical Optics approach. In this 
case, this product represents the normalized photon beam emittance. When 
P Lf, as in many undulator cases, one may compare, for rough estimations, 
the electron beam emittance and the radiation wavelength as we have done 
before. However, in the case of a bending magnet one may typically have /3 of 
order 10 m and Lf of order 10~^ -r- 10~^ m. The ratio ^ 1 now consti- 

tutes an extra large parameter of the problem. In this case, even if the electron 
beam emittance is two order of magnitude smaller than the wavelength, due to 
diffraction effects one can still apply a Geometrical Optics approach, because 
max[27re/3/(LyA), 1] •max[27reLj/(/3A), 1] » 1, i.e. the photon beam emittance 
is much larger than the wavelength. As a result, dimensional analysis suggests 
that bending magnet radiation may be treated exhaustively in the framework 
of Geometrical Optics even for third generation light sources. 

As discussed above, when /3 ~ L/ a large electron beam emittance (com- 
pared with the radiation wavelength) is a necessary and sufficient condition 
for the Geometrical Optic approach to apply. In spite of that, when j3 ^ Lf 
or /9 <^ Lj, a large electron beam emittance is a sufficient, but not neces- 
sary condition for the Geometrical Optic approach to be possibly used for 
source characterization. Let us prove this statement with undulator sources 
in mind . It is enough to prove that the wider class of quasi-homogeneous 
sources, which includes situations when the electron beam emittance is not 
larger than the wavelength, can be described in terms of Geometrical Optics. 

Let us then consider the class of quasi-homogeneous virtual sources for undu- 
lator devices. The cross-spectral density of the virtual source (positioned at 
2; = 0, i.e. at the virtual source plane) can be written as in Eq. (2), that we 
rewrite here for convenience in terms of coordinates fx,y and ^fx,y'- 

Go{fx, fy, Af^, Ary) = i (fa;, fy) ^'(Af^, Afy ) . (196) 

As usual, the Fourier transform of Eq. (196) with respect to all variables will 
be indicated with 



Note that, even though in the case of undulator sourees one often has P Lf, 
there are situations when f3^LfOTf3<^Lf and when the electron emittance is 
of order of the wavelength. However, in the undulator case, very large values of the 
ratio /3/Lf of order 10^ 10^, typical of the bending magnet case, are unrealistic. 
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— oo — oo — oo — oo 

X exp[2i(^,Af; + dyAf'^)] exp[2i(Ae;f^ + A0>;)] . 

(197) 

The two quantities i{f^J-y) = Go{fx,fy, 0,0) and r{6^,6y) = ^o(^x, ^y, 0, 0) 
are always positive, because, by definition of Gg, they are ensemble averages 
of quantities under square modulus. 

Let us now introduce the Fourier transform of Eq. (196) with respect to Afx,y'- 



^o{f,,rv,0,,Oy)= J dAf', J dAf'y Go{f,,fy,Af'^,Af'y) 

— oo — oo 

X exp[2i(^^Af; + eyAr'y)] . (198) 

Accounting for Eq. (196), i.e. in the particular case of a (virtual) quasi- 
homogeneous source, Eq. (198) can be written as 



$o(fx, fy, 9„ 9y) = i (f,, fy) f(4, 9y) , (l99) 

having recognized that r{6x,0y) = ^o(^a;, ^y, 0, 0) is the Fourier transform of 
the spectral degree of coherence g. The distribution $0, being the product of 
two positive quantities, never assumes negative values. Therefore it may always 
be interpreted as a phase space distribution . This analysis shows that quasi- 
homogeneous sources can always be characterized in terms of Geometrical 
Optics. It also shows that, in this particular case, the coordinates in the phase 
space, fx,y and O^^y, are separable. 

Eq. (198) is the definition of a Wigner distribution. In the case of quasi- 
homogenous sources, as we have just seen, the Wigner distribution is never 
negative and, therefore, can always be interpreted as a phase space distri- 
bution. In the case of non quasi-homogeneous sources one may still define a 

It should be remarked that this result has been obtained only on the ground of 
mathematical basis, i.e. without ascribing to I and T any physical meaning. In other 
words, we simply considered the cross-spectral density Go as the ensemble-averaged 
product of E{ri) and E*{f2) without ascribing to the function E any physical 
meaning. Physically, as has been said in the Introduction, in the quasi-homogeneous 
case r can be identified with the radiant intensity of the virtual source (compare 
with Eq. (6)). This follows from a statement similar to the van Cittert-Zernike 
theorem for quasi- homogeneous sources (see [3]). Note that the intensity and the 
Fourier transform of the spectral degree of coherence are obtained back from the 
phase space distribution, Eq. (199), by integration over coordinates 9x,y and fx,y 
respectively. 
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Wigner distribution using Eq. (198). The integral of the Wigner function over 
its coordinates must still be finite . However the Wigner function itself is 
not always a positive function. As a consequence it cannot always be inter- 
preted as a phase space distribution. On the one hand, quasi-homogeneity is 
a sufficient condition for the Geometrical Optics approach to be possibly used 
in the representation of the source. On the other hand though, necessary and 
sufficient conditions for $o to be a positive function are more difficult to find. 

One may observe that Bochner's theorem may be used to investigate whether 
the Wigner function can be interpreted as a phase space distribution in the 
case of non-homogeneous sources. In particular, it is necessary and sufficient 
to look for non-negative definite cross-spectral density functions. However, in 
general, it is not trivial to investigate wether a function is non-negative definite 
(see footnote 14) and therefore this observation does not seem to constitute 
a simplification. We will simply leave the search for necessary and sufficient 
condition for $o to be a positive function as an open question. We did not rule 
out, for undulator setups, the possibility of having a positive Wigner distribu- 
tion in the non quasi-homogeneous case. At first glance it may look like such a 



This ensures that $o has finite integral over its variables. $o is the Fourier trans- 
form of a correlation function (the cross-spectral density) of electromagnetic fields. 
The fields being physical quantities can carry only a finite amount of energy and 
they are limited in spatial extent. As a result the cross-spectral density must be an 
integrable function over its variables and so must be, by definition, 

Bochner's theorem " in its elementary form asserts that every non- negative defi- 
nite function of a broad class has a non-negative Fourier transform and, conversely, 
that the Fourier transform of every non-negative function of a broad class is non- 
negative definite. This class includes functions which fall off sufficiently rapidly to 
infinity to ensure that their Fourier transforms are continuous functions" [cited 
from paragraph 1.4.2. of reference [3]]. It should hereby be stressed the difference 
in the mathematical language between a positive function f{^) ^0 for every real 
value ^ and a non-negative definite function. The function h is said to be non- 
negative definite when "for an arbitrary set of N real numbers ^i, S,2,---,S,n and N 
arbitrary complex numbers ai, 02, a^r, YliLiYlf=iK^j^i^j ~ ^i) ^ 0." [cited 
from paragraph 1.4.2. of reference [3]]. Based on the assumption of a quasi- 
homogeneous source authors of [3] use Bochner's theorem instead of our previous 
discussion to demonstrate that the spectral degree of coherence g in Eq. (196) 
is necessarily non-negative definite. In fact, since the intensity / is a positive 
function, the sign of J2i^=iJ2f=i^i^j9{^'f^x iji^'^^y ij) is the same of the sign of 
E£i EJLi a*a,ji{f,^ ij,fy ij)g{Af:, ij, Avy „•) = YlZi E^Li < {E{rj)E*{fi)) = 



^ 0, quantum erat demonstrandum. This demonstration is more 



involved than ours, even though it is based, as the ours, on the positivity of the 
square modulus of quantities. The reason for this complexity is that it uses a more 
general theorem, i.e. Bochner's theorem. However, as we have seen, it is not neces- 
sary to invoke Bochner's theorem in the quasi-homogeneous case. 
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case brings advantages in the formulation of the imaging theory, because the 
source can be described in terms of Geometrical Optics. On the one hand, in 
the case the Wigner distribution is positive, the evolution of the radiation in 
free space can be described by a ray-tracing approach, as the Wigner distri- 
bution can be interpreted as a phase space distribution. On the other hand 
though, such fact is almost irrelevant because it is not of help when optical 
elements are considered. As we will see later on, there are two particular con- 
ditions at the basis of a simplified formulation of the imaging theory based 
on the incoherent point spread function of the optical system. The first is the 
separability of the cross-spectral density is the product of a factor depending 
on f and a factor depending on Ar. The second is a transverse dimension of 
the source much larger than the transverse coherence length. As we have al- 
ready seen these two conditions, together, define a quasi-homogeneous source. 
Quasi-homogeneous sources are necessarily characterized by a positive Wigner 
function. However, the positivity of the Wigner function alone is not sufficient 
to obtain a simplified formulation of the imaging theory in terms of inco- 
herent point-spread function. A similar remark holds for a particular kind of 
sources often considered in literature also in connection with undulator radi- 
ation (see [14, 15]). These sources, characterized by a cross-spectral density 
G- 



I I{ri)\J I{r2)g{ri — r2) are called Shell sources (in particular, in [14, 15] 
Gaussian-Shell sources are discussed which assume gaussian a profile for both 
/ and g). They exhibit separability of the cross-spectral density, but are not 
quasi-homogeneous because the transverse dimension of the source fails to 
be much larger than the transverse coherence length: a simplified formula- 
tion of the imaging theory does not hold in this case either. Moreover, as it 
will be more extensively discussed here below and in Section 15, the Shell 
model (and, in particular, the Gaussian- Shell model) may be useful for de- 
scribing light sources other than undulator-based or for educational purposes, 
but does not describe any practical realization of an undulator source. 



We are now interested to find, in particular, equivalent conditions for quasi- 
homogeneity in terms of the electron beam sizes Nx,y and divergences D^^y that 
apply to our case of interest, i.e. third generation light sources. In order to do 
so we start deriving an expression Q for the Fourier transform of the cross- 
spectral density at the virtual source position. This is given by calculating the 
limit of Eq. (124) for Zo^l and taking advantage of Eq. (127). Aside for an 
inessential multiplicative constant we obtain: 



^?(0,^„^^,A^„A^^) =exp -2A^,A^^ 

(4 + ^x) 



X exp 



exp 
- \ 2- 



2D. 



exp 



2D,, 



oo — oo 
2" 
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xsmc 



"(4-A^;)2+(0,-A^;)^ 



xsmc 



(4+Ai)2+(0,+A^;)^ 



(200) 



We have said that the quasi-homogeneity of the virtual source is equivalent to 
(i) separability of the cross-spectral density G in the product of two factors 
respectively depending on O^^y and AO^^y and (ii) a large characteristic scale 
of dx,y with respect to the characteristic scale of 1^9 x.y From condition (i) 
follows that the virtual source is quasi-homogeneous only if it is possible to 
factorize the Fourier transform of the cross-spectral density, Q in Eq. (200), 
in the product of two factors separately depending on 9x,y and AO^^y. Such 
factorization, for third generation light sources, is equivalent to a particular 
choice of the region of parameters for the electron beam: ^ 1, Dx ^ 1 and 
either (or both) 1 and ^ 1 . In this case, the second condition (ii) 

is automatically verified as one can verify inspecting Eq. (200). 

An intuitive picture in the real space is given by a (virtual) quasi-homogeneous 
source with characteristic (normalized) square sizes max(A'^j;j^, 1) and charac- 
teristic (normalized) correlation length square of order Dx^y, 1). As al- 
ready remarked before, in the quasi-homogeneous situation the horizontal and 
the vertical directions can be treated separately, because Eq. (200) factorizes 
in the product of factors separately depending on the horizontal and on the 
vertical coordinates. This corresponds to a large number of independently 
radiating sources given by the product 

Mx,y = max(7Vj;,2/, 1) max(L'j;,2/, 1) ■ (201) 



The number M^^y is, in other words, an estimation of the number of coherent 
modes in the horizontal and in the vertical direction . The number M~l 

It should be remarked here, that these conditions describe the totality of third 

generation quasi- homogeneous sources. In fact, while a purely mathematical analysis 
indicates that factorization of Eq. (200) is equivalent to more generic conditions 
{Nx 3> 1 and iV^ 1, or Dx 3> 1 and Dy 3> 1), comparison with third generation 
source parameters reduces such conditions to the already mentioned ones. 

This is in agreement with an intuitive picture where the photon-beam phase space 
reproduces the electron-beam phase space up to the limit imposed by the intrinsic 
diffraction of undulator radiation. Imagine to start from a situation with Nx^y 3> 1 
and Dx^y 3> 1 and to "squeeze" the electron-beam phase space by diminishing Nx^y 
and Dx,y On the one hand the characteristic sizes of the phase space of the electron 
beam are always of order Nx^y and Dx,y On the other hand the characteristic sizes 
of the phase space of the photon beam are of order max(Nx^y, 1) and max{Dx^y, 1): 
diffraction effects limit the "squeezing" of the phase space of the photon beam. 
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is the accuracy of Geometrical Optics results compared with Statistical Op- 
tics results or, better, the accuracy of the quasi-homogeneous assumption. It 
should be noted that, as M^^y approaches unity, the accuracy of the quasi- 
homogeneous assumption becomes worse and worse and M^^y cannot be taken 
anymore as a meaningful estimation of the number of modes: it should be 
replaced by a more accurate concept based on Statistical Optics. To complete 
the previous statement we should add that Mx,y completely loses the mean- 
ing of "number of modes" when Geometrical Optics cannot be applied. For 
instance when both Ny and Dy are of order unity (or smaller), one can state 
that the Geometrical Optics approach fails in the vertical direction because 
the phase space area is getting near to the uncertainty limit. In this case it is 
not possible to ascribe the meaning of "number of modes" to the number My 
simply because the Geometrical Optics approach in the vertical direction fails. 
However, when Ny and Dy are of order unity (or smaller), but both iVa, 3> 1 
and ^ 1, the cross-spectral density admits factorization in the horizontal 
and in the vertical direction and the source in the horizontal direction can be 
still described, independently, with the help of Geometrical Optics. 

Up to now we discussed about the roles of Geometrical and Statistical optics 
in the characterization of the source only. However, as already remarked, the 
specification of the source constitutes only part of the solution of the imaging 
problem. One has, in fact, to track information regarding the source through 
the optical beamline up to the observation plane. Depending on the situation 
Geometrical Optics may be used or not. For instance, a quasi-homogeneous 
source may well be described in terms of a phase space distribution, but if 
diffraction effects dominate the photon beam transport to the observation 
plane, one cannot use ray-tracing techniques to calculate the intensity profile 
at the observation plane. However, as we will see in the next Section, if the 
virtual source is quasi-homogeneous, the intensity at the observation plane can 
always be expressed as a convolution product between the impulse response 
of the optical system and the intensity which would be recovered at the obser- 
vation plane in the case of an ideal optical system (i.e. one with no aberration 
and non-limiting pupil apertures). In this case, the entire line may be studied 
with the help of ray-tracing programs if and only if the impulse response of 
the system can be recovered by means of Geometrical Optics techniques. 

In Geometrical Optics, a Hamiltonian description of the optical system holds 
so that interaction with optical media (i.e. the system evolution) is conve- 
niently modelled in terms of symplectic transformations. A given symplectic 
transformation S acts on point = {Tox^Sox^^oy^doy) of the phase space $0 
at 5o = and maps it to a point p = {f^, 0^, fy, 6y), of the phase space $5 at 
Zo — z according to 

P = S{~Po) . (202) 
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The phase space distribution is therefore transformed according to 



(203) 



According to Liouville's theorem, the area of the phase space is conserved 
during this process. In the particular case of hnear transformations, one can 
use a matrix formahsm. N successive hnear transformations are represented 
by N matrices Li ... Ln and the resulting transformation is represented by N 
successive matrix multiplications, which give the matrix L = Lk • ... ■ Li. The 
action of L on an element of the phase space is then naturally represented by 
multiplication. The variables p in phase space characterize a ray with a certain 
direction and offset with respect to the optical axis. The task of calculating 
the phase space distribution after a given number of optical elements through 



(p) = l>c 



(204) 



where t indicates transposition, or through the more general Eq. (203) can 
be solved by ray-tracing programs. Once $| is known, these codes usually 
integrate it over the variable 6x and 9y to give the intensity distribution 



oo oo 



Hz^r^^ry) = j j dOy (r-x,0:c,ry,9y) . (205) 



-oo — oo 



However, the same programs may also be used to calculate the Fourier trans- 
form of the spectral degree of coherence through 



oo oo 



r[z,0x,0y)^ J dfx J dry^s(fx,ex,fy,0y)) . (206) 



-OO — oo 



In particular, in free space, Eq. (204) becomes 

(fx, Ox, fy, Oy) = i (O, fx - zOx, fy " zOy^ V (O, O^, By) , 

while Eq. (205) and Eq. (206) reduce to convolutions: 



(207) 



oo oo 

i (z, fx, fy) ^ j dOx j dOy i (O, fx - zdx, fy - zdy) f (O, 4, Oy) (208) 

— oo — oo 

and 
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t (^z,9^,9yj = j dr^ j dry i (o,r:c - z9:c,ry - z9y^T (o,9^,9yj 

— oo — oo 

^r(0,9,,9y) . (209) 

Note that T calculated at 5 = has direct physical sense as the intensity 
distribution in the far zone, i.e. the angular spectrum. Then, Eq. (209) tells 
that, at arbitrary distance z, the angular spectrum does not vary. 

The intensity recovered at the image plane in the case of an ideal optical system 
is a scaled copy of that at the virtual source, regardless of the source. Generally, 
although as we will see exceptions apply, such correspondence between the 
intensity of the source and the observed intensity is only true in the case 
the observation plane is the image plane. In the case of a quasi-homogeneous 
virtual source. Geometrical Optics as well as Statistical Optics techniques can 
be employed to recover the intensity at the observation plane. Results from 
the Geometrical Optics and from the Statistical Optics approach must then 
coincide. Let us prove this fact considering the particular case N ^ 1 and 
£> » 1 in a given direction and showing that we are able to recover Eq. (151) 
by means of Geometrical Optics techniques, namely by means of the matrix 
formalism employed in ray-tracing codes. 

In this particular situation, the photon beam can be modelled as if a Gaussian 
photon beam was present at i = with the same horizontal phase space of 
the electron beam. This is an ansatz on the virtual quasi-homogeneous source 
based on the phase space picture described above since strictly speaking it does 
not make sense to talk about a Gaussian photon beam inside the undulator, 
i.e. within the radiation formation length. If, however, this ansatz is made, we 
can describe the optical equivalent of the Twiss matrix at z — 0. Let us first 
introduce the notion of normalized Twiss parameters as: 



ax = otT , 

lT = Lu,lT , 

e ={uj/c)e, 



(210) 



where Ot, and 77- are the Twiss parameters and e is the emittance per- 
taining the photon beam . In the case under study they are identical to the 
analogous electron beam parameters. We have 



^ Pt{0) -ariO) 
y-ariO) 7t(0) 




(211) 



The Twiss parameters are the second moments of the phase space distribution of 
the photons divided by the emittance. 
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For this exemplification we will assume a non-limiting pupil aperture. Then, 
the linear transformation mapping a phase-space point in f = to a phase- 
space point in z = Z2 is represented by the matrix L. In our particular case of 
interest we have 




\^ Z2/{zt - Z1Z2) Zi/{Zi - Z2) j 




(212) 



As one can see from Eq. (42), L describes a free-space fiight followed by a 
focusing element and a second free-space fiight. A point (/o, ^o) of the photon 
beam phase space at i = is transformed, at i = ^2 into 




L 



Vo 



(213) 



while the Twiss parameters for the photon beam aX, z = Z2 are described by 
the matrix 



=Z2 \z=\j 



(214) 



with 



eaT{z2) 



(i + D) 

Ab+(A + Df Q^zl 



i + (i + id) Qh] [d+{A + b) Qzi 



Z2 } , 



(215) 



e/3r(^2) 



Ab+{A + DfQ^zl 



{Zl - Z2f 



A + D 



(216) 



and 
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AD 



eiT{z2) = - +iA + D) Q^zl + 
A + D ^ ' 



(a + D) z\ 2 (i + L>) Qzi 



(^1 - i-i)' 



+ 



Z\ — Z2 



(217) 



It should be recalled that parameters Q and A have been defined in Eq. (148) 
and Eq. (149). 



The photon phase space distribution at i = ^2 is described by 



27re 



exp 



1t{z-2)11 + 1aT{z2)kfii + PT{z2)fil 
2e 



(218) 



The relative intensity is derived from Eq. (218) integrating over the 771-coordinate, 
which gives 



I{z2) = exp 



2/2 



2e(3T{z2) 



(219) 



Finally, substitution of the expression for 6/3^(^2) obtained in Eq. (216), yields 
back Eq. (151), as it should be. Similar conclusions may be obtained for the 
spectral degree of coherence integrating over the /i-coordinate and inverse 
Fourier transforming the result. 

In spite of these results, we should stress again that Statistical Optics is the 
only mean to deal with the stochastic nature of Synchrotron Radiation in 
general. Only in particular cases Synchrotron Radiation can be treated in 
terms of Geometrical Optics. As we have just discussed, one of these cases is 
constituted by second generation fight sources, when Nx,y 3> 1 and D^^y ^ 1. 
Experiments in this region of parameters can take advantage of ray-tracing 
code techniques. 

To conclude this Section, we would like to make a much stronger statement: 
there are practical cases of interest when the Statistical Optics approach must 
be used even for second generation fight sources. This should not sound too 
awkward since, as we have stated before, the impulse response of an optical line 
may not be treatable in terms of Geometrical Optics. Gonsider, for instance, 
the setup illustrated in Fig. 17. This consists of an entrance slit, a grating and 
an exit slit, that is a grating monochromator. The grating equation, which 
describes how the monochromator works, relies on the principle of interference 
applied to the light coming from adjacent grooves. Such principle though, can 
only be applied when phase and amplitude variations of the electromagnetic 
field are well defined across the grating, that is when the field is perfectly 
transversely coherent. If the transversely coherent spot of the radiation is 
smaller than the grating, not all the grating is taken advantage of, resulting 
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entrance 
slit 



grating 



I 




exit 
slit 




Fig. 17. Illustration of a grating monochromator. 
in a decrease of resolution in wavelength. To better explain this point, with 



reference to Fig. 17, let us indicate the width of the entrance slit with d and 
the angle of incidence of the incoming radiation with a. Moreover, let Dg be 
the typical dimension of the grating, Ng the total groove number, m the order 
of diffraction and z the distance between the entrance slit and the grating. 
The maximal (relative) resolution which can be obtained with a particular 
grating is given by {mNg)~^ . Qualitatively, to obtain such maximal resolution 
we must have a transverse coherence area of at least the size Dg cos(q;). If it is 
smaller, not all the grating is used. We now need to transform this qualitative 
requirement into a quantitative requirement. 

In practical situations, the grating is placed in the far zone with respect to 
the entrance slit. This is because the radiation spot size at the grating should 
be at least of order DgCos^a), which is much larger than the slit aperture d. 
If we assume the slit uniformly illuminated, we can consider the slit itself like 
a quasi-homogeneous source with rectangular profile. Then, the van Cittert- 
Zernike theorem applies at the grating position in the far zone. As a result, 
the modulus of the spectral degree of coherence l^^l in the far zone is equal 
to the modulus of the Fourier transform of the intensity profile at the slit, 
which is a rectangular profile. The following expression for \g\ is found in the 
dispersion direction: 



A quantitative requirement for the coherence property of the radiation at the 
grating can be given imposing that \g\ varies within a fixed interval. For in- 
stance, one may require 0.8 < |(?(Ar)| < 1. This requirement may be changed, 
and is somewhat subjective. However it corresponds to a quantification of the 
transverse coherence properties of the radiation on the grating. In particular, 
if the criterium 0.8 < |5f(Ar)| < 1 is chosen, the argument of the sine func- 




(220) 



84 



tion in Eq. (220) is allowed to vary a certain range [— X, X] with X ~ 1.13. 
Moreover, the maximal value for Ar is DgCos{a), i.e. the grating dimension. 
Putting all together we obtain the following condition for the monochromator 
setup parameters: 

1I^£^^ ^X^ 1.13 . (221) 
Xz 



The same calculation may be repeated with a different choice for the minimal 
allowed value of |g'(Ar)|. This would lead to a different value of X. 

Our result is in agreement with the conclusion that one may draw considering 
the following relation [16]: 

AA ^ (dD,cos{a)\ J_ 

\ \ Xz ) mNg ^ ' 



Eq. (222) describes how the entrance slit width hmits the resolution in wave- 
length according to. The second factor on the right hand side of Eq. (222), 
i.e. {mNy)^^ is, again, the maximal relative resolution. This resolution can 
be obtained by setting the first factor to unity. This yields a result in para- 
metrical agreement with Eq. (221). The right parametric dependence in the 
condition for the maximal resolution can also be obtained in another way. If 
the radiation sent through the slit has, at the grating, a spot size equal to 
DgCos{a), the condition for transverse coherence is given by the space-angle 
product: 

^^^dDgCOsja) (223) 
2z 2n ^ ^ 



Again, the qiialitative estimation in Eq. (223) is in parametrical agreement 
with the quantitative calculation in Eq. (220). 

The grating works with a resolution near to the theoretical limit (miVg)~^ only 
with transversely coherent radiation. We may say that the purpose of the en- 
trance slit is to supply a transversely coherent radiation spot at the grating 
in order to allow the monochromator to work with a certain resolution. This 
fact must hold for any light source, and in particular for second generation 
light sources. The bottom line is that this monochromator setup cannot be de- 
scribed in terms of Geometrical Optics even in the case of a second generation 
light source: transversely coherent radiation means that the image on the exit 
slit is close to the diffraction limit. Therefore, in this case. Geometrical Optics 
can only be used for approximate estimations, while a correct treatment must 
involve the application of Statistical Optics techniques. 
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8 Imaging of quasi-homogeneous undulator sources: effect of aper- 
ture size 



We will now consider, with the help of Eq. (135), the effects of a pupil in 
the one-dimensional case (that can be practically realized with the help of a 
slit aperture and a cylindrical lens) when aberrations are not present. First, 
in Section 8.1, we will consider Gaussian quasi- homogeneous sources. Then, 
in Section 8.2 we will see that the arguments for Gaussian sources can be 
generalized without modifying their substance with the help of some notational 
change to treat the case of non-Gaussian quasi-homogeneous sources as well. 



8.1 Quasi-homogeneous Gaussian undulator sources 



In the case for ^ 1 and S> 1 we can treat the horizontal and the vertical 
direction separately. Then, also the function P and V can be expressed as the 
product of factors in the horizontal and in the vertical direction. In particular 
we may consider the pupil function 



P(r) 



1 




if |f| < d 
otherwise 



(224) 



where f may represent either the variable x or the variable y. According to 
the definition in Eq. (101) this gives: 



V{u) = 2asinc(a'u) . 



(225) 



We can use Eq. (135) and Eq. (156) to describe the case when the lens is in 
the far zone, that is when condition (136) is satisfied. From Eq. (146) we can 

estimate the characteristic size of the source, that is of order and of the 

correlation length at the source, that is of order 1 / \/^. According to condition 

(136), the lens is in the far zone when / Zi <^ \/^, in agreement with the 
far zone limit of Eq. (142), which is obtained for A <^ D . In this limit, Eq. 
(135) and Eq. (156) give 

Condition (136) is usually not discussed in textbooks describing thermal sources. 
In fact, for perfectly incoherent thermal sources, the far zone is defined by z ^ 6"o, 
(Jo being the source transverse size, i.e. when the paraxial approximation is valid. 
Therefore, in this case, the pupil is always in the far zone. 
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Gp{zi,ri, Afi) =4a^exp 
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2im 



r,;Af,; 



X du dAu exp 



2N 



exp 



-2Dzl{Auf 



xsmc < a 



xsmc < a 



m 



— (vi + Afj) — u — Ait 



m 



— (rj - Afi) - + Axi 
^1 



(226) 



According to the far field fimit of Eq. (142), the quantity Dzf is the square 
of the radiation spot size on the pupil, while zf/N is essentially the square of 
the coherence length on the pupil. Two interesting limiting cases of Eq. (226) 
can be obtained comparing these two characteristic scales with a^, that is the 
square of the pupil size. 

First, let us consider the case zf/N < <C Dzf. As we will demonstrate 
later on, in all situations when the quasi-homogeneous assumption is verified, 
the exponential function in Aii inside the integral in Eq. (226) behaves like a 
5-Dirac distribution, and one obtains 



Gp{zi, ri, Afi) = 4d^ exp 



xsmc < a 



xsmc < a 



2im 

Zl 

m 

Zl 

m 



-uAfi 



Zl 



du exp 

fi + Afi) — u 
'fi - Afi) - u 



zfu^ 
2N 



(227) 



This corresponds to a relative intensity 



Ip{zi,ri) = ^ I du exp 



zju'^ 
2N 



smc 



a \ — — u 
\di I 



(228) 



where the normalization constant C is given by 



C— J d'U' 6xp 



zfu'^ 
2N 



|sinc (du) 



(229) 



Eq. (228) expresses the image as a convolution of the geometrical image with 
the slit diffraction pattern (in two dimensions this would be the Airy pattern). 
It is valid for values zf/N < d^ <^ Dzf and also for zj / N <ti d,'^ <^ Dzf: the 
difference between these two cases is that in the first the pupil influence is 
significant, while in the second it is not. Inspection of Eq. (228) and use 
of Eq. (31) yields the ratio between the size of the diffraction pattern and 
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the geometrical image: y\M\y N a/ dij = \^ayN/zij . When zj/N < 

such ratio is comparable with unity, i.e. the diffraction pattern significantly 
influences the image formation process. When, zf/N <^ a^, this ratio is much 
smaller than unity. As a result, the pupil influence is not significant and the 
image is given by the geometrical image. The ratio between the size of the 

diffraction pattern and the geometrical image, Zi/a\/^, gives the resolution 
of the image due to diffraction effects. It is also interesting to note that in 
the limiting case when zl/N <^ a? <^ Dzf, Eq. (227) presents the asymptotic 
behavior 



Gp{zi, ri, Afi) = Aa? exp 



2im 



exp 



2 —2 

m rf 



2N 



X J du sine i d 

— oo 



in . ^ 

— Arj — u 



smc < a 



— Arj + u 



(230) 



The convolution theorem yields the following expression for the spectral degree 
of coherence: 



gp{zi,Afi) = ^ I df'|P(r')|'exp 



-2i 



f'Afj 



(231) 



the normalization factor T> being given by 



oo 

r>= df'\p(f') 



(232) 



After the substitution u — ^ u'y/N/zi we may rewrite Eq. (228) as a function 
of ^ = dmfi/zi with the help of the only parameter p — dy/N / zi, that is easier 
to plot 



-^p(C) = ^ I du' exp 



u 



12 



|sinc ({ -Pm)|^ , 



(233) 



where C can explicitly be calculated as: 



F 1 
2^ 



(-1 + exp \-2p^\ I + - erf [V2p 

L L JJ p L . 



(234) 
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Fig. 18. Image intensity for a quasi-homogeneous source, Ip, as a function of 
^ = amfi/zi, calculated with Eq. (233), for different values of the parameter 

p = a\/^/zi. The plot illustrates the one-dimensional image formation problem 
(slit aperture, cylindrical lens). 

The function erf(-) in Eq. (234) indicates the error function. Note that the 

variable ^ = amfi/zi may also be written as ^ = {a/di)fi, where di = Zi — Zi 
and |M| = di/zi = m~^. It is interesting to remark that the ratio di/d is 
the dimensionless characteristic size of the Fresnel zone, i.e. \di/{2'Ka). Ip is 
plotted in Fig. 18 for several values of the parameter p. In Fig. 18 we also plot 
the asymptotic behaviors of Ip for small values of p, i.e. Ip — sinc^(^), that is 
the diffraction pattern from a slit, and for large values of p, i.e. exp[— (^/p)^/2]. 

On the one hand, the integral in du in Eq. (228) is a convolution. On the 
other hand, aside for numerical factors, the Fourier transform of \V{u)\'^ ~ 
|2asinc(a'u)p can be given in terms of the triangular function tri(-), defined 

IQ 

as ^ 



It should be noted that x = 2Af. The reason why we introdTiccd the new variable 
X is to keep a certain homogeneity of notation when comparing with reference [8]. 
Since our definition of Af differs for a factor 2 with respect to that in [8], it is 
somewhat convenient to introduce % = 2Af. 
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Fig. 19. Upper plots: profile of the pupil function P (slit aperture), together with 
V. Lower plots: autocorrelation function of the pupil, P, together with J^iV)- 
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l + x/(2a) if-2a<x<0 
l-x/(2a) if 0<x<2a 
otherwise . 



By means of the convolution theorem, Eq. (228) can be written as 



— oo 

The normahzation factor J\f is given by 



tri(|^,exp 



m 



J dx 
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2zl 
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(235) 



(236) 



'2N d 



(237) 
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Eq. (236) can be written as an analogous of Eq. (5.7-10) of [8], that sometimes 
goes under the name of Shell's theorem. Let us define the autocorrelation 
function of the pupil as 



oo 



(238) 



In Fig. 19 we plot the profile of the pupil function P, together with V and the 
autocorrelation function of the pupil P, together with its Fourier transform 
T{P). With the help of Eq. (238), Eq. (236) can be written as 



1 °° 



2zl 



P{X) exp 



m _ 



(239) 



that is Eq. (5.7-10) of reference [8]. 

A second limiting case of interest is found when <C zj/N <^ Dzf. The 
pupil is coherently and uniformly illuminated. In this case both the exponential 
functions inside the integral in Eq. (226) behave like 5-Dirac functions yielding 



Gp{zi,ri,Ari) =4a exp 
xsinc < a 



2im 



smc < a 



Zl 

"m 

— (rj - Afi 

-Zl 
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Zl 



(240) 



corresponding to a relative intensity 



Ip{zi,ri)^smc'^ 



am 



Zl 



1 



V ( —fi 

Zl 



(241) 



Also, using the definition of the spectral degree of coherence given in Eq. (129) 
one sees that \g\ = 1, i.e. the pupil is coherently illuminated. Eq. (241) is the 
analogous of Eq. (5.7-14) in [8]. 

Results obtained here deal both with the cross-spectral density Gp and the 
relative intensity Ip in the presence of the pupil. We classified these results 
comparing the square of the extent of the pupil with the square of the 
radiation spot size Dzf and of the coherence length zf/N at the pupil location. 

Let us now consider the intensity Ip only. In the case of quasi-homogeneous 
sources, there is a more general method taken from the theory of linear systems 
to account for the pupil presence. This method can be extended to account 
for aberrations as well, and allows a more compact treatment of the pupil 
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effects because it does not depend on how the extent of the pupil scales with 
respect to the coherence length and to the radiation spot size. It is based 
on the concept of line spread function. A linear time-invariant system in two 
dimensions can be characterized by the knowledge of the point spread function 
(or impulse response) h{x,y). Given a certain input f{xi,yi), the output at 
any point {x, y) is given by the convolution of the point-spread function h and 
the input /. The point spread function h{x, y) is the response to a 5-Dirac 
signal at position (0,0), i.e. 6{xi,yi). The line spread l{x) is, instead, the 
response obtained from a line input S{xi), which is independent of yi, and can 
be calculated by integrating the point spread function h{x, y) with respect to 
the y- variable (see, for instance, [17] Sec 6.2.). 

In the quasi- homogeneous case, the intensity Ip can be written as a convo- 
lution between a suitable line spread function I ^° and the intensity / which 
does not account for the pupil presence, that is 

ip{zi,fi) = [i*i]{z,,n) (242) 



The line spread function / acts as the passport of the imaging system, and 
depends on the properties of the lens only. In this case we are only considering 
the effect of a finite pupil dimension, i.e. we are accounting for diffraction 
effects from the pupil. More in general, I may depend on lens apodization or 
aberrations too. 

In the case under study here, the line spread function of the system is given 

by 



l{fi) ^ smc" [tf}j . (243) 

For instance, it is straightforward to see that substitution of the (magnified) 
input signal / = exp[— ■u^/(2|MpiV)] with the input signal / = 5(-u/|M|), that 
is a line input, in Eq. (228) gives back Eq. (243). 

It should be emphasized that the resolution due to diffraction effects is of order 
Zi/{a\[^). In all cases when this resolution is better than (i.e. zi/{a\[^) is 

'^^ Here we are treating a two-dimensional system, but we are considering the case 
when we have separability properties for both the source and the pupil. For the 
source this means that Nx ^ 1 and 2> 1, while for the pupil it means that the 
pupil is rectangular. This case is practically realized with the help of slit apertures 
and cylindrical optics. The line spread function I is, then, the proper tool to consider. 
If one wants to consider the situation when no separability property for the pupil 
is present, one should take advantage of an approach based on the point spread 
function. 
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smaller than) the resolution of the ideal image (which does not account for 
the pupil presence) the pupil docs not play any role and, with the accuracy 
of the calculation of the ideal intensity (see Section 15), the / function in Eq. 
(243) cannot be distinguished from a 5-Dirac. In our study case, the ideal 
intensity / = ex.p[—zlu^ / {2N)] is calculated with an accuracy which is much 
worse than the quasi-homogenous accuracy (see Section 15 for details) and 

is of order max(l/\/^, 1/^^). In general, it is important to compare the 
resolution due to diffraction effects with the accuracy of the calculation of 
the ideal intensity. For example, on the basis of such comparison, one may 
conclude that the case for zf/N <^ < Dzf can be calculated assuming that 

the line spread function in Eq. (243) is a (5-Dirac. In fact zi/{a\[^) ~ 1/\J ND 

since a? < IDzl and therefore Zi/{a\[^) <^ max(l/\/^, \/\[^). This kind of 
reasoning can be used to treat any quasi- homogeneous case, and the / function 
can be modified to include aberrations and apodization effects as well. In the 
most general case, accounting for the effects of the pupil in one dimension, 
Eq. (135) yields the relative intensity 
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where the normalization factor S is given by 



du dAu exp 
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(245) 



In the quasi-homogeneous case (including cases when the source is not Gaus- 
sian, see Section 6), if the lens is in the far field, one can write 



G i Zi, —u, —Ail ] ^ I izi, —u ] g [zi, —Au ] exp 
\mm / \ m/ \m / 



2izi 



m 



u ■ Ail 



. (246) 



Further on, within the accuracy of the quasi-homogeneous approximation, the 
spectral degree of coherence g behaves hke a Dirac 5-function in the calculation 
of both the intensity and the cross-spectral density. In fact, the accuracy of the 
incoherent impulse response \V\'^ is also the accuracy of the quasi-homogenous 
assumption, and this is the accuracy with which we can substitute g with a 
Dirac 5-function on the image plane. As a result, in analogy with Eq. (227) 
one has 
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Gp{zi,ri,Ari) = 



—oo 



xV — {fi + Arj) - it V* —{fi- Afi) - it 




(247) 



while the relative intensity can be written as 




(248) 



— oo 



I) being defined in Eq. (232). The line input response is obtained by setting 
I{zi, ziU/m) to a Dirac 5-function, thus obtaining the line spread function: 



The line spread function is therefore the inverse Fourier transform of the au- 
tocorrelation function of the pupil. The autocorrelation function of the pupil, 
i.e. the Fourier transform of the line spread function, is also known as the Op- 
tical Transfer Function (OTF). Other relevant quantities introduced in linear 
system theory are the phase of the Optical Transfer Function and its modulus, 
which is known as the Modulation Transfer Function (MTF) [17] . 

Eq. (248) and, consequently, the line spread function approach, constitutes 
a universal description of the intensity in all quasi-homogeneous cases. In 
literature the line spread function is used to describe perfectly incoherent 
sources only. Note that, in general, radiation produced by an electron beam in 
an undulator is similar to an incoherent sum of many independent laser-like 
beams. Yet, it cannot be considered as an incoherent sum of point sources 
because, as we have seen in Section 3, a single electron cannot be considered 
as a point-like radiation source. Radiation produced by a single electron is 
similar to a laser beam. If no influence of focusing is present in the undulator, 
this laser-like beam has a waist located in the center of the undulator. At 
the waist the radiation wavefront is plane and the radiation spot size is much 
larger than the wavelength. We extended the use of the line spread function 
approach to the realm of quasi-homogeneous sources. 

In the case of third generation light sources the line spread function method 
can almost always be applied in the horizontal direction. However, it fails in 
the vertical direction, where third generation light sources are seldom quasi- 
homogeneous. If the source is not quasi-homogeneous, the cross-spectral den- 
sity cannot be factorized in the product of the intensity and of the spectral 
degree of coherence, or the coherence length is not short (compared with the 
size of the source). As a result, the incoherent line spread function I cannot 



2 



,/ N ^ ni 

l{fi) = V -fi 



(249) 
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be used to describe the system. The function V is known as the coherent 
hne spread function and must be used in its place. In fact, when the source 
starts to exhibit a high degree of transverse coherence (i.e. in the non quasi- 
homogeneous case), the coherent hne spread function, V acts on the field at 
the image plane analogously to the way the incoherent line spread function 
acts on the intensity at the image plane. To see this, it is sufficient to inspect 
Eq. (100). 

It is interesting to compare this viewpoint with what can be found in literature. 
For instance, in [6] , where a condenser system is discussed, one may read: " The 
intrinsic divergence of the extreme ultraviolet (EUV) undulator considered 
here is 6*0611 = SO/^rad, which is larger than the beamline acceptance 6'accept of 
48/irad. Therefore it is evident that the incoherent source approximation holds 
here and the term incoherence source is used accordingly in this paper". In 
the following Sections of their work, authors of [6] use a point-spread function 
approach to accoimt for aberration effects: in their paper, the intensity at 
the source is used, instead of the cross-spectral density, in order to evaluate 
both the intensity and the degree of coherence at the image plane. Such an 
approach is justified in the passage above, where they state that the source is 
incoherent. 

The statement in [6] about the incoherence of the source is a misconception. 
According to such statement, perfectly coherent undulator radiation produced 
by an electron beam with zero emittance should exhibit incoherent properties 
when the radiation divergence is larger than the acceptance of the optical 
system. In contrast with the assertion made in [6], the coherence properties 
of the source are independent of the beamline elements which follow. In order 
to discuss about the coherence properties of the source one has to refer to the 
radiation field at the virtual plane location only. In particular, the fact that 
the source is coherent (or not) does not depend on how the beam acceptance 
angle scales with the intrinsic (single particle) divergence of the undulator 
radiation. Our conclusion is that the only parameters which describe wether 
a source is quasi-homogeneous or not are (in the vertical direction) Ny and 
Dy. If the source is quasi-homogeneous, a point-spread function approach can 
be used. If not, the more general results described in Section 14.1 should be 
considered. In particular, in the case of [6], the vertical rms dimension of the 
source is Uy — 16 //m and the radiation wavelength is A = 13.4 nm, while the 
undulator (see [7]) is composed of 55 periods, each one 8 cm long. This means 
Lyj = 4.4 m. Moreover, the vertical emittance at ALS is ey ~ O.lnm, while 
the vertical beta function for beamline 12 is f3y = 4.2 m ~ Lyj. As a result 
both Ny ~ 0.1 and Dy ~ 0.1 and the source is non-homogeneous. We conclude 
that, in this case, approximations like Eq. (247) or Eq. (248) cannot be used. 
Eq. (376) in Section 14 should be considered instead. 
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As a final remark we should stress that, even in cases when the virtual source 
is quasi-homogeneous, one should verify the assumption that the lens is in 
the far zone, before applying a point spread function formalism. In contrast 
to this, note that in the usual framework of Statistical Optics, the radiant 
intensity from thermal sources is distributed over a solid angle of order 27r, 
and optical elements can always be considered in the far zone. 



8.2 Quasi-homogeneous non-Gaussian undulator sources 



In the present Section 8.2 we will extend results obtained in Section 6. Results 
obtained in the previous Section 8.1 apply for a quasi- homogeneous Gaussian 
undulator source only. In particular, under the assumptions ^ 1 and 
Dx ^ 1 the cross-spectral density can be factorized in a horizontal and in a 
vertical contribution, and results in Section 8.1 can be applied in the horizontal 
direction. Note that, \l 1 and 1, one has, automatically, ^ 

and £>-r » 1 and the same results in Section 8.1 can be separately apphed 
in both the horizontal and the vertical directions. Here, with the help of Eq. 
(135), we will include the effects of a pupil in the one-dimensional case when 
the source is still quasi-homogeneous, but non-Gaussian. In particular, we 
will still assume ^ ^ ^-nd ^ ^ and concentrate our attention on 
the vertical direction. First we will study the case when Ny ^ 1 and Dy 
is arbitrary and, then, the case when Ny is arbitrary and Dy ^ 1. We will 
see that the reasoning applied in the case of quasi-homogeneous Gaussian 
sources also holds in the case for quasi-homogeneous non-Gaussian sources 
as it relies on the separabihty of the cross-spectral density only. As a result 
we will present practical examples of how, with minor substitutions, we can 
extend our analysis of the pupil effects to the case of non-Gaussian sources. 



8.2.1 Source with non- Gaussian angular distribution in the vertical direction 

Let us start considering the case when Dy is arbitrary and Ny ^ 1. The pupil 
function and V are given by Eq. (224) and Eq. (225). The r-direction should 
be now substituted with the y-direction. 

We can use Eq. (135) and an asymptotic expression of Eq. (168) to describe 
the case when the lens is in the far zone, that is when condition (136) is 
satisfied. Prom Eq. (162) we can estimate the typical size of the source that 

is of order ^f^, and of the correlation length at the source, that is of or- 
der min[l/y^, 1]. According to Eq. (136), the lens is in the far zone when 
\[^y/zx < max[Y^, 1]. In this limit, Eq. (135) and Eq. (168) give 
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Prom Eq. (159), one sees that zl max[£)j,, 1] is of the order of the square of the 
radiation spot size on the pupil, while z\lNy is of the order of the square of 
the coherence length on the pupil. Note that the limiting expression obtained 
from Eq. (250) for Dy ^ 1 is Eq. (226). As before, two interesting limiting 
cases of Eq. (226) can be obtained comparing these two scales with a^, that 
is the square of the pupil size. 

First, let us consider the case z\jNy < <^ zf max[l, Dy]. As we have already 
discussed, in all situations when the quasi-homogeneous assumption is verified, 
the exponential function in Ait inside the integral in Eq. (250) behaves like a 
5-Dirac distribution. As in Eq. (227) one obtains 
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Setting Ayi — one yields the intensity 
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(252) 



where the normalization constant C has already been defined in Eq. (229). The 
result in Eq. (252) is equivalent to the intensity already given in Eq. (228). It 
can also be written as in Eq. (236), that is 



Ip{zi,yi)^^ I dxexp 
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where J\f is defined in Eq. (237). After introduction of P as in Eq. (238), Eq. 
(253) can be rewritten as Eq. (239), 
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(254) 



In this case, the expression for the intensity is the same as in the case Dy ^ 1. 
Also note that the results obtained for the case zf/Ny ^ <^ zfmax[l, Dy] 
are valid in the limit zf/Ny <^ <^ zlmsix[l, Dy] as well. In this case Eq. 
(252) can be simphfied to 



U{zi,yi)^exp 
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2N 



(255) 



The second limiting case that we can mention here for comparison with what 
has been done in the Gaussian case is when cP' <^ zl/N <^ max[£)j,, l]zl the 
pupil is coherently and uniformly illuminated. In this case one recovers the 
same results in Eq. (240) and Eq. (241). 

To sum up, we obtain, in all situations, the same intensity as in the case 
Dy->1 . 



8.2.2 Source with non-Gaussian intensity distribution in the vertical direc- 
tion 

Let us now study the case when Ny is arbitrary and Dy^l. 

We can use Eq. (135) and an asymptotic expression of Eq. (192) to describe 
the case when the lens is in the far zone, that is when condition (136) is 
satisfied. Prom Eq. (188) we can estimate the typical size of the source that 

is of order maxf^^T^, 1] , and of the correlation length at the source, that is 
of order l/\fDy. According to Eq. (136), the lens is in the far zone when 
max[Y^, l]/zi < y/hy. In this limit, Eq. (135) and Eq. (192) give 



Gp(zi,yi,Ayi)^Ad exp 



2im 

Zl 



oo 

X j d^yBXp 
—oo 



J du dAii 

y + 



2M, 



B(^2/)exp -2byzl{Auy 



xsmc < a 



xsmc < a 



m 



— {yi + A^j) -u- Au 

Zl 



m 



— {yi - ^m) -u + Au 

Zl 



}■ 



(256) 



According to Eq. (171), the quantity zfDy is of the order of the square of the 
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radiation spot size on the pupil, while zl/ max[iVy, 1] is of the order of the 
square of the coherence length on the pupil. Note that the limiting expression 
obtained from Eq. (256) for Dy ^ 1 is Eq. (226). We will study again two 
limiting cases of Eq. (256), which can be obtained comparing these two scales 
with a^, that is the square of the pupil size. 

First, let us consider the case z\/ max[iVj,, 1] ^ ^ zlD. As before, because 
of the quasi-homogeneous assumption is verified, the exponential function in 
Ati inside the integral in Eq. (256) behaves like a 5-Dirac distribution. One 
obtains 
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This corresponds to the intensity 
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where 
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In analogy with Eq. (239), Eq. (258) can also be written as: 
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where the normalization factor J\f is defined as 



99 




(261) 

Note that results obtained in the case zl/ max[iVj^, 1] < <^ zlD are also 
valid in the asymptote for z\/ max[iV'j^, 1] <^ <^ zlD. In this case, Eq. (258) 
is simplified to 
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where 




(263) 



The second limiting case that we will mention here for comparison with what 
has been done in the Gaussian case is for ^ zf/ max(iVj^, 1) <C Dyzf the 
pupil is coherently and uniformly illuminated. In this case one recover the 
same results in Eq. (240) and Eq. (241). 

However, in general, the expression for the intensity is different from that for 
7V„ > 1 . 



9 Aberrations and imaging of quasi-homogeneous sources 



Up to now we have discussed cases when no aberrations are present. Although 
there are widespread treatments of Aberration Theory in literature, we will in- 
troduce our own here, so that this work is self-consistent. In particular, we will 
focus on the one-dimensional case, which has not been treated widely in books 
and monographies, aside for some exception (see [18]). In the present Section 
9 we will assume that the virtual undulator source is quasi-homogeneous. Al- 
though we will begin introducing the Optical Transfer Function (OTP) for 
the system, we will mainly be concerned with the line spread function of the 
system. In addition to that we will discuss the case of severe aberrations, 
presenting new analytical results for this asymptote and comparison with nu- 
merical calculations. Finally, we will intensively discuss defocusing aberrations 
and present general analytical results for this case too. Our particular consider- 
ation of the defocusing case is justified by the fact that this is a privileged kind 
of aberration in the framework of Fourier Optics. In fact it shows a quadratic 
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Fig. 20. The line spread I function in the presence of defocusing aberration. 



dependence of the phase error along the pupil aperture, that is quite a natural 
dependence in the Fourier Optics approach. Quadratic phase factors, which 
can be interpreted as defocusing aberrations, are present even in the case of 
a pinhole camera setup (that will be treated in the next Section 10) and of 
imaging from ideal lenses in an arbitrary plane of interest behind the lens (as 
it will be seen in Section 11 and 12). 



9.1 Aberrations and the effect of aperture size 



In general, aberrations may be accounted for substituting the pupil function 
P, the coherent line spread function V, the autocorrelation function P and the 
Fourier transform of the autocorrelation function J^{P) with new functions, 
respectively PayVa, Pa and J-'{Pa), that account for aberrations. Since from the 
very beginning of the present Section 9 we considered only the one-dimensional 
case, we will consider one-dimensional aberration theory only. Generalization 
is possible, although calculations would become more cumbersome. 

Mathematically, the presence of aberrations in one-dimension modifies the 
pupil function by means of a phase error $(r) leading to 



101 



<l.(f) = b3 (f/a)3 




^ r.a/d. 

10 I I 



Fig. 21. The line spread I function in the presence of coma aberration. 



cD(f) = b4 (f/a)4+b2 (r/a)2 




b4=3/2 n ■ b2=0 



Fig. 22. The hne spread / function in the presence of both spherical and defocusing 
aberrations. 
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Fig. 23. The transfer function in the presence of defocusing aberration. 
P„(r) = P(r) exp [i^ (f)] . (264) 

It is customary (see, for instance, [18]) to consider phase errors of the form: 
^(r) = ^n(^0 ' ^^^^^ 

n being an integer number, and bn being the maximum phase error at the edge 
of the aperture. The value n = corresponds to a constant phase error and 
has no effects, as it cancels out when one calculates the pupil autocorrelation 
function. The value n = 1 contributes to the autocorrelation function for 
a phase term linearly varying with the position. Its only effect is to shift the 
image position. We will not deal with the cases n — and n — 1. We will focus 
instead on the values n — 2 corresponding to defocusing, n — 3 corresponding 
to coma and n — 4 corresponding to spherical aberrations. 

As we have seen in Section 8, the presence of the pupil can be dealt with in- 
dependently of how the pupil extension a scales with respect to the radiation 

spot size and coherence length at the pupil location. Every quasi-homogeneous 
imaging system can be described by means of a line spread function that de- 
pends on the physical characteristics of the pupil only. As has been shown 
in Section 8, the knowledge of the hne spread function I and of the inten- 
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Fig. 24. The transfer function in the presence of coma aberration. 

sity at the image plane without the pupil influence allows one to reconstruct 
the actual image by means of a convolution operation. In the present case of 
interest, the line spread function must be calculated accounting for the pres- 
ence of aberrations. With in mind the purpose of calculating the line spread 
function in presence of aberrations, and with the only assumption of a quasi- 
homogeneous source, we start considering Eq. (137). The quasi-homogeneous 
assumption allows to represent the cross-spectral density at the virtual source 
plane as G{0, —ziu, —ziAu) = /(O, —ziu)g{0, —ziAu). By definition of quasi- 
homogeneity, and with the accuracy of the quasi-homogeneous assumption 
(max[l, Ny] ■ max[l, Z)^])"^/^, the spectral degree of coherence g plays the role 
of a Dirac 5-function in the calculation of the intensity. Such calculation begins 
from Eq. (137). Accounting for the pupil influence, the following expression 
for the intensity at the image plane is therefore found: 



Ip{zi,ri) 



V 



j du 











/(O, -ziu) , 



(266) 



where we recall that V is deflned in Eq. (232). Note that Eq. (266) is obtained 
from Eq. (137) under the only assumption of quasi-homogeneity. Also, when 
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Fig. 25. The transfer function in the presence of both spherical and defocusing 
aberrations. 



|'^aP(0 — ^ ^i') we obtain back an ideal lens and ip{zi,ri) — /(O, — mrj). 

By definition, the line spread function can be obtained by letting /(0,f) 
S{r) in Eq. (266). As a result we have 



/ dx Paix) 



exp 



■ X- 
di 



(267) 



The proof of the last equahty in Eq. (267) is based on the autocorrelation 
theorem, which states that if the (two-dimensional) Fourier transform of a 
function w{x, y) with respect to variables and ay is indicated by w{ax, ay), 
then the Fourier transform of the two-dimensional autocorrelation function of 
w{x, y) with respect to the same variables and ay is given by \w{ax.i oty)\'^- 
In formulas, after definition of the autocorrelation function 



oo oo 



A[w]{x,y)= j dj] j d^w{r] + x,^ + y)w*{r],^) , 



(268) 



— oo — oo 
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which is equivalent to 

oo oo 

A[w]{x,y)^ I dri I d^wir] + x/2,C + y/2)w*iri - x/2,C - y/2) , (269) 

— oo — oo 

the autocorrelation theorem states that 

oo oo 

J dx J dyex.p[i{axX + ayy)]A[w]{x,y) ^ \w{ax,ay)\'^ . (270) 

— oo — oo 



Eq. (267) could have been written down immediately just recalling that the line 
spread function is the Fourier transform of the pupil autocorrelation function. 

It is convenient to make the change of variable x — ^ Axa. In fact, because 
of the definition of 6„ in Eq. (265), only the ratio Ar/a, or x/a = 2Af/a is 
important. Aside for an unimportant factor a, the line spread function becomes 



oo 

l{fi) — j dAx Pa(aAa;) exp 



—i—AxTi 
di 



(271) 



The autocorrelation function of the pupil Pa, which is the Fourier transform 
of I, is known under the name of Optical Transfer Function (OTF) of the sys- 
tem. The Optical Transfer Function and the hne spread function are obviously 
equivalent and their knowledge solves the problem of accounting for aberra- 
tions. The Optical Transfer Function can be written explicitly modifying Eq. 
(238) to include a phase error in the expression for the pupil function, that is 



\{dAx) ^ j dx P 

X exp < i$ 



/ Ax\^ 










P 
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d\x + 



Ax\ 



Ax\ 

a\x 

2 J 



(272) 



Accounting for the definition of the phase error $ in Eq. (265) one obtains 
from Eq. (272) the following expression for the Optical Transfer Function: 




d[x -\- 



Ax 
~2' 

Ax 



a X 



Ax\ 



- \ x- 



Ax 



(273) 
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Here P is the pupil function with no aberrations or apodizations and only an 
aberration term of order n has been considered. If x is outside the interval 
[—1,1], at least one of the P functions in the integral gives zero value. As 
a result we may substitute the integration limits in Eq. (273) to obtain the 
following expression for the Optical Transfer Function: 



1 



dx P 



a\x+ —j 



X exp < ibn 



P 

AxV 



a X 




(274) 



With the help of Eq. (265) it is possible to directly calculate how different 
aberrations modify the expression for the autocorrelation function and its 
Fourier transform. These computations can be carried out by means of nu- 
merical techniques for any value of n and 6„, in analogy with what has been 
done in [18]. 

For completeness we will now calculate the Optical Transfer Function in sev- 
eral situations, which may also be found in [18]. In addition, we will also 
present the line spread function typical of these situation, which cannot be 
easily found in textbooks. Other kind of aberrations can be treated in the 
same fashion, and pupils with different shape may be selected. We are inter- 
ested in the case when the influence of the phase error is comparable with the 
influence of diffraction effects on the pupil, i.e. when ~ 1. Here we will 
consider several cases for defocusing (n — 2), coma (n = 3), spherical aber- 
rations (n = 4) and a combination of defocusing and spherical aberrations as 
well. This last situation is per se interesting, because it illustrates how it is 
possible to improve the quality of a lens with spherical aberration by further 
introducing a defocusing aberration (in the case under study, i.e. for ~ 1). 
In Fig. 20, Fig. 21 and Fig. 22 we plot the hne spread functions describing 
these aberration cases. As an aside it is worth to anticipate here that Fig. 20 
is strictly related with the resolution of a pinhole camera setup. This will be 
demonstrated in the next Section 10. As has been already said, the knowledge 
of the line spread function is completely equivalent to the knowledge of the 
transfer function (OTF) of the system. In Fig. 23, Fig. 24, and Fig. 25 we plot 
the transfer functions relative to the same cases treated in Fig. 20, Fig. 21 and 
Fig. 22. 

In closing, it is interesting to deal with the limit when if/TV <^ <^ Dzf. 
In this case the coherence length at the pupil is much smaller than the char- 



acteristic pupil size, or, equivalently, Pa{fi)P*{f2) — l-Pa(^i) 



as 



presented in Eq. (7.2-15b) of [8]. In order to retain this last simplification when 
aberrations are present, we must require that the phase of the pupil function 
Pa (that is now a complex object) is not appreciably different when fi and 
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r2 are separated by a distance of order of the coherence length or smaller. 
This is equivalent to the requirement that the characteristic scale of the lens 
imperfections is much larger than the coherence length. Mathematically, this 
means that < 1. Under this assumption, aberrations cannot affect the 
cross-spectral density in the limit zj/N <^ <^ Dzf. This fact is known from 
a long time and, as reported in [8], it was first discovered by Zernike [19]. 
The same limiting case can be presented in the line spread function formal- 
ism. One should recall that the resolution due to diffraction effects is of order 
5i/(a\/iV), as has already been seen in the previous Section 8. This means 
that in the limit zj/N <^ <^ Dzj it makes sense to account for diffraction 
effects from the pupil, because in this case the resolution due to diffraction 
effects is worse than that related with the quasi-homogeneous approximation 

However, one may choose to worsen the resolution 

of the calculations from to the resolution due to diffraction effects, 

zi/{a\[^) <^ 1. This is equivalent to neglect diffraction effects. In this case, 
since < 1, the autocorrelation function of the pupil can be substituted 
with unity or, equivalently, the line spread function / plays the role of a Dirac 
5-function in the calculation of the intensity. Therefore, aberrations cannot af- 
fect the intensity distribution at the image plane. Moreover, the expression for 
the spectral degree of coherence Eq. (231), remains vahd for |6„| < 1 because, 
as already discussed, Pa{fi)P*{f2) — |-Pa(^i)|^ — |-Pa(^2)P- One concludes that 

in this limit, and with resolution zi/{a\[^) <^ 1, aberrations cannot affect the 
coherence properties on the image plane. 

9.2 Severe aberrations 

It is now interesting to discuss an analytical treatment valid in the case for 
^ 1, which exploits the simplifications arising from the large parameter 
Under this constraint, aberrations will be considered severe. 

9.2.1 Physical Optics prediction of the line spread function 

Under the approximation \hn\ ^ 1, it is possible to present an analytical cal- 
culation for the line spread function / which characterizes the imaging system 
in the case of quasi-homogeneous sources. Then, once the line spread function 
is known, one obtains the intensity at the image plane by convolving the line 
spread function and the intensity from an ideal system (i.e. without accounting 
for the pupil infiuence). 

Let us now focus on that term in the phase factor of Eq. (274) which is linear 
in Ax, i.e. on nft^Axx""^. We will assume that the integrand contributes to 
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the integral for all values of x inside the interval [—1,1], otherwise the au- 
tocorrelation function would be suppressed, as the effective integration range 
would be smaller than [—1,1]. Then, a typical scale of the autocorrelation 
function is obtained in terms of Ax by imposing Axx""-^ ~ 1. In fact, 
as n\hn\^xx^~^ > 1 the integrand starts to exhibit fast oscillatory behavior, 
thus suppressing the integral. Thus, the characteristic scale ^Xtyp ~ l/(n|6„|) 
is found. Since we assumed ^ 1, we can state that, with accuracy l/|6n|, 
the functions P inside the integrand can be substituted with unity and the 
nonlinear phase factors in Ax'^ with k — 2,3... can be neglected, at least 
for reasonable orders of n, as they would give rise to typical scales of order 
'^l\bn\^^^ ^ 1/l^'nl- As a result we obtain the following major simplification: 



1 

\{al^x) = j 



dxexp 



inbnAxx 



n-l 



(275) 



Eq. (275) can be integrated analytically for all values of n (therefore including 
defocusing, coma, spherical or higher order aberrations). After definition of 



T{n,bn,Ax) 



(n — 1) [—inbnx] "-i 
n 

X 



(n- i)r 



-r 



n — IJ \n — 1 
r(s, z) being the incomplete Euler gamma function 



—inbn^x 



, (276) 



oo 

r(s,^) ^ j dt f~^ex.p[-t] , 



(277) 



we have the following result: 

Pa(dAx) = Re [T (n, bn, Ax)] + U{n) ■ Im [T (n, bn, Ax)] , 



(278) 



where n(n) is the parity of n, i.e. n(n) = if n is even and n(n) = 1 if n is 
odd. Eq. (278) is valid for any value of n > 1. 

It is interesting to compare the shape of the autocorrelation function obtained 
in the limiting case 3> 1 with that obtained with numerical calculations 
which do not exploit the simplification based on the large value of the pa- 
rameter \bn\. They rely, instead, on the exact formula for the autocorrelation 
function, Eq. (274) . This gives a visual idea of the accuracy of the asymptotic. 
Fixing bn = 97r we plot the autocorrelation function for defocusing aberrations, 
with n — 2, in Fig. 26 and Fig. 27. For coma aberrations, with n = 3, we plot 
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the real part of the autocorrelation function in Fig. 28 and Fig. 29, while the 
imaginary part is plotted in Fig. 30 and Fig. 31. The function T should be 
truncated as Ax > 2. 



Explicit substitution of Eq. (275) in Eq. (271) gives 



l{ri) — j dAx j dxexp inbnAxx 

— oo —1 



n-1 



exp 



.a . 

—l — AxTi 

di 



(279) 



Finally, exchange of the integration order and calculation of the integral in 
dAx yields the following expression for the line spread function, provided that 
\bn\ > 1: 



-1 



dx S I nbnx" ^ — -Ti 
\ di 



(280) 



Eq. (280) may be exphcitly evaluated with the help of the new integration 
variable y — Care must be taken in separating the cases when n is 

even and when n is odd. 



When n is even we obtain 



nbn 



2(n-l) 



dy 5\y - -^ri 



d _\ f abs(y)' 



n-2 
' n-1 



-nbn 



2(n-l] 



-rect 



^ 2dinbn , 



di ) V "^^n 

d abs (f .j)\ 
diubr,. I 



n-2 
n-1 



(281) 



Here the function rect(a;) is defined, as before, following [8], and is equal to 
unity for ^1/2 and zero otherwise. When n is odd we have 



^ nbn / " \ 

/(f.)=;^/ci,5h/-lf. ( 



n-1 



-rect 



di 

dfi 1 



, diJibr, 



y_ 

.nbr, 

dfi 
, diubr, 



n-2 
n-1 



n-2 
'n-1 



(282) 



Introduction of the new variable r" = dri/ {diubn) allows to write Eq. (281) 
and Eq. (282) in a more compact way, 



n-2 
1 



(283) 
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when n is even, and 



1 / 1\ "-2 

l{f") = ^— rect [f" --)■ {f"y— (284) 
Til 1 V 2 / 

when n is odd. Note that we have normahzed Eq. (283) and Eq. (284) in such 
a way that integration of I in df" over the real field gives unity. 

As we will see later on, Eq. (283) and Eq. (284) can be found with the help of 
Geometrical Optics alone. We will refer to such derivation as the Geometrical 
Optics prediction of the fine spread function. We plotted Eq. (283) or Eq. (284) 
as a function of f" for different aberrations. In Fig. 32 we plotted the case of 
defocusing aberration , in Fig. 33 the case of coma and in Fig. 34 the case 
of spherical aberrations. Also, in these figures, comparison with numerical 
calculations is shown for the severe aberration cases 62 = Stt, 63 = Ott and 
64 = Qtt. Note that Eq. (283) is symmetric in f" (when n is even I is symmetric) , 
while Eq. (284) is not (when n is odd I is not symmetric). This is consistent 
with the fact that the Optical transfer function is real in the case n is even, 
while it has a non-zero imaginary part in the case n is odd (see, for instance. 
Fig. 23, Fig. 24 and Fig. 25 and, later on. Fig. 26, Fig. 28 and Fig. 30). The 
same behavior is also found in Fig. 20, Fig. 21 and Fig. 22. Furthermore it 
should be noted that the line spread functions in Fig. 33 and Fig. 34 are not 
convergent for values of f" near zero. However, the meaning of the line spread 
function is that of the imaged intensity from a line input, but a line input is 
not physical, and must be represented in terms of a generalized function, a 
5-Dirac function. It is not surprising that such an object may lead to a singular 
result. 



9.2.2 Physical Optics and Geometrical Optics 

In the previous Section 9.2.1 we presented a calculation of the line spread 
function / for the case ^ 1 based on Physical Optics. This gave us Eq. 
(280) for the line spread function or, equivalently, Eq. (281) and Eq. (282), that 
are an exphcit evaluation of Eq. (280) in the case n is even or odd. We have 
seen that the analytical treatment in Section 9.2.1 follows mathematically from 
a major simplification of Eq. (274), arising from the large parameter ^ 1. 
Physically, this condition means that effects of diffraction from the pupil can 
be neglected. Consider the expression for the pupil autocorrelation function, 
Eq. (272). On the one hand, when x ~ the phase term due to aberration 
effects in Eq. (272) becomes comparable to unity, thus leading to oscillatory 
behavior of the integrand. On the other hand, the pupil finite aperture limits 
the integration in Eq. (272) for a; ~ 1. Therefore, ^ 1 diffraction effects 
can be neglected, while they start to become important when < 1. It 
follows that it must be possible to obtain results in Section 9.2.1, which have 
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Fig. 26. The transfer function in the presence of defocusing aberration (n = 2, 
62 = Ovr). Numerical techniques have been used to calculate the exact autocorrela- 
tion function that is compared to the analytical evaluation of the autocorrelation 
function in the severe aberration asymptote. 

been derived with the help of physical Optics, with the help of Geometrical 
Optics only. In the next Section 9.2.3 we discuss how this can be done. 

Before doing that it is worth to discuss the relationship between Geometrical 
and Physical Optics, which may otherwise be misleading. We should make 
clear that when we discuss about Geometrical or Physical Optics we are talking 
about possible ways of calculating the line spread function of the system 

In Section 8 we have introduced the concept of line spread function I and 
we have demonstrated that, under the assumption that the virtual source is 
quasi-homogeneous, / constitutes a sort of passport for a given lens. It relates 
the intensity from any quasi-homogeneous source imaged with a perfect lens 
and the intensity obtained by using a particular non-ideal lens. The intensity 
from a specific optical system can be recovered as a convolution of the intensity 
obtained in the case of an ideal optical system and the line spread function. 
The intensity from an ideal system at the image plane is, by definition, the 
intensity at the virtual source. Therefore, when the quasi-homogeneous ap- 
proximation is applicable, it is always possible to break the imaging problem 
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Fig. 27. An enlarged version of Fig. 26. Numerical techniques have been used to 
calculate the exact autocorrelation function that is compared to the analytical eval- 
uation of the autocorrelation function in the severe aberration asymptote. 

into two separate problems. First, specify the intensity distribution of the 
source. Second, specify the optical system through the I function. 

When diffraction effects are neghgible with respect to aberration effects, i.e. 
when \bn\ ^ 1, the line spread function of the lens can be calculated by both 
Physical Optics considerations (as in Section 9.2.1) and Geometrical Optics 
considerations, as in the next Section 9.2.3. In this finite aperture size 

does not influence the calculation of the line spread function. It is responsible 
for the quantity of the total energy transmitted only. Once the source charac- 
teristics are specified, one may use a ray-tracing code to get the image intensity 
from a non-ideal system or, equivalently, one may calculate the / function and 
convolve with a scaled version of the intensity on the virtual source. When 
diffraction effects are not negligible anymore, / can be evaluated with the help 
of Physical Optics considerations only. In this case, use of ray-tracing codes 
to solve the imaging problem makes no sense. Yet, the quasi-homogeneous 
approximation allows one to use a line spread function approach. Convolu- 
tion of / with a scaled version of the intensity on the virtual source solves the 
imaging problem. From this viewpoint the quasi-homogeneity of the source is 
an a priori condition with respect to the possibility of applying Geometrical 
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Fig. 28. The real part of the transfer function in the presence of coma aberration 

{n = 3, 63 = Ovr). Numerical techniques have been use to calculate the exact autocor- 
relation function that is compared to the analytical evaluation of the autocorrelation 
function in the severe aberration asymptote. 

Optics for the solution of imaging problems from a non-ideal setup. 

A final word of caution should be spent regarding notations historically used to 
represent aberrations (sec, for instance, [20]). On the one hand, the knowledge 
of the phase error $(f) = 6„(f/a)" is equivalent to the knowledge of the 
surfaces of equal phase, i.e. of the wavefronts of the electromagnetic field. In 
this case one usually talks about "wave aberration". On the other hand, the 
knowledge of the derivative d$(r) / dr (or, in more dimensions, of the gradient 

— » — * 

V$(f)) is equivalent to the knowledge of the vector field orthogonal to the 
wavefronts of the electromagnetic field. In the case diffraction effects are not 
present, one may identify d^{r)/dr with the extra angular displacement of a 
ray and recover the deviation of the transverse coordinate of a ray on the image 
plane. In this case, usually, one talks about " Geometrical aberration" [20] . It 
should be clear though, that the presentations in terms of "Wave aberration" 
and " Geometrical aberration" are completely equivalent from a mathematical 



In reference [21] the term "Rays aberration" is used in place of "Geometrical 
aberration" 
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Fig. 29. An enlarged version of Fig. 28. Numerical techniques have been use to calcu- 
late the exact autocorrelation function that is compared to the analytical evaluation 
of the autocorrelation function in the severe aberration asymptote. 

viewpoint, regardless the value assumed by the parameter with respect to 
unity. Therefore, one needs to clearly distinguish between the language used 
(Geometrical or Wave aberration) and the possibility of applying Geometrical 
Optics to calculate the line spread function. When < 1 this is not possible, 
and one has to rely on Wave Optics predictions for the line spread function 
I only. When \bn\ ^ 1 one may rely both on Wave Optics predictions (see 
Section 9.2.1) or Geometrical Optics predictions (see Section 9.2.3), and the 
two predictions must coincide. 



9.2.3 Geometrical Optics prediction of the line spread function 

As has been already said in Section 9.2.2, when diffraction effects are negligible 
with respect to aberration effects, i.e. when ^ 1, the line spread function 
of the lens can be calculated by both Physical Optics considerations (as in 
Section 9.2.1) and Geometrical Optics considerations. In this finite 
aperture size does not influence the calculation of the line spread function, 
and is responsible for the quantity of the total energy transmitted only. 
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Fig. 30. The imaginary part of the transfer function in the presence of coma aber- 
ration (n = 3, 63 = Ovr). Numerical techniques have been use to calculate the exact 
autocorrelation function that is compared to the analytical evaluation of the auto- 
correlation function in the severe aberration asymptote. 

In the Geometrical Optics limit, Maxwell equations can be replaced by the 
simpler Eikonal equation, which should be solved for surfaces of equal phase. 
In our one-dimensional case of study, these surfaces are indeed lines on the 
plane f — z and may be indicated with the family 0p(f, 5), where p identifies 
a particular value of the phase. Once the functions (j)p are known, one can 
recover the usual ray-tracing techniques remembering that rays are, at any 
point, normal to surfaces with equal phase. The following ray equation holds: 

s(f, z) = V(t)p{f, z) , (285) 



where s indicates a vector field tangent to the rays expressed in normalized 
units. If the lens is ideal one recovers, at the image plane, the intensity profile 
of the source reversed and magnified. In particular, a line input would be 
mapped to a line on the image plane. However, aberrations modify the surfaces 
of equal phase in Eq. (285), because a phase error is to be added to 0p . As a 
result, when calculating .?by means of Eq. (285), one obtains an extra angular 
displacement in normalized units for each ray dependent on the transverse 
position of the ray, that is 
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Fig. 31. An enlarged version of Fig. 30. Numerical techniques have been use to calcu- 
late the exact autocorrelation function that is compared to the analytical evaluation 
of the autocorrelation function in the severe aberration asymptote. 

a(f) = ^ = 1,Jty" . (286) 
dr a \a I 



It follows that the coordinate of the ray going through the lens at transverse 
position f has a transverse position, at the image plane, given by 

fi{f)^^bJty ' . (287) 
a \a I 



The relation between the coordinate of a ray at the pupil and its transverse 
position at the image plane for n = 2 (defocusing) , n — 3 (coma), and n — A 
(spherical aberration) is plotted in Fig. 35. 



If the lens is hit by a finite number of rays, the output of the system at the 
image plane is constituted by a finite sum of Dirac-(5 functions. Each of these 



can be represented in the implicit form S 



n-1 



. Therefore, if the 



lens is homogeneously illuminated by an infinite number of rays from the input 
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Fig. 32. The line spread function versus the reduced coordinate on the image plane 

f" = ari/ {diubn) in the case of severe defocusing aberration(n = 2, 62 = 97r) and 
comparison with the geometrical optics prediction. 

line source, we obtain the line spread function 



a 

m) = / 



df5 




n—l 



(288) 



Finally, using the new integration variable x = f/a and normalizing I so that 
the integral of / gives unity yields 

l{n) = I dx5 U - ^^n^""') ' (289) 



that is equivalent to Eq. (280). From a ray-tracing viewpoint, the problem of 
calculating the line spread function reduces to the problem of transforming a 
uniform distribution of rays into a non-uniform distribution related to the non- 
linear transformation Eq. (287). When < 1 instead. Geometrical Optics 
cannot be used to calculate the / function, and the Eikonal approximation 
fails. As we have seen in Section 9.2.1 this is equivalent, in the language of 
Physical Optics, to a situation when the simplification in Eq. (275) does not 
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Fig. 33. The line spread function versus the reduced coordinate on the image plane 
f" = afi/{dinbn) in the case of severe coma aberration (n = 3, 63 = 9tt) and 
comparison with the geometrical optics prediction. 

hold. 



10 Pinhole optics 



Taking advantage of Eq. (133), we will now study the case when images of a 
virtual source can be obtained with the help of a pinhole (i.e. a pupil without 
a lens), without further lenses or mirrors. When a pinhole can be treated 
as an imaging system, people refer to it as an (X-ray) pinhole camera. Here 
we will consider the geometry in Fig. 36. The study of this relatively simple 
setup will be helpful to reach a better understanding of Section 11, dedicated 
to imaging in the focal pane. Section 12, where we will describe imaging in 
any plane behind the lens, and Section 13, that will deal with the depth of 
focus of an imaging system. Moreover, it will also suggestively show how a 
problem apparently not related with the theory of aberrations (in the pinhole 
camera setup there is not even a lens) can formally be treated like a defocusing 
aberration problem. This is due to the appearance of a quadratic phase factor 
in the equation for the intensity at the image plane. 
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■ Numerical calculation: n=4, b =97i 
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Fig. 34. The line spread function versus the reduced coordinate on the image plane 
f" = afi/ {diuhn) in the case of severe spherical aberration (n = 4, 64 = 97r) and 
comparison with the geometrical optics prediction. 

An X-ray pinhole camera has the same properties as the more familiar visible 
light pinhole-camera. The main advantage of such a lensless imaging system is 
that a pinhole is easier to fabricate than a lens. Pinhole cameras can be com- 
bined with X-ray Synchrotron Radiation sources and detectors for a number 
of relatively speciahzed applications [22] . 

The conditions under which the pinhole can be treated as an imaging system 
are non-trivial, and are not always satisfied. With the help of Eq. (133) we can 
investigate the properties of the image in the limiting case for / — > 00, i.e. 
when there is no lens. We will restrict ourselves to the one-dimensional case, 
thus simplifying the vectorial notation in Eq. (133) to scalar notation. The 
assumption of separability of the cross-spectral density in the horizontal and 
in the vertical direction (A^^^ 3> 1 and Dx ^ 1) suggests, in fact, to discuss 
horizontal and vertical directions separately. In general, one can see that the 
following conditions must be satisfied in order to form an image of the source 
(in one dimension): 

1. The pinhole must be in the far field. In this case, using Eq. (127), the 
cross-spectral density on the pupil plane is given by 
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Fig. 35. Amount by which a geometrically traced ray departs from a desired location 
in the image formed by an optical system. The ordinate for each curve is the height 
at which the ray intersects the image plane. The abscissa is the coordinate of the 
ray at the pupil plane, n = 2 represents defocusing, n = 3 represents coma, and 
n = 4 represents spherical aberration cases according to Eq. (287). 
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G r, Af) 



1 



exp [2ir • Ar/fi] ^ (0, — r, — Af) . 



(290) 



2. In the integrand of Eq. (133) two specific phase factors appear. These 
are exp[2if'- Af'/(i2 — Zi)] and exp[2if'- Af'/^i], the latter appearing through 
Eq. (290). Both must be neghgible. 

3. The pinhole size must be larger than the coherence length on the pinhole 
plane, or the pupil functions P in Eq. (133) will modify the dependence of 
G{z2) on f (i.e. the image will not have a good resolution). 

If conditions from 1. to 3. above are satisfied, the pinhole camera works as an 
imaging system forming an inverted image of the source. From Eq. (133) one 
may sec that the image is magnified of a quantity |M| = ^2/^1, because of the 
same reason as the lens. 

We will now give a physical interpretation of conditions 1. through 3. stated 
above. We begin considering the limiting case when ^ 1 and ^ 1, i.e. 
a Gaussian quasi-homogeneous virtual source. Further on we will see up to 
what extent this assumption can be relaxed. We assume a large magnification 
constant |M| ~ d/zi ^ 1, where d — Z'z — z^. Here this assumption will be 
accepted for simplicity and relaxed, later on, to an arbitrary value of |M| . 

The second condition requires that two distinct phase factors in Eq. (133) 
may be neglected. The assumption d ^ ii leads to the single requirement 
fAf/ii <^ 1. On the one hand, f is limited by the presence of the pinhole, 
i.e. we must impose f < a. On the other hand, Af is limited by the coherence 

length at the pinhole, i.e. Af < ii/y^, because otherwise the cross-spectral 
density G in Eq. (133) drops to zero. As a result we obtain that the second 
condition given above can be expressed in mathematical terms by 



It is possible to give a clear interpretation of condition (291) in terms of 
Geometrical Optics. In fact, on the one hand, the minimal geometrical spot 
size from a line source is given, at the image plane, by |M|a. On the other 



hand, the size of the image is of order |M| yiV, by definition of magnification 



|M|. Then, in order to have a good resolution, we must require that the image 
size of a point source be much smaller than the image size of the object, i.e. 



^jN/a > 1, that is condition (291). 

The third condition given above requires that the pinhole size be larger than 
the coherence length on the pinhole plane, otherwise the pupil functions P in 




(291) 
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Eq. (133) would modify the dependence of 6(^2) on r (i.e. the image would 
be influenced by the pupil). This can be mathematically stated by requiring 
that 



d > 



(292) 



Condition (292) has a natural explanation in terms of diffraction theory. In 

fact, the diffraction spot due to the presence of the pupil can be estimated as 
d/d. In order to have a good resolution, we should impose that the diffraction 

spot be much smaller than the image size of the object, i.e. |M| / (d/d) ^ 1, 
that is condition (292). 

Finally, the first condition given above requires that the cross-spectral density 
at the pinhole position be Eq. (290), i.e. the pinhole must be in the far zone 
region. This fact can be alternatively stated by requiring that the radiation 

spot size at the pupil be dominated by the angular divergence D, i.e. z\\[t) ^ 

\[^ . In the case N < D, one should require in any case that ^ 1. This can 
be mathematically expressed by the requirement 



Zl ^ max 









D 



(293) 



Note that condition (293) and the initial assumption of a quasi-homogeneous 
source are equivalent to the requirement that the pinhole be far enough for 
the van Cittert-Zernike theorem to apply. In the case of a perfectly incoher- 
ent object (e.g. thermal hght), the source radiates over an angle 27r. Then, 
the validity of the van Cittert-Zernike theorem (i.e. the requirement that the 
radiation spot size at the pupil be dominated by the angular divergence) is 
equivalent to the condition that the transverse dimension of the source be 
much smaller than the distance between the source and the pupil. This is the 
same condition for the paraxial approximation to be applicable. As a result, 
condition (293) is always considered satisfied in usual treatments describing 
pinhole setups in the presence of incoherent objects. 

The three conditions (291), (292) and (293) can be summed up in the following: 



iV > a > 



Zl 



^ max 



N 



1 1 
-J 



(294) 



Let us now discuss about the resolution of the pinhole camera. Consistency of 
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Fig. 36. Geometry for a pinhole camera setup. S indicates the virtual source plane, 
condition (294) requires that 



iV > . 



(295) 



Although condition (294) requires that condition (295) be satisfied, it does 
not pose any constraint on the relative magnitude of with respect to Zi. 
This observation suggests the presence of another characteristic scale of the 
problem, a^/zi. This scale is linked with the resolution of the system. When 
a^/zi ^ 1 or a'^/zi ^ 1 wc arc in the presence of an extra large or small pa- 
rameter and, thus, we have two asymptotic regimes. In order to systematically 
consider this issue we start writing the expression for the intensity profile for 
the pinhole camera image with the help of Eq. (133) and Eq. (290), that is 



i(z2, r) = — / df'dAf'l g{0, -f', -Af') exp 



xP(f' + Ar')P{f' - Af') \ exp 



+ f'Af' 
. \zi d 

-—tAt 
d 



(296) 



As usual, the normalization factor Z is chosen in such a way that 7(^2, 0) = 1. 
The phase in parenthesis {...} must be neghgible under condition (291), and 
the infiuence of the pupil function P must also be negligible under condition 
(292). Here we retain both these factors, because we are interested in studying 
the resolution of the pinhole camera, i.e. the accuracy of our calculations. 



The phase contribution in 2/d is always much smaller than the phase contri- 
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bution in 2/5i, since we assumed |M| 3> 1. The integrand contributes to the 
integration results only for those values of f' and Af' for which the phase factor 
is not much larger than unity, otherwise the integrand exhibits fast oscillatory 
behavior and it effectively averages to zero. As a result we may neglect, with 
some accuracy over the accuracy of the integral, the phase term in l/o?, that 
will start oscillating for higher values of f' and Af'. 

Note that, under the accepted assumptions iV 3> 1 and Z) 3> 1, the expression 
for the Fourier transform of the cross-spectral density at the virtual source 
position is given by 



^(0,f', Af')=exp 



2iVAf 
^1 



12 



exp 



2zlD 



(297) 



in agreement with Eq. (145). 



Analysis of Eq. (297) shows that the exponential function in f' exhibits a char- 
acteristic scale of order Zi\[t>. Conditions (291) and (293) require z\\[t) ^ 

^ a. As a result, since zi\[t) ^ a, we can approximate the exponential 
function in f' in Eq. (297) with unity to obtain an equation still suitable for 
investigating the resolution of the pinhole camera, that is 



I{^2,r) = ^ / dAf'exp 



-oo 
• oo 



2iVAf'2 



exp 



2^ - A 

-Ti\T 

d 



X • 



exp 



i\ — \ r'Af' 



P{f' + Ar')P{f' - Af ') 



(298) 



The quantity in parenthesis {...} in Eq. (298) is the autocorrelation function 
of the pupil. It accounts for a phase error, exactly as in the case of aberrations 
of the second order (defocus). Condition (292) states that the characteristic 
scale of Af ' in Eq. (297) is small compared to a, that is the characteristic scale 
of the pupil function P. If also a'^/zi ^ 1 we see that the characteristic scale 
of P is large compared with the scale imposed by the phase in parenthesis 
{...} in Eq. (298), because r'Ar'/zi > 1 at Af ~ a and r' '-^ a. Therefore, 
under the assumption / ii ^ 1 we may neglect the dependence of P on Af ' 
in Eq. (298) and obtain 



i{z2,r) 



oo 

i/ 



dAr'P{Ar') exp 



2iVAf 

^1 



/2 



exp 



2^ -A ^; 

-—rAr 
d 



(299) 
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where 



P(Af') = J df' \P{f')\'^ exp 

— oo 

or, equivalently : 



■-2 

^— r Ar 



(300) 



P(Af') = / cir'exp 



z— r Ar 



(301) 



Except for an unessential multiplicative factor d, Eq. (301) is formally equiv- 
alent to Eq. (275) with n=2 and 



b2^ — :$>l. (302) 

This means that we may study the problem of the resolution of the pinhole 

camera as an aberration problem: in particular, a defocusing aberration. The 
parameter range when d'^/zi ^ 1 leads to equations similar to the case of 
severe aberrations when ^ 1, treated in Section 9.2. The autocorrelation 
function of the pupil, that is the Optical Transfer Function of the system, is 
then obtained by integration of Eq. (301) and reads: 



P(Af') — sine 




(303) 



where Eq. (301) has been used and an unessential multiplicative factor 2d has 
been neglected. Substitution of Eq. (303) in Eq. (299) yields 



I{z2,r) 



oo 

i/ 



aAr smc — Ar exp 



27VAf'2 

Zl 



exp 



2^ -A- 
-— r Ar 

d 



(304) 



Eq. (304) allows an estimation of the resolution by taking the ratio of the 
width of the sine and of the exponential function in Ar'. On the one hand. 



We assume, as done before, that no apodization is present. 
^^It should be noted that we switched back from notations x and Ax, used in 
Section 9, to our usual notation f' and Af', with f' = x and 2Af' = aAx. We 
remind that the reason why we used notations x and Ax in Section 9 was for the 
reader's convenience, as these notations aUow direct comparison with aberration 
theory developed in standard textbooks. 
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the width of the exponential function is of order zi/{2yN). On the other 
hand, the width of the sine function is of order zi/{2a). As a result, when 

a^/zi » 1 the resolution of the camera is of order a/\[^. This procedure 
is justified by the fact that Eq. (304) can be interpreted as a convolution of 
a rectangular profile with a (new) Gaussian function, and that the width of 
these two functions can be obtained by taking the inverse widths of the sine 
function and the Gaussian function in Eq. (296) and by multiplying them by 
the factor d. 

Let us now deal with the case when / Zi <^ 1. Going back to Eq. (298), we see 
that a is narrow compared with the scale imposed by the phase in parenthesis 
{...} in Eq. (298), because f'Af'/ii ^ 1 at Af' ~ a and f' ~ a. Therefore 
we can neglect the phase factor in parenthesis {...}, which corresponds to a 
case with no aberrations. The autocorrelation function of the pupil can now 
be written as a triangle function 



P(Af') = tri l^^j , (305) 

that should be substituted into Eq. (299) to give the analogous of Eq. (304) 
in the limit a? / Zi <C 1, that is 



. _ 1 7 Mr'' 
7(^2, ~ '2 J "^^^'^^^ I I ^^P 
—00 



a 



27VAf'2 

^1 



exp 



■ 2i _ ^ 
— -r/\r 

■ d 



(306) 



Eq. (306) can also be interpreted as a convolution of a sinc^(-) profile with a 
Gaussian profile. 

Similarly to Eq. (304), Eq. (306) allows an estimation of the resolution by 
taking the ratio of the width of the triangular function and of the exponential 
function in Af' . As before, on the one hand the width of the exponential 

function is of order zi/ (2\/i^). On the other hand, the width of the triangular 
function is of order 2a. As a result, when d^/zi <^ 1 the resolution of the 

camera is of order zi/{d\[^). 



When d? / zi ^ 1 the resolution of the camera is of order d/yN. One has 
better resolution as the pupil becomes smaller and smaller, but the condition 
d'^/zi ^ 1 becomes less and less satisfied. When d'^/zi ^ 1 the resolution of 

the camera is of order ii/(a\/iV). In this case one has better resolution as the 
pupil becomes larger and larger, but the condition d'^/zi -C 1 becomes less and 
less satisfied. As a result there must be an optimum for pupil apertures of order 
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~ ii. This optimum depends on the object considered (in this discussion, for 
instance, we assumed a Gaussian source) and on the definition of the width of 
a function, that may vary depending on circumstances. However, starting from 
Eq. (298), we may present an expression for the Optical Transfer Function of 
the pinhole camera, which is the quantity in parenthesis {...}. Such quantity 
can be written as 



Ppc(Af', a, zi) = j df'P{f' + Af')P{f' - Af') exp 

—a 

After introduction of ^ = f'/a and Q = a^/zi we can write Eq. (307) as 



i2- 



r' Af' 



Zl 



(307) 



Note that the limit fl — > oo of Eq. (308) is a sinc(-) function, while the limit 
Q — > is a tri(-) function, as it should be. Eq. (308) enters in the expression 
for the intensity as the term in {...} in Eq. (298), that can be reinterpreted as a 
convolution product. Then, the Fourier transform of Eq. (308) with respect to 
Af' /a is one of the factors in this convolution product, and its width is related 
with the resohition of the pinhole camera. This Fourier transform is the line 
spread function for the pinhole camera. After introduction ol y = afi/d and 
fj = Af /a we can write the line spread function for the pinhole camera as 



Af 



(308) 



/pe (fi, y) = I I / di^P (e + r/) P (e - r/) exp [^2r/ {i^ - y)]^ , (309) 

where we accounted for the fact that the Fourier transform integral in d^ 
is limited by the presence of the pinhole to the range [—1,1]. On the one 
hand, the integral in parenthesis {...} in Eq. (309) is the Fourier transform of 
P{i + v)P{C ~ v) calculated with respect to as a function of 2{^Q, + y). On 
the other hand, the function P{C, + ti)P{C, — 1]) is a window function similar 
to the pupil function P. It is equal to unity for values of |?7| < 1 — |^|, and it 
is zero elsewhere. Therefore, the quantity in parenthesis {...} in Eq. (309) can 
be calculated analytically yielding 



1 



^p,(n..)^/d rP«""!><' (310) 

J E\l + y 



Note that if Q, — > we have 
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Fig. 37. The line spread function Ipc with Q = 0.5. In this case the asymptote for 
— > 0, Eq. (311), is well-matched to the numerical evaluations. 



Ipci^^y) - I d^sm[2y{l-m = 2sinc^(|/) , 

y 



(311) 



as it should be, while if Q — > oo, after the change of variable ^ — > ^' — flC, 
we can rewrite Eq. (310) as 



-a 



e+y 



{Si[2{y + n)] - Si[2{y - m , 



(312) 



where Si(-) indicates the sin integral function. It can be seen that, as — > oo, 
the function defined by Eq. (312) approximates more a more a rectangular 
function which is constant for —ft < y < Q and equal to zero elsewhere, as it 
should be. 



We may use a new variable y' = y/V^ in Eq. (310), thus obtaining: 
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Fig. 38. The line spread function /pc with O = 2. In this case the asymptote for 
Q — > 0, Eq. (311), starts to diverge from numerical evaluations. The width of Ipc 
is close to its minimum. 



sm 



2{^n + y'VQ){l-\^\) 



/pe {n, y') = / di ^ ^ . (313) 



The reason for this is that, now, both the asymptotes Eq. (311) for O ^ 1 and 
Eq. (312) for ^ 1 present a characteristic width 1 / and respectively. 
In Fig. 37, Fig. 38, Fig. 39 and Fig. 40 we present various shapes of the /pc 
function and its asymptotic limit for different values of Jl. Once a definition 
of the width of the function /pc is chosen, the optimal operation point for the 
pinhole camera may be set requiring that the width of /pc be minimal. It should 
be remarked that the definition of the width of /pc is somewhat subjective. 

One may, for instance, define the width of /pc to be the Pull Width Half 
Maximum (FWHM). In Fig. 41 we present different plots for /pc normalized 
to /pc(0) for different values of f2 = 1,2,3,4,5,6. The minimal value of the 
FWHM of the line spread function for the pinhole camera happens to be 
located somewhere between = 4 and f2 = 5. It may be estimated to be 
about 1.5. As a result one concludes that the optimal aperture for the pinhole 



130 




Fig. 39. The line spread function Ipc with = 10. An intermediate situation between 
the two hmits for Q — ^ and Q, — ^ oo. 

is given by 



~ 4.5ii 



(314) 



or, in dimensional units, 



a ~ 4.5 — . 



(315) 



We can estimate the resolution 6r of the pinhole camera in dimensional units 
imposing 6y' to be equal to the FWHM of the line spread function, i.e. 1.5, 
and taking advantage of the definition of y'. Therefore we obtain 



Sr ~ 1.5- 



-d 



(316) 



Results in this Section have been obtained under the assumption |M| 3> 1. A 
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Fig. 40. The line spread function Ipc with 0, = 100. The asymptote in Eq. (312) 
starts to match the numerical calculations. In the limit O — > oo they both merge 
to a rectangular function. 



generalization can be presented introducing 



Zid 

zi + d 



(317) 



and the analogous dimensional value Zee = ZesLw to be used in place of zi 
and zi in Eq. (314), Eq. (315) and Eq. (316). In dimensional units we have 



a ~ 4.5 



(318) 



and 



Sr ~ 1.5- 



UJZeS 



(319) 



Before proceeding we should discuss about the applicability of Eq. (313). Up 
to now, in fact, we discussed about a Gaussian quasi-homogeneous undulator 
source (TV » 1 and D ^ 1). Eq. (313) can be applied in a wider variety of 
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Fig. 41. Relative line spread function Ipc/lpd^) for different values of il. 

cases. In particular, Eq. (313) can be applied when (a) the virtual source is 
quasi-homogeneous, (b) the pinhole camera is installed in the far zone and 
(c) <S max[l,£)]i^. When conditions (a), (b) and (c) are satisfied, the 
line spread function for the pinhole camera, /pc can be calculated from Eq. 
(313). At the beginning of the present Section 10 wc stated three conditions in 
order to have an image of the source, besides the Gaussian quasi-homogeneous 
assumption (that we relaxed, here, to the more general case (a)). Condition 
1. given at the beginning of the present Section 10 can be identified with (b). 
Condition 1. and 2., together, are responsible for (c). The other requirement 
was formulated in order to obtain a good quality image of the source. However 
/pc can be calculated regardless of it and, in different situations, one may have 
a good quality or a bad quality of the image, depending on how the width 
of /pc scales with the width of the ideal image. The intensity profile at the 
observation plane located at z = Z2 is given, in any case, by the convolution 
product between the line spread function /pc and the ideal image. In particular, 
for an undulator source with Gaussian intensity profile at the virtual plane 
such a product is given by /pc * exp[— f^/(2M^iV)], where |M| = d/zi. When 

3> a 3> Zi/ \[^ a good quality image is formed at the observation plane. 
In other situations instead, the width of /pc is of order of, or even larger than. 
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the width of the ideal image, and a bad quality image will be formed. To 
sum up, in order to calculate the intensity distribution at the observation 
plane located at position z = 2.2 behind the pinhole, we should first calculate 
Q — a} jZfS^. Then use Eq. (313) to calculate the line spread function l^^. 
Finally, convolve with the ideal image from the source. The minimal width 
of the line spread function, which correspond to the best accuracy for the 
image, is for a^/5cff — 4.5. Note that conditions (a), (b) and (c) above do not 
require the source to be an undulator source. In particular, (a), (b) and (c) are 
automatically satisfied for thermal sources, or perfectly incoherent objects. As 
a result, our theory can be applied for visible light imaging as well. 

These remarks simplify our discussion and drastically widen the applicability 
region of our results. In particular they can describe bending magnet sources. 
For instance, at ESRF (European Synchrotron Radiation Facility) a pinhole 
camera setup has been developed for electron beam diagnostics [22]. Syn- 
chrotron Radiation from a bending magnet is imaged through an X-ray pin- 
hole camera setup. Subsequent analysis of the pinhole camera image allows 
one to retrieve the electron beam sizes. In the situation studied in [22] the 
diffraction angle is of order of 7"^ ~ 10~^ rad. The system is reported to 
generate an image of the source when the pinhole is moved at about Z\ — 23 
m from the dipole magnet. At this distance, and at the critical wavelength, 
the spot size due to divergence of the single particle radiation is estimated 
to be about 7^^,21 ~ 2mm ^ (^x.y^ cFx,y being the transverse electron beam 
sizes. This guarantees that the pinhole is in the far field, i.e. condition (b) is 
satisfied. The critical wavelength is of order A ~ 3 • 10~^^ m, since one may 
read, in [22], that "the contribution from difl'raction assumes a typical photon 
energy of 40 keV". The radiation diffraction size can be estimated to be of 
order A7/(27r) ~ lO^^/um ^ (yx.y- which demonstrates the quasi-homogeneity 
of the source, i.e. condition (a). In reference [22] one can also find c/j = 11.5 
m. We obtain an optimal pinhole aperture of a ~ 12//m. Since the pinhole 
aperture is much smaller than the spot size due to divergence of the single 
particle radiation, condition (c) is satisfied as well. As a result our theory can 
be applied to the situation treated in [22]. Our results should be compared 
with the actual choice in [22]. From Fig. 2 in that reference we can conclude 
that and 2a ~ 50/xm in the horizontal direction. This choice is not optimal, 
since corresponds to the value Vl ~ 10. The line spread function for the pin- 
hole camera in the horizontal direction is illustrated, in this case, in Fig. 39. 
Estimation of the resolution is not easy because, as already said, is related to 
the definition of the width of /pc, which is quite subjective in the particular 
case depicted in Fig. 39. However, knowing the profile of the fine spread func- 
tion, we may use it in order to deconvolve experimental results and extract 
the electron beam size. From the same Fig. 2 in [22] we can also conclude that 
2a ~ 25/xm in the vertical direction. Although the definition of the width of 
/pc is somewhat subjective, we can conclude that the choice made in [22] is 
near the optimal value in the vertical direction. The resolution in the vertical 
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direction turns out to be, from Eq. (319), 5r ~ 14/xm. As a result the res- 
olution in [22], that is about 26/im, is somewhat underestimated. It should 
be clear from the previous discussion that decreasing the pinhole dimension 
beyond the optimal size not only will decrease the photon flux, but will also 
worsen the camera resolution. 



11 Imaging in the focal plane 



In this Section we will investigate the intensity distribution on the focal plane 
due to a quasi-homogeneous source. We will first consider the case of a Gaus- 
sian quasi- homogeneous source {N ^ 1, £) ^ 1) and subsequently generalize 
our conclusions, as done before for the more comprehensive case of generic 
quasi-homogeneous sources. The characteristics of the intensity distribution 
in the focal plane can be treated in formal analogy with the pinhole camera 
setup treated in Section 10 and with the physics of aberrations, in particular 
defocusing aberrations, described in Section 9. This may seem counterintu- 
itive, since, at first glance, we are treating completely different systems from a 
physical viewpoint. The formal analogy between these situations is a demon- 
stration of the power of the combined Statistical and Fourier Optics approach, 
which allows one to unify study cases otherwise completely distinct. Such uni- 
fication can be seen from the expression for the cross-spectral density at any 
distance from the source in the presence of a pupil, Eq. (133). As usual, we will 
assume that the lens is in the far field, i.e. Dzf » N. Under this assumption 
we may use Eq. (127) to characterize the cross-spectral density at the lens 
position. In the image plane, the far field assumption allows cancellation of 
all phase factors in the integrand of Eq. (133) and leads to Eq. (226). This 
cancellation does not hold anymore for the focal plane. Use of the focal plane 
condition f — Z2 — zi and of Eq. (127) yields the intensity (Af / = 0) 



7(^2, rf) = ^J df'dAf'lg(0, -f', -Af') exp 



i — r'Af' 





2i 


jexp 


— -TfAr 
[ f \ 



(320) 



As one can see, the phase factor in the expression for the cross-spectral den- 
sity, Eq. (127), survived in Eq. (320) due to the choice of the focal plane as 
the observation plane. Prom a formal viewpoint Eq. (320) is identical to Eq. 
(296). In fact, as the reader will remember, the phase f'Af'/di in Eq. (296) is 
negligible for \M\ ^ 1 and can be retained without change in the formalism 
in the case \M\ is not much larger than unity by defining z^s according to 
Eq. (317). As before, the normalization factor Z is chosen in such a way that 

/>2,0) = 1. 
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Under the (for now) accepted assumptions N ^ 1 and Z) 3> 1, the expression 
for the Fourier transform of the cross-spectral density at the virtual source 
position is given by Eq.(297). Substitution of Eq. (297) in Eq. (320) yields 
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(321) 



Note that when the pupil influence is negligible (that is the case when ^ 
Dzl > zl/N), Eq. (321) reads 



H^f^ff)^^ I dAr'ex-p 
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(322) 



The reader may check that evaluation of Eq. (322) gives back Eq. (190), as 
it should be. For large non-limiting apertures, the extra phase imposes an 
extra Fourier transformation of the integrand, which gives the usual result. 
The intensity in the focal plane is a scaled version of the Fourier transform 
of the spectral degree of coherence on the virtual source plane. Yet, there are 
situations when one may recover an image of the virtual source at the focal 
plane. This happens when parameters arc such that the only influence of a 
flnal pupil aperture is to make the phase factor in parenthesis {..} in Eq. 
(321) negligible. Looking for a region in parameter space where this situation 
is realized is equivalent to what has been done in Section 10. Three conditions 
were given such that the phase factor in f'Ar' in Eq. (296) could be neglected. 
In that case, for a pinhole camera, we had image formation at all positions 
after the pinhole. 

In particular, when Dzf ^ N ^ ^ zj/N not only we can neglect the 
phase factor in Eq. (321) but we can also see that the width of the exponential 
function in Af' is much narrower than that of the pupil function P. Therefore, 
the dependence on Af' in P can be neglected. Moreover, the width of P in f' 
is much narrower than the width of the exponential function in f ', so that the 
latter can be neglected as well. As a result we simplify Eq. (321) to obtain 
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dr'dAf' exp 



^1 



P(f')rexp 



2i _ ^ 

— ^Tf Ar 
/ 



(323) 



The integral in dr' yields an unessential multiplication constant to be included 
in Z, and one is left with 



Hzf,rf) = - I dAf'exp 



2iVAf'2 

^1 



exp 



2? - A ^/ 

■— rf Ar 
/ 



(324) 



Including another unessential constant in the normalization factor Z, Eq. 
(324) can be written as 



As a result, when Dzf ':$> N ^ a,'^ ^ zf/N, wc obtain, in the focal plane, a 
scaled image of the virtual source. This is exactly what a pinhole camera does 
in the parameter region when conditions 1. to 3. in Section 10 are satisfied. 
A correspondence between completely different problems has thus been estab- 
lished thanks to the combined power of Statistical Optics and Fourier Optics 
approach. 

Such a correspondence can be pursued further, up to a complete formal iden- 
tification between the pinhole camera and the focal imaging system. To this 
purpose we restrict our analysis to the case <C Dzf. In the opposite limit 
we would obtain the already treated result for negligible pupil influence. In 
this situation we can neglect the exponential function in f'^ in Eq. (321). Eq. 
(321) can thus be simplified to 



zfrj 
2pN 



(325) 
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TfAf' 



P{f' + Ar')P{f' - Ar') 



(326) 



Once the substitutions / — > d, Zf — > Z2 and f/ — > f are made, Eq. (326) 
is identical to Eq. (298). This means that, in the limit a? <^ Dzl studying the 
pinhole camera setup is completely equivalent to studying the focal imaging 
detup. As in Eq. (298), the quantity in parenthesis {...} constitutes an Optical 
Transfer Function, and its Fourier transform yields a line spread function for 
the system. Note that in the intermediate region for ~ Dzl retain 
the same formalism: in this case though, the exponential function in f'^ in Eq. 
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(321) cannot be neglected and results would be different, depending also on 
the source parameter Dy which introduces a factor formally identical to lens 
apodization. 

However, in the case <S Dz\, we can proceed in perfect parallelism with the 
study of the pinhole camera setup in the previous Section 10. Starting from 
Eq. (326), we may present an expression for the Optical Transfer Function of 
the pinhole camera, which is the quantity in parenthesis {...}. This quantity 
can be written as 



a 

Pip(Af', a, zi) = J df'P{f' + Af')P(f' - Af') exp 



i2- 



r'Af' 



Zl 



(327) 



After introduction of ^ = r'/a and Q = a^/zi we can write Eq. (327) as 





Af'' 



. (328) 



Similarly as before, the Fourier transform of Eq. (328) with respect to Af' /a 
is the line spread function for the system. Since it refers to the focal plane it 
will be indicated with Zfp. After introduction oi y — afi/ f and r] = Af' /a we 
can write such Fourier transform exactly as Eq. (309), that is 



Zfp (Q, y) = I I / dr^P {^ + v)P(^- v) exp [i2rj (^Q - y)] | , (329) 

where we accounted for the fact that the Fourier transform integral in is 
hmited by the presence of the lens to the range [—1, 1]. On the one hand, 
the integral in parenthesis {...} is the Fourier transform of P{$^ + r))P{^ — rf) 
calculated with respect to as a function of 2{^VL + y). On the other hand, the 
function P{i + r])P{^ — r]) is a window function similar to the pupil function P. 
It is equal to unity for values of I77I < 1 — 1^| and is zero elsewhere. Therefore, as 
in the previous Section 10, the quantity in parenthesis {...} can be calculated 
analytically yielding back Eq. (310). Moreover, as in the previous Section 10 we 
may use a new variable y' = y/v'Ti in Eq. (329), thus obtaining the following 
final expression for the line spread function: 



} sin 2(cn + (/v^)(l - lel) 
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Eq. (330) is identical to Eq. (313). Therefore, the shapes of the function /fp 
and of its asymptotic hmit for different values of Q are the same as those for 
Ipc given in Fig. 37, Fig. 38, Fig. 39 and Fig. 40. As before, once a definition for 
the width of Zfp is chosen, the optimal operation point for the pinhole camera 
may be set requiring that such width be minimal. Defining the width of Zfp 
to be the full width half maximum (FWHM) one concludes that the optimal 
lens aperture for the imaging in the focal plane is given by 

~ 4.5ii (331) 

or, in dimensional units 

~ 4.5^ . (332) 



We can estimate the best resolution in the focal plane 5rf in dimensional units 
requiring that Sy' be equal to the minimal FWHM of the hne spread function, 
i.e. 1.5, and taking advantage of the definition of y', which gives 

Sr, . . (333) 



Note that results in this Section have not been obtained under the assumption 
|M| » 1 as those in the last Section. Therefore, substitution of Zi — > Zes is 
not required in this case. 

Similarly as before, we should now discuss the applicability of Eq. (330). Up 
to now we discussed about a Gaussian quasi-homogeneous undulator source, 
but Eq. (330) can be applied in a wider variety of cases. In particular, Eq. 
(330) can be applied when (a) the source is quasi-homogeneous, (b) the lens is 
installed in the far zone and (c) -C max[l, D]zf. When conditions (a), (b) 
and (c) are satisfied, the line spread function for the lens, kp can be calculated 
from Eq. (330). Depending on the situation, one may have a good quality or a 
bad quality of the image. This is related with how the width of Ifp scales with 
the width of the ideal image. The intensity profile at the focal plane (located 
at i = ^1 -|- /) is given, in any case, by the convolution product of the line 
spread function /fp and the ideal image. In particular, for an undulator source 
with Gaussian intensity profile at the virtual plane, such product is given by 

Ifp * exp[-f7(2M2iV)], where |M| = f/zi. When //^ > a > zi/\/^ a good 
quality image is formed at the observation plane. In other situations instead, 
the width of /fp is of order of the width of the ideal image, and a bad quality 
image will be formed. To be specific, for values of ^ 1 the line spread 
function and, therefore, the intensity distribution at the focal plane tends to a 
stepped profile, while for values of Q <^ 1 one obtains a sinc^(-) profile. To sum 
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up, in order to calculate the intensity distribution at the focal plane behind 
the lens, we should first calculate Q = / Zi. Then use Eq. (330) to calculate 
the line spread function /fp. Finally, convolve with the ideal image from the 
source. The minimal width of the line spread function, which corresponds 
to the best image resolution, is when a? / z\ ~ 4.5. Note that conditions (a), 
(b) and (c) above do not require the source to be an undulator source. In 
particular, (a), (b) and (c) arc automatically satisfied for thermal sources, or 
perfectly incoherent objects. As a result, our theory can be applied for visible 
light imaging as well. These remarks simplify our discussion and drastically 
widen the applicability region of our results. 



12 A unified theory of incoherent imaging by a single lens 



In the present Section 12 we will develop a unified theory which is applicable 
to imaging in an arbitrary plane behind the lens. In the previous Section 11 
we discussed the possibility of imaging the source at the focal plane. There 
we recognized that under conditions (a) the source is quasi-homogeneous (b) 
the lens is installed in the far zone and (c) a? <^ max[l, D]zl, the intensity on 
the focal plane can be calculated as a convolution product of the line spread 
function of the system and the ideal image, which is a scaled version of the 
intensity distribution on the source plane. Moreover, it was shown that the 
characteristics of the intensity distribution in the focal plane can be treated 
in formal analogy with the pinhole camera setup in Section 10 and with the 
physics of aberrations, in particular defocusing aberrations, described in Sec- 
tion 9. We have seen that when the lens is in the far field (condition (b)) 
Eq. (127) can be used to characterize the cross-spectral density at the lens 
position. In the image plane, this assumption alone allows cancellation of all 
phase factors in the integrand of Eq. (133) and yields Eq. (226). In the previ- 
ous Section 11 we showed that this cancellation does not hold anymore for the 
focal plane. From this viewpoint the focal plane is not privileged in any way 
with respect to other observation planes. It is this last remark which suggests 
a generalization of the previous results. For a generic observation plane posi- 
tioned at z — Z2 one may use Eq. (127) to get an expression for the intensity 
(Af = 0) from Eq. (133), in analogy with Eq. (320) and Eq. (296), that is 



(if'dAf'|^?(0, -f', -Af') exp 
xP(f' + Ar')P{f' - Af') \ exp 



2i{~ + ^ + i]f'Af' 

V / dj 

2^ , 
— -T'/Xr 
. d 



(334) 



where d = Z2—Z\. As one can see, the phase factor in r'Af is more complicated 
than those Eq. (320) and Eq. (296) due to the complete arbitrariness of the 
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observation plane. Nevertheless, Eq. (334) can be recast to the same form of 
Eq. (320) simply with the help of with a new definition of Zes, that is 



ZeS f Zi d 



(335) 



Accounting for Eq. (335) we may rewrite Eq. (334) as 



i(z2, f)^J df'dAf'!^g{0, -f', -Af') exp 
xP(f' + Af')P{f' - Af')| exp 

which is equivalent to Eq. (320). 



i- — r'Af' 

■ ZeS 

2^ -A ^/ 

— -rAr 
d 



(336) 



All that is left to do is now follow, step by step, the previous Section 11. This 
leads to the following results. When conditions (a), (b) and (c) are satisfied one 
may introduce, as before, a line spread function of the system, characteristic 
of the observation plane z — Z2. This fine spread function is formally identical 
to Eq. (330): 



k in, y') = j di i , (337) 



with Q, — a^/ieflf and y' = yj ^JVl, y' being defined by y = df/d. Here similarity 
techniques have been employed. The five dimensional parameters A, a, /, Zi 
and Z2 have been reduced to the only parameters Vt and y' , so that one is left 
with the calculation of a dimensionless function in y' depending on the single 
parameter f2, Eq. (337). As before, depending on the situation, one may have 
a good quality or a bad quality of the image. This is related with how the 
width of scales with the width of the ideal image. The intensity profile at 
the observation plane located sX z = Z2 is given by the convolution product of 
the line spread function and the ideal image. In particular, for an undulator 
source with Gaussian intensity profile at the virtual plane, such product is 

given by U * exp[-f2/(2M2iV)], where |M| = d/zi. When > a > Zi/\[fj 
a good quality image is formed at the observation plane. In other situations 
instead, the width of is of order of the width of the ideal image, and a 
bad quality image will be formed. To sum up, we presented here an algorithm 
to calculate the intensity distribution of radiation at any observation plane 
located at position z = Z2 behind the lens, given an arbitrary value of d 
and a system satisfying conditions (a), (b) and (c). First one should calculate 
Q = a^/zes- Then use Eq. (337) to calculate the line spread function l^. Finally, 
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convolve with the ideal image from the source. As before, when $7 ^ 1 the line 
spread function tends to a stepped profile and, as a result, the intensity profile 
reproduces a stepped profile too. Note that, being arbitrary, we cannot 
give a relation in terms of a and Zi corresponding to an optimal line spread 
function. For instance, when Z2 — > Zi, i.e. when we consider the asymptotic 
limit for the image plane, we obtain from Eq. (335) that ZeS — ^ oo and 
fl — > 0. Then, from Eq. (337) one obtains 1^ — ^ smc'^{afi/d), exactly as in 
Eq. (311). In this case no defocusing aberration is present, quadratic phase 
term having being cancelled by the particular choice of the observation plane. 
As a result, we cannot give a criterium for an optimal line spread function: we 
only have diffraction effects so that the larger the aperture a (always within 
the constraint imposed by condition (c), i.e. ^ max[l, D]zf) the better the 
quality of the image. In closing, it is worth to mention that the asymptote 
for / — > oo corresponds to the pinhole camera setup already discussed in 
Section 10. As a result, this particular case can be treated in terms of our 
unified theory as well. 



13 Depth of focus 



According to [23] "the depth of focus of a lens is the permitted displace- 
ment, away from the focal or image plane, for which the intensity on axis is 
diminished by some permissible amount" . In particular, when plane wave illu- 
mination is considered on a perfect circular lens, the focal plane corresponds 
to the plane where the radiation assumes the minimal spot size, and the in- 
tensity on axis reaches a maximum at that point. It can be shown that, in 
this case, the on-axis intensity decreases of about 20% when "the observation 
plane is displaced from the ideal focal plane [...] by an amount" (see [23]) 



lA'zl = , (338) 



where the quantity NA indicates the numerical aperture of the lens. NA = 
sin ^, 9 being the "half angle measured from the optic axis at the focus back to 
the lens" [23]. The concept of depth of focus described in Eq. (338), describes 
a case of coherent illumination of the lens by a plane wave, and the lens is 
treated, here, as a condenser. The object to be imaged is, in fact, the radiation 
source itself. Since we are dealing with coherent illumination we may say that, 
in this case, the depth of focus is diffraction limited. 

The depth of focus as described in Eq. (338) is parametrically related to an- 
other concept of depth of focus which is used, for example, in Optical Lithog- 
raphy [23, 24]. In this case one needs to illuminate a wafer with a demagnified 
image of a given pattern on a mask. Here the lens (actually the collection of 
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lenses) is no more treated as a condenser: its function is to produce the demag- 
nified image of the mask. The mask itself must be illuminated by means of a 
condenser system instead, and, as remarked in [23], "the ability to print fine, 
high contrast features is significantly affected by the degree of coherence within 
the optical system. If there exists a high degree of spatial coherence, diffrac- 
tion from adjacent mask features will interfere in the image plane, significantly 
modifying the recorded pattern". The mask should therefore be illuminated 
by incoherent light and should be considered as a quasi-homogeneous source 
itself. For a lithography setup, the resolution R is the width of the diffraction- 
limited point spread function of the system, while the distance X over which 
the image is in proper focus is, more quantitatively, the distance over which 
R is increased by some permitted amount. Although a lithography setup con- 
stitutes a completely different setup with respect to the condenser system 
illuminated with coherent plane waves, R and X turn out to be respectively 
[24]: 

R = h^; X = k,:^. (339) 



The distance X is therefore parametrically related to the distance \A'z\ defined 
in Eq. (338) and even though these two quantities refer to very different setups, 
involving quite different physics, also X is named, as \A'z\, depth of focus. 

To complicate the situation further one is frequently interested in the depth 
of focus I A^l for a condenser system when the lens is not illuminated by plane 
waves but by other kind of non quasi-homogeneous or quasi-homogeneous 
sources, which is a different situation from both cases considered in the pre- 
vious discussion. In this Section we will treat the case of a condenser lens 
illuminated by a quasi-homogeneous source characterized by AT^ 3> 1 and 
Dx ^ 1. Separate treatments of the x and y direction are thus allowed, which 
simplify to one- dimensional cases. Our study applies to the horizontal direc- 
tion only. However, if also A^^ ^ 1 and Dy ^ 1 the same results for the 
horizontal direction can be applied to the vertical direction as well. Our defi- 
nition of depth of focus will be relative to the plane of smallest spot size of the 
radiation. Therefore, in the present Section 13, the depth of focus of a lens is 
defined as the permitted displacement, away from the waist plane, for which 
the intensity on axis is diminished by some permissible amount. 

Before starting to discuss about the depth of focus, it is necessary to derive an 
expression for the point where the radiation spot size is the smallest. This can 

be obtained with the help of Eq. (216) in Section 7. In fact, since 2\J e^{z2) is 
the characteristic width of the Gaussian radiation spot size at the observation 
plane located at i = ^2, it is sufficient to look at the point -22bcst where the 
derivative of 6/3(^2) is zero to obtain the position of the waist plane. With the 
help of Eq. (216) one finds 
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^ _ (/ + zi)n- bfzj + ml _ (/ + ii) Azj - (/ - h) bzj 

^2best — ^^/^ n2 ^ ^/-s x2 

7V + L'(/-Zi) Azl + D[f-z^) 

(340) 

It is interesting to study two limiting cases of Eq. (340) for A ^ 1) and for 
A — > respectively. Consider first the case when A^ D. From Eq. (340) one 
may see that, in this first case, i2best — ^ % = Z+^i- This means that the waist 
plane asymptotically goes to the focal plane of the lens. Consider now the case 
when A — > 0. Again from Eq. (340) one may see that, in this second case, 
52best — ^ Zi = zf/{zi — /). This means that the waist plane asymptotically 
goes to the image plane of the virtual source (located at i = 0). 

Prom the analysis of the two limiting cases ior A ^ D and for A — )■ 0, we 
conclude that one always has Zi + f < i2best < ^i- In other words, the waist 
plane is always located between the focal and the image plane. This results 
may seem counterintuitive. In fact, for our Gaussian virtual source, the waist 
is located at Zg — 0, which is imaged, reversed and magnified, at the image 
plane. Therefore, at first glance, the smallest spot size should be located at 
the image plane. The reason why this is not the case is due to the presence of 
the magnification factor, which linearly increases with Z2. On the one hand the 
source has a waist located at is = and its size increases symmetrically as one 
moves away from z^ = in both the positive and the negative direction. On 
the other hand the magnification factor increases linearly with Z2 and is not 
characterized by the same symmetric dependence on the displacement from 

In the following Section 13.1 and Section 13.2 we study the problem of the 
depth of focus in a condenser system with a Gaussian quasi-homogeneous 
source. Effects from the pupil width will be neglected in Section 13.1, Finite 
pupil dimensions will be accounted for in Section 13.2. 

13.1 Large non-limiting aperture 

We will discuss, for simplicity, the case when |M| = di/zi <^ 1. In this case, 
from Eq. (30) we obtain di ~ /. Since we assumed as a starting point |M| <^ 1 
it follows that / <^ £i, and that the distance between the image and the focal 
plane is A = — / ~ j/zi — |M|/. Moreover, the radiation spot size at 

the image plane is known to be of order |M|\/i^ from Eq. (157), while the 
radiation spot size at the focal plane is obtained from Eq. (154) and is of 

order \[^f. As has been discussed above, the position of the waist goes from 
the focal plane to the image plane as we pass from the near to the far zone. 
This can also be seen by comparing the radiation spot sizes at the image and 
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at the focal plane in the near and in the far zone. In the near zone D. It 

follows that V^/ ^ y/^f/^ij the radiation spot size at the focal plane is 
much smaller than that at the image plane. On the contrary in the far zone 

A<^ D. It follows that ^/^f » \[^f/zi ~ \[^di/zi, i.e. the radiation spot 
size at the focal plane is much smaller than that at the image plane. We will 
analyze these two cases separately. 



13.1.1 Far zone 

In this case A <^ D and the waist is near the image. On the one hand, as has 
been already remarked, the radiation spot size on the image plane is about 

|M|'\/7^ and does not depend on D. On the other hand, the radiation spot size 
on the non-limiting pupil aperture is defined by the beam divergence and can 

be estimated as \[t>Zi. This means that the rate of change of the radiation 

spot size from the lens to the waist can be estimated as \[t)z\l f ., since di ~ /. 
The parametric dependence of the depth of focus can be found requiring that 

the spot size increase due to a displacement Ai, that is {^\[t)z\j be 



of order of the radiation spot size at the waist, i.e. |M|yA'". This yields the 
following Geometrical Optics prediction for the depth of focus: 



^ . (341) 



One may separately estimate the diffraction size related with an aperture 
of size \[t>zi (the radiation spot size at the non-limiting pupil aperture) at 
the observation plane. Such estimation yields a diffraction size di/{\[t)zi) ~ 
f /{\[t)zi). Now, requiring that the spot size increase due to a displacement 
A'5, that is (^\[t>zi/ A'i, be of order of the diffraction size f /{\[t)zi) one 
may estimate the diffraction limited depth of focus as 

A'z ~ ^ . (342) 
Dzi 



The following comparison between the diffraction limited depth of focus A'z 
and the Geometrical Optics prediction for the depth of focus Az can be pre- 
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sented: 



Eq. (343) allows one to conclude that diffraction effects can be neglected within 
the accuracy of the quasi-homogeneous approximation, and that Eq. (341) is 
a correct estimation of the depth of focus in this case. Note that the definition 
in Eq. (338) has no meaning in the case of a condenser system with a quasi- 
homogeneous source. Finally, it is interesting to remark that the depth of focus 
is much shorter than the distance A between the waist (image) and the focal 
plane. In fact 



1 



A \d \M\f \DzI 



N 



< 1 , (344) 



because of the far zone assumption {A <^ D). 



13.1.2 Near zone 



Similar estimations can be made in the near zone, i.e. assuming A ^ D., 
when the waist is near the focal plane. On the one hand, as has been already 

remarked, the radiation spot size at the focal plane is about \[t)f and does not 
depend on TV. On the other hand, the radiation spot size on the non-limiting 
pupil aperture is defined by the electron beam waist and can be estimated as 

. This means that the rate of change of the radiation spot size from the 

lens to the waist can be estimated as / f since the waist is near the focal 
plane. The parametric dependence of the depth of focus can be found imposing 

that the spot size increase due to a displacement Az, that is ^V^//^ Az be 

of order of the radiation spot size at the waist, i.e. \/^f. This yields the 
following Geometrical Optics prediction for the depth of focus: 

Az ~ -^^^ . (345) 



One may separately estimate the diffraction size related with an aperture of 
size (the radiation spot size at the non-limiting pupil aperture) at the 
observation plane. Such estimation gives a diffraction size di/\[^ ~ f /\[^. 
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Now, imposing that the spot size increase due to a displacement Ai, that is 
^ViV //^ be of order of the diffraction size f /\[^ one may estimate the 
diffraction hmited depth of focus as 

A'zc^^ . (346) 



The following comparison between the diffraction limited depth of focus A'i 
and the Geomterical Optics prediction of the depth of focus Ai: 



Eq. (347) allows one to conclude that diffraction effects can be neglected within 
the accuracy of the quasi-homogeneous approximation, and that Eq. (345) is 
a correct estimation for the depth of focus in this case. Finally, it is interesting 
to remark that the depth of focus is much shorter than the distance A between 
the waist (focal plane) and the image plane. In fact 



Ai _ 1 

AT- ']M\f~\~J 



Dzi 



< 1 , (348) 



because of the near zone assumption [A^ D). 



13.2 Effect of aperture size 



We will now analyze the influence of a finite aperture on our analysis of the 
depth of focus. For simplicity we will consider the far field all our 

examples involving effects of finite aperture size have been relying on the far 

field assumption. At first glance, substitution of the radiation spot size \[^zi 
with the finite aperture size a solves all issues, as it allows estimations of the 
rate of change of the radiation spot size and of diffraction effects. This would 
lead to the following equations for the Geometrical Optics prediction of the 
depth of focus: 



M ~ ^-^^ , (349) 
Zia 
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for the diffraction limited deptfi of focus: 



^'z-St , (350) 



d2 



for the ratio between A'z and Az: 
A'z p zict, Zl 

rsj . _ 



(351) 



and for the ratio between the depth of focus and the distance A: 
Az /^V^ 1 ^ 

-r- — . . • ^ . (352) 

A zia |M|/ a ^ ' 



We have seen in Section 11 that there exist particular situations when the 
phase factor in parenthesis {...} in Eq. (321) can be neglected. 

In one of these cases, when Dz\ ^ N ^ a? ^ zf/N one recovers a scaled 
image of the virtual source at the focal plane, where the scaling factor //^i ^ 
|M| (see Eq. (325)). Thus, in this situation, it is not correct to state that the 
radiation spot size at the image plane is smaller than the radiation spot size 
at the focal plane, because in both these planes we have the same, identical 
image of the virtual source with approxiniatively the same scaling factor. It 
should be remarked that in this case the depth of focus Az is longer than the 
distance between focal and image plane. A, as one may see from Eq. (352) 
when N ^ a^. Another critical scale to account for is zf/N. As we have seen 
in Section 11, when <^ zf/N one obtains, at the focal plane, the diffraction 
pattern of the pupil (aberrated or not, depending on how scales with respect 
to Zl). This case corresponds to coherent illumination of the pupil. In this 
situation, the diffraction limited depth of focus A'z is much longer than the 
Geometrical Optics prediction of the depth of focus Az, as one may see from 
Eq. (351). In this case, Eq. (350) should be used instead of Eq. (349) in order 
to estimate the depth of focus of the system. Summing up, Eq. (349) can be 
used without any remarks in all the cases when ^ zf/N. When ^ and 
a"^ ^ zf/N one should remember that the depth of focus becomes longer than 
the distance between the image and the focal plane, this case corresponding the 
situation when an image of the virtual source is formed in the focal plane. In 
all cases when cP' <^ zf/ N instead, we have coherent illumination of the pupil. 
The diffraction limited depth of focus becomes longer than the Geometrical 
Optics prediction of the depth of focus and, for estimations, one should use 
Eq. (350) in place of Eq. (349). 
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14 Solutions to the image formation problem for non- homogeneous 
undulator sources 



We will now consider non-homogeneous undulator sources. The next Section 
14.1 will deal with the case of a horizontally quasi-homogenous and vertically 
non-homogeneous source. Special emphasis will be given to the study of this 
situation, which is relevant for the majority of Synchrotron Radiation sources 
of the third generation. In the following Section 14.2 we will discuss, instead, 
the case of a horizontally non-homogeneous and vertically diffraction limited 
source. These two cases practically deal with all third generation light sources, 
from the soft to the hard X-rays. 

Our treatment is based on the assumption that the beta functions have minima 
in the center of the undulator. Placing the virtual source in the center of the 
undulator, we will obtain a particularly simple expression for the cross-spectral 
density of the source and its Fourier transform. In particular, calculations will 
yield real results. Statistical Optics methods conjugated to Fourier Optics 
will allow us to give explicit presentations of the cross-spectral density and 
its Fourier transform at the virtual source position. Such presentations will 
account for all complexities of the source, intrinsic properties of undulator 
radiation and electron-beam phase space distribution. Based on our previous 
work [2], we will first calculate the cross-spectral density G in the far zone 
limit. Then, with the help of Eq. (127), we will be able to recover the Fourier 
transform of the cross-spectral density Q at the virtual source position. By 
inverse transforming that expression we will finally recover the cross-spectral 
density as well. 

In general, we cannot give an explicit expression for the cross-spectral den- 
sity at any observation plane. In other words, as before, Eq. (133) cannot be 
calculated explicitly for any value of z\ and Z2. However, once the center of 
the undulator is fixed as the virtual source position, there exists a privileged 
observation plane. This is the plane where the center of the undulator is im- 
aged and depends on the position zi chosen for the lens. Throughout this 
paper we simply called it the image plane. Starting from the expression for 
the cross-spectral density and its Fourier transform at the virtual source we 
will be able to calculate, for any choice of ^i, the cross-spectral density on the 
focal {z2 = Zf) and on the image {z2 = Zi) plane. The procedure that we will 
use to perform these calculations is similar to what has been proposed in the 
case of homogeneous sources and will take advantage of Eq. (131) and Eq. 
(132). 

As usual, we will begin our investigations neglecting pupil effects. Results for 
a large non-limiting aperture on the image plane will be generalized to include 
effects of the pupil with the help of Eq. (135). Finally, in Section 14.3 we will 
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Fig. 42. Plot of the universal function M, used to calculate the cross-spectral density 
at the focal plane when I, Dx ^ I, Ny and Dy are arbitrary. 

critically discuss the assumptions made throughout this paper and see how 
and up to what extent they can be relaxed. 



i^.i Horizontally quasi-homogenous and vertically non-homogeneous source 



14.. 1.1 Large non-limiting aperture 

Let us first retain the assumption I and iVa; :§> 1, but allow Dy and Ny 

to assume arbitrary values. Eq. (141) allows the reconstruction of the vertical 
cross-spectral density in the far zone. Note that, when both A^^^ < 1 and 
Dy < 1 the source is non-homogeneous. In this case the far zone limit is for 
values 5o ^ 1, and one obtains 



G{zo,9y,My)^e-x.^ 



i26yZoA6y 



exp 



-2NyA§l 
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{9y + 4>y) 



2A, 



(353) 



where the function M{a, S) is the normahzed version of a universal function 
first defined in [2] and reads: 



M{a,d) 



80F 



rSmc 



smc 



+ (a + 5y 



(354) 



A plot of the M function is given in Fig. 42. The M function is a real function. 
Another remarkable property of M is its invariance for exchange of a with 5. 
Also, M is invariant for exchange of a with —a (or 5 with —5). The following 
relations between universal functions hold: 



is{a)^M{a,Q) , 



m 



1 °° 



(355) 
(356) 



and 



7(xj 



1 ^ 

— I deexpH2xe]M(e,0) . 



(357) 



Using Eq. (127) and Eq. (353) we obtain the Fourier transform of the cross- 
spectral density at io = at the virtual source position 



g (0, u, Au) = exp \-2NyAu^ 



,exp 



2D,. 



(358) 



Inverse transforming Eq. (358) we obtain the cross-spectral density at the 
virtual source position 



G'(0,y,Ay) =exp[-2L>^Ay'] j ^exp 



{y + 0)' 

27V„ 



M{Ay,(t>). (359) 



In analogy with Eq. (126), we defined M. as 
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Fig. 43. Plot of the universal function A4, used to calculate the cross-spectral density 
at the image plane when ^ 1, Dx ^ ^, and Dy are arbitrary. 



M{x,y) = I "^"^ I dSe^p[-2i{ax + 5y)]M{a,S) , (360) 

— oo —oo 

where /C is given in Eq. (174). A plot of Xi is presented in Fig. 43. Note that 
the symmetry of M for exchange of a with —a or S with —S impUes that 
7W is a real function. This fact, together with Eq. (359), implies that G at 
the virtual plane, that is an ensemble- averaged field product, is real as well. 
This is a consequence of the fact that single particle sources produce a field 
characterized by a plane wavefront at the virtual source position. It can be 
seen comparing Eq. (36) with Eq. (94) or, directly, by inspecting Eq. (97). 
Moreover, note that Ai is invariant for the exchange of x with y. 

From Eq. (131) and Eq. (358) we obtain the cross-spectral density on the focal 
plane in the vertical direction, that is 
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G{zf,yf,Ayf)=exp 
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(361) 



In the vertical direction, the relative intensity on the focal plane is therefore 
given by 



iVf/f 



o. 



X 



exp 



2D„ 



2D„ 



-1 



(362) 



because Isi^y) — M{(f)y, 0). In the limit Dy <^1, Eq. (362) reduces to 



i{zf,yf)^is[H 



(363) 



The modulus of the spectral degree of coherence on the focal plane in the 
vertical direction reads 



li/(%,y/,^y/)l=exp 
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(364) 



In the limiting case for Dy <^ 1 one obtains the simplified expression 



\9{zf,yf,^yf)\^ew 



(365) 
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Fig. 44. Plot of the universal function x-, used to calculate the modulus of the 
spectral degree of coherence on the focal plane when iVj; S> 1, -Da; ^ 1; -^i/ ^ 1 and 

Dy < 1. 

where, as in Eq. (145) of [2], 



M(«,5) 



1/2 



Is{a - 8) Is{a + 5) 



1/2 



(366) 



A plot of X is given in Fig. 44. Eq. (364) reduces even more in the case when 
both Dy and TV^^ <^ 1. In this case we have 



\9{zf,yf,Ayf)\=x\ 



(367) 



On the image plane, Eq. (132) and Eq. (359) give the following cross-spectral 
density in the vertical direction: 
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With the help of Eq. (360), Eq. (356) and Eq. (174) one sees that 



M{0,y) 



27r2/C 



d6 



exp[i{-2S)y] / daM{a,S) 



B{y) . (370) 



The solution of the image formation problem is thus constituted by a con- 
volution product between a Gaussian function and the universal function B, 
which admits the analytical representation given in Eq. (186). The intensity 
on the image plane is independent of the value of Dy. In the limit Ny^l Eq. 
(369) gives back Eq. (169). Instead, in the limit iVj, -C 1 we obtain 



(371) 



The modulus of the spectral degree of coherence in the vertical direction is 
given by 



\g{zi,yi,Ayi)\= / d^exp 



xA^(mA^j,0)exp -2Dym^Ayi 
{myi + mAyy + 




M{0,4>) 



-1/2 



-1/2 



(372) 



In the case Dy and Ny <^ 1 one obtains 
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Fig. 45. Plot of the universal function X, used to calculate the modulus of the 
spectral degree of coherence on the image plane when S> 1, -Da; S> 1, ^ 1 
and Dy <^ 1. 

\g{zi, Vi, Ayi) \ = X (mAyj, myi) (373) 



where, in analogy with Eq. (366), we define 



M{x,y) 



1/2 



B{x — y) B{x + y) 



1/2 



(374) 



Note that when A^j^ S> 1 and Dy assumes arbitrary values, the modulus of the 
spectral degree of coherence, Eq. (372), simplifies to Eq. (170). A plot of the 
universal function X is given in Fig. 45. 



14.. 1.2 Effect of aperture size 

We will now include, with the help of Eq. (135), the effects of pupil in the 
one-dimensional case. The pupil function and V are given by Eq. (224) and 
Eq. (225). The r-direction should be now substituted with the y-direction. 
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We can use condition (136) and Eq. (368) to describe the case when the 
lens is in the far zone. From Eq. (368) we can estimate the order of the 

source dimension, that is max[l, \//^], and the order of the coherence length, 
that is min[l/-y^, 1]. The lens is in the far zone when max[l, \[^y\/ z\ <^ 

ma.x[\J Dy, 1]. Note that several particular cases of interest are automatically 
included in this condition: the case for Dy ^ 1 and Ny ^ 1, that gives 
zi » ^ Ny/Dy, the case for Dy ^ 1 and TV^^ < 1, that gives zi ^ 1/^fWy, 

as well as the case for Dy <1 and Ny 1, that gives Zi ^ \//V^- All these 
situations have been previously discussed in Section 5. The new situation left 
to consider is for A^^ ~ 1 and Dy ~ 1. In this case, the lens is in the far field 
when ii » 1. Prom Eq. (135) and Eq. (368) we obtain 



Gp{zi,yi,Ayi 
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(375) 



According to Eq. (353), zf max[l, Dy] is of order of the square of the radiation 
spot size on the pupil, while zf/ max[7Vj^, 1] is of order of the square of the 
coherence length on the pupil. It is interesting to see that if Dy S> 1 the 
exponential function in Au has a very narrow characteristic length with respect 
to unity, and one can make the substitution Ai{ziAu, (py) — > J\4{0,(j)y) = 
B{(l)y), thus getting back Eq. (256). If iV^ :§> 1 instead, 4>y can be set to zero 
in the exponential function in (py + ziu, and the entire exponential function 
can be taken out of the integral in d(t)y. Then, it is possible to show that the 
surviving integral in d(l)y is equal to 7(^1 Art), giving back Eq. (250). Finally, 
if both TVj, > 1 and Dy :$> 1 we get back Eq. (226). 

Eq. (375) should be considered as a limiting case when no aberrations are 
present. More generally, when aberrations and apodization are present, Eq. 
(375) should be substituted by the following expression: 



Gp{zi,yi,Ayi)=4:d^exp 
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(376) 



Eq. (376) is very general and is valid under the only assumption that the 
lens is in the far field region. The use of old coordinates yn and yi2 instead 
of yi and A^j somewhat clarifies the meaning of Eq. (376). Eq. (376) states 
that the cross-spectral density accounting for the pupil influence (aberrations 
and diffraction effects) is a double convolution product of Eq. (368), i.e. the 
cross-spectral density in the ideal case, and Va, which can be written as 



Gp{zi,yn,yi2)= {Wii)E*{u2))*A{ui)*V:{u2) {Vium) 



E*Va {yn) 



E* * v: 



{yi2)) ■ 



(377) 



The intensity at the image plane is found setting yn = yi2 = yi, which gives 



Ip{zi,yi)^ {E{u,)E*{u2)) *Va{ui) *V:{U2)\ {yi,yi) 

E*Va\m)) ■ (378) 



In particular, in the case of a completely coherent source, we may neglect the 
ensemble average. 

Note that in order to evaluate the intensity at the image plane, it is no more 
enough to know the ideal intensity and to convolve with a line spread function. 
One has to know the cross-spectral density in the ideal case and convolve twice 
with Va, which is known as the amplitude line spread function of the system 
(or the amplitude point spread function in the two-dimensional case) and is 
a more general identifier of the system characteristics. Even in the non quasi- 
homogeneous case one may continue to use, for evaluating the intensity at the 
image plane, an algorithm based on the calculation of ideal characteristics and 
further convolution with a function characterizing the system. The difference 
with respect to the quasi-homogeneous case is that the amplitude line spread 
function must be used in place of the line spread function, and that the cross- 
spectral density must be used in place of the intensity. 

Note that, in our approach, the Wigner function plays no role, as it con- 
stitutes an artificial quantity in the image formation problem, whereas the 
natural quantity to consider is the cross-spectral density. Even if the virtual 
source can be characterized by a phase space distribution (i.e. by a positive 
Wigner function) one has no simplification in the imaging problem. The only 
simplification of Eq. (378) takes place in the quasi-homogeneous case when, 
due to the separability of the cross-spectral density variables and to a short 
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coherence length (compared with the size of the source) , one obtains the usual 
incoherent line spread function formalism. 



14-2 Horizontally non-homogeneous and vertically diffraction limited source 



We will now relax the assumption of a large horizontal electron beam size 
and divergence and deal with the non-homogeneous case when assumes 
arbitrary values, while <^ 1. This is a rather exotic range of parameters, 
and we will discuss the case for a large non-limiting aperture only. In this case, 
for third generation light sources we have automatically iVj, <^ 1 and Dy <^ 1 
because Cy <^ e^^* . In this case, in the limit Zo ^ 1, from Eq. (124) we obtain 



G = exp (i29^A9^Zo) exp [-2N^A9l 
'{9,-A9,f + {9y-A9y) 



xsmc 



exp (i29yA9yZo^ 

{9, + A9,)^ + {9y + A9yf 



smc 



(379) 

Using Eq. (127) and Eq. (379) one obtains the Fourier transform of the cross- 
spectral density at Zg — 0, i.e. at the virtual-source position, that is given 

by 



g(0,9^,A9^,9y,A9y)^exp 



-2N,A9i 



smc 



{9,-A9^f + {9y-A9yf 



xsmc 
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(380) 



Inverse transforming Eq. (380) it is possible to express the cross-spectral den- 
sity at the virtual source position as 



(7(0, X, Ax, y, Ay) = j d9^dA9^d9ydA9y exp [-2N^A9l 



This treatment can be easily generahzed to the more exotic situation for Nx 
arbitrary, <C 1, Ny arbitrary and Dy <C 1. However, here we will be concerned 
with third generation light sources only. Assuming reasonable values for Px ^ 10 L^, 
and €y = x^x with x — 10^^ we see that Dx <SC 1 implies both Ny <^ 1 and Dy <^ 1. 
Therefore we will avoid to make generalizations which are not pertinent to the case 
under study, e.g. exotic case for (3x > ^0 Lyj, (3y < 10~^ Lyj or (3y > 10 L^. 
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From Eq. (131) and Eq. (380) we obtain the cross-spectral density on the focal 
plane 
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(382) 



The relative intensity on the focal plane is therefore given by 



A%>^/>Z//)=sinc^ 



(383) 



This is just the relative intensity on the focal plane from a single electron, i.e. 
Eq. (91). It is interesting to note that the modulus of the spectral degree of 
coherence on the focal plane depends on Axf only, and can be written as 



|5(%,a;/,A%,2//,A^j)|=exp 

(384) 

In the limit iVa; <^ 1 one recovers the deterministic case of a single particle. In 
this limit \g\ reduces to unity, and the wavefront is perfectly coherent. 

As regards the image plane, we should note that Eq. (381) is not easy to 
manipulate analytically in the most general case. However, when Ax = and 
= 0, one can calculate the intensity of the virtual source. With the help 
of Eq. (175) and Eq. (181) one obtains 
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i{Q,x,y)= j d(f)xdexp 



(385) 



where we have set 



(386) 



The function ^ was ahcady defined in Eq. (95). Using Eq. (132) and Eq. (385) 
we can now give the following expression for the relative intensity: 



^(0a;,niyi) 




In the limit -C 1 we have 



^(0.,O) 



(387) 



i{zi,Xi,yi)^m^'^{mXi,myi) , 
in agreement with Eq. (96). 



(388) 



14-3 General imaging considerations 



In the present Section 14.3 we discuss a general algorithm for the solution to 
the image formation problem for undulator sources based on our Statistical 
Optics approach. 

Eq. (376) is an expression for the cross-spectral density on the image plane 
in the case of a non-homogeneous undulator source and of a lens with an 
arbitrary pupil function (i.e. a lens with aberrations, apodization and finite 
aperture size). However, we assumed that the electron beam has (i) a Gaussian 
transverse profile and (ii) a large horizontal emittance compared with the 
radiation wavelength (A^^: 3> 1 and S> 1), that (iii) the radiation frequency 
is tuned at perfect resonance with the fundamental frequency of the undulator, 
i.e. C <C 1 ^'^ , that (iv) the minimal beta functions in both horizontal and 



This means that monochromatization is good enough to neglect finite bandwidth 
of the radiation around the fundamental frequency, as weU as electron beam energy 
spread. 
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vertical directions are located at f = 0, that (v) there is no influence of 
focusing inside the undulator, that (vi) the lens is placed in the far zone and, 
finally, that (vii) the observation plane is located at position z = Zi, where 
the virtual source (that we assume at i = 0) is imaged. Assumptions (i), 
(n), (in), (iv) and (v) are related to the form of the cross-spectral density 
at the virtual source plane, Eq. (359). They are very often, but not always 
verified. Moreover they do not depend on the particular imaging setup related 
with a given photon beamline. Assumptions (vi) and (vii) instead, are related 
with the imaging setup i.e. with how the lens and the observation points are 
positioned. 

The majority of the assumptions from (i) to (vii) are often verified for third 
generation light sources. As a matter of fact, our theory is specifically built to 
deal with third generation light sources. However, a generalization to include 
the case of spontaneous undulators installed in XFEL facilities (see e.g. [25], 
[26]) is certainly desirable. Of all restrictions from (i) to (vii), (v) is the more 
difficult to be relaxed. In order to do so, one needs to modify the expression 
for the single particle field to account for the infiuence of focusing inside the 
undulator. The other assumptions may be more easily relaxed, to give a more 
general algorithm for the calculation of the cross-spectral density in the case 
of arbitrary position of the lens (near or far zone). The case of spontaneous 
undulators installed in XFEL facilities is a particular one when it is needed to 
deal with both horizontal and vertical electron beam emittances comparable 
or smaller than the radiation wavelength. In fact, e-^ cri Sy ~ 0.3A. Results 
based on a Gaussian model of the electron beam with generic N^^y and D^^y 
may be presented. To this purpose one may use Eq. (124) to calculate the 
cross-spectral density in free space. As one may see inspecting Eq. (124), it 
is no more possible to separate the horizontal and the vertical direction. Eq. 
(124) is still subject to assumptions (i), (iii), (iv) and (v). Our theory is built 
by exploiting many simplifying assumptions. When some of them fail, one 
should go back to the point where the invalid simplification is exploited, and 
use a more generic expression in its place. In particular, while assumption (i) 
is quite realistic in the case if storage ring sources, it is to be regarded as a 
conventional assumption when treating sources based on linear accelerators. 
Our most generic expression, Eq. (114) should be used in place of Eq. (124) 
if one wishes to relax assumption (i), as well as (iii) and (iv). Once the cross- 
spectral density in free space is known through Eq. (124) or Eq. (114), one 
may account for a generic observation plane, thus relaxing assumption (vii). To 
this purpose one needs to put attention on the fact that any observation plane 
located at position z — Z2 is related with a certain plane at position z — Zg in 
front of the lens through the lens- maker equation, Eq. (21). The next step in 
our algorithm consists in finding the cross-spectral density ai z = Zg, which 
will be imaged at our chosen observation plane, located at i = ^2- To this 
purpose, Eq. (128) may be used. In fact, Eq. (128) describes the propagation 
of the Fourier transform of the cross-spectral density, in free space. As we 
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have seen, the fact that the lens is placed in the far field allows cancellation 
of one phase factor in the relation between the cross-spectral density on the 
image plane and the cross-spectral density on the virtual source plane. This 
can be seen using condition (136) in Eq. (132). Such phase factor should be 
retained if the lens is in the near zone (that is when condition (136) is not 
satisfied). In this way, assumption (vi) can be relaxed as well. At this point, 
the cross-spectral density G{zs,f,Ar) is known. The final step consists in the 
calculation of the amplitude point spread function Va of the system. The cross- 
spectral density at position z = Z2 is found convolving twice the product of 
the cross-spectral density at the source plane and the extra phase-factor due 
to the failure of condition (136) with the amplitude point spread function. 
This gives 



where m is now defined as m = {zi — Zs)/{z2 — zi), according to Eq. (22). 

It should be noted that, in this paper, together with a theory for third gen- 
eration light sources we also developed a particular language that can be 
applied for a wider range of problems. Up to now, research works dealing 
with transverse coherence properties of Synchrotron Radiation have used the 
standard language developed to treat Statistical Optics problems. Such a lan- 
guage has a very limited scope because Statistical Optics has mainly dealt 
with thermal-like sources. As a result, a time domain approach has often been 
used. Quasi-stationary approximation and ergodicity are usually assumed, so 
that time averages are used instead of ensemble averages. Then, the concept 
of cross-spectral purity [3, 8] must be forcefully evoked in order to separate 
longitudinal and transverse coherence effects, which are described through 
the mutual intensity function. Such language though, is not suitable to de- 
scribe Synchrotron Radiation experiments, where many radiation pulses are 
collected and results are averaged over an ensemble. Our approach starts from 
the very foundation of Statistical Optics, thus avoiding inconvenient assump- 
tions. A language based on a frequency domain description and on ensemble 
averages over an ensemble of radiation pulses (each corresponding to a differ- 
ent electron bunch) has been developed. In our paper we aimed at presenting a 
satisfactory description of the physics involved in the characterization of light 
sources from a statistical viewpoint, trying to exhaustively explain where the 
main ideas come from and where they lead to. We restricted our attention 
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to a quantitative treatment of third-generation light sources but we also laid 
the foundations to describe other kind of radiation sources, e.g. spontaneous 
undulator sources installed in XFEL facilities. In fact, many of the features 
of the relatively specialized setup considered in this work are common in the 
general theory of undulator sources. 



15 Supplementary remarks on quEisi- homogeneous undulator source 
asymptotes 



In the last Section 14 we presented quite general results for vertically non- 
homogeneous sources. Under the only assumption of a large horizontal emit- 
tance {D^ ^ 1 and A^^, ^ 1), Eq. (353) specifies the cross-spectral density 
of the radiation in the far zone and in the vertical direction without con- 
straints on Ny and Dy. Under the same conditions, Eq. (359) specifies the 
cross-spectral density of the virtual source in the vertical direction. Note that 
Eq. (359) is a convolution. When the dependence of G on can be isolated 
in a single factor, the source is forcefully quasi-homogeneous, while there are 
no cases when it is described by a more generic kind of Shell model (i.e. 
G ~ I (yi)^ I {1/2)0 (Ay)). It is important to stress this fact in connection 
with several research works [14, 15] devoting particular attention to the rela- 
tion between the Gaussian-Shell model and undulator sources. As we remarked 
in Section 7, Shell models (and, in particular, Gaussian-Shell models) may cer- 
tainly be useful for describing light sources other than undulator-based ones 
and for educational purposes, but they do not describe any practical realiza- 
tion of undulator radiation sources. 

Prom Eq. (353), setting A^^^ = 0, one obtains the intensity distribution in the 
far zone 



,exp 
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having used Eq. (355). From Eq. (359), setting Ay = 0, one gets the intensity 
distribution of the virtual source in the center of the undulator , 



l{0,y) = / d^exp 



iVy + 
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(391) 



^® We are always assuming that the minimal beta function of the electron beam is 
located at the center of the undulator. 
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having taken advantage of Eq. (370). Is in Eq. (390) and B in Eq. (391) 
are the universal functions given in Eq. (160) and Eq. (186). Both Eq. (390) 
and Eq. (391) are convolutions, and are valid regardless the values of A'j^ and 
Dy, i.e. regardless the fact that the source is quasi-homogeneous or not. In 
the case of a large non-limiting aperture and an ideal lens, other two exact 
results can be found in Section 14, which are independent of the values of 
A^j, and Dy. In fact, Eq. (362) and Eq. (369) prove that the exact expression 
for the intensity distribution is, both on the focal and on the image plane, a 
convolution between Gaussian and universal functions. 

When the source is quasi-homogeneous, with an accuracy scaling as the inverse 
number of modes, l/y^maxfTYj^, 1] max.[Dy, 1], we may take the approximation 
A^(Ay,0) ~ 7(Aj/)B(0) in Eq. (359). 

This fact may be demonstrated as follows. First, let us introduce a normalized 
version of the one-dimensional inverse Fourier transform of the function M, 
that is 



1 °° 

M'(Ay, A^;) M (Oy, AOy) exp [-2iAy^, 



— oc 
oo 



A 



M {Ay, y) exp 2iMyy dy . 



(392) 



The normalization factor A in Eq. (392) is defined as 



A= / M{Q,y)dy , 



(393) 



so that M'(0, 0) = 1. The cross-spectral density in Eq. (359) can therefore be 
written as 



G'(0,y,Ay)=exp[-24Af' 

oo 

X / (iwexp [— 2my] exp 



-2NyU^ 



M'{Ay,u) , 



(394) 



having used the convolution theorem. Under the quasi-homogeneous assump- 
tion, we can approximate M'{Ay,u) ~ M'{Ay,0)M'{0,u). To show this, let 
us represent M'{x,y) using a Taylor expansion around the point (0,0). One 
obtains 
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M\x,y) = l + Y. 



k\ 



X 



y 



x=0 



dy>' 



y=o 



+0{xy) , 



(395) 



where the normahzation relation M'(0,0) = 1 has been taken advantage 
of. Similarly, one may consider the following representation of the product 
M'{x, 0) M'(0, y) also obtained by means of a Taylor expansion: 



M'{x, 0)M'(0,y) 



X 



M'(0,0) + £ 
k=l 

OO 

M'(0,0) + ^ 



x'' d''M'{x, 0) 



^ k\ dx^ 
yUm'{Q,y) 



i=i 



OO 1 



7n /\ A' 



a; 



„ci-M'(x,0) 



2/=0. 



a;=0 



+ 0{xy) , 



(396) 



having used M'(0,0) = 1. Comparison of the last cquahty in (396) with 
the right hand side of Eq. (395) shows that M'{x,y) ~ M'{x,0)M'{0,y) up 



to corrections of order xy ~ 1/y max[A^j,, 1] max[Z)y, 1], that is the quasi- 
homogeneous accuracy. Using this approximation in Eq. (394) yields 



6(0, y, Ay) = exp [-2byAf] M'{Ay, 0) 

OO 

X J du exp [—2iuy] exp 



M'(0,u) . 



(397) 



Finally, recalling the definitions of 7 and B we can write Eq. (397) as 



00 

G{0,y,Ay)^exp[-2DyAf] ^{Ay) j #exp 



Eq. (398) is valid in any quasi-homogeneous case. 



{y + 4>f 

2K, 



B{$) . (398) 



Note that Eq. (398) accounts for diffraction effects through the universal func- 
tions 7 and B. This may be traced back to the use of the inhomogeneous wave 
equation to calculate the cross-spectral density for the virtual source, from 
which Eq. (398) follows. Deriving Eq. (398), we assume a large number of 
modes, and this justifies the use of phase space representation as an alterna- 
tive characterization of the source, in place of the cross-spectral density (i.e. 
Eq. (398) itself). 
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Setting Ay = 0, Eq. (398) gives the exact intensity distribution at the vir- 
tual source, i.e. Eq. (391). The spectral degree of coherence on the virtual 
source is then recovered using the definition of quasi-homogeneous source 
G — I{y)g{Ay). Since the source is quasi-homogeneous, the Fourier trans- 
form of the spectral degree of coherence g{Ay) yields the intensity in the 
far zone. Remembering that 7 and Is form a Fourier pair, we conclude that, 
starting from Eq. (398) it is possible to reproduce the exact result for the 
intensity in the far zone, Eq. (390). Quite remarkably, Eq. (398), which is 
derived under the quasi-homogeneous approximation and is related to an ac- 
curacy l/y^maxfiVy, 1] max[Z)y, 1], yields back two results, Eq. (390) and Eq. 
(391) which are vahd regardless the fact that the source is quasi-homogeneous 
or not. Moreover, in the case of perfect optics and non-limiting pupil aperture, 
and independently of the quasi-homogeneous assumption, the intensity profile 
in the virtual plane reproduces the intensity profile in the image plane, while 
the intensity profile in the far zone reproduces the intensity profile in the focal 
plane. Therefore, we can also conclude that Eq. (398) gives both the intensity 
in the focal and in the image plane for an ideal lens. Note again that also these 
two results, Eq. (362) and Eq. (369), have perfect accuracy. They are exact and 

are not subject to the quasi- homogeneous accuracy l/-ymax[/Vj^, 1] max[Z)j,, 1]. 

Let us now consider the particular quasi-homogeneous case when both /V^ :§> 1 
and ^ 1. In this case, the number of modes along the virtual source is of 
order \J NyDy and the normalized coherence length can be estimated as ~ 

\j \j^y <^ 1. In the case of non-ideal optics, once a line spread function I for 
the system is found, one may obtain the intensity distribution of the radiation 
by convolving / with the ideal image. The accuracy of these calculations is now 

the accuracy of the quasi-homogeneous assumption, 1 / \J NyDy. In the case a 
pupil with aperture a is present, we concluded in Section 8.1 that it makes 

sense to account for diffraction effects when a <^ \[T)zi. In this case in fact, 
the ratio between the width of the line spread function and the width of the 

ideal image is of order zi/ {a^JWy) 3> 1/ \j NyDy. As a result, when Eq. (398) is 
used to calculate the ideal intensity, it makes sense to account for diffraction 
effects. Worsening the accuracy in the calculation of the cross-spectral density 
of the source, we may reduce Eq. (398) to 



G(0,y,A^) =exp -2DyAy 



exp 



2N, 



(399) 



that is Eq. (146). Note that neglecting the product with the 7 function can be 

done with an accuracy 1 / \/^y, while extraction of the exponential function 
in y from the convolution product with the B function can be done with an 

accuracy 1/JWy. In our study case when Dy ^ 1 and Ny ^ 1, the overall 
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accuracy of Eq. (399) (or Eq. (146)) can be estimated as max(l/ y Dj^, 1/ y Ny), 
that is the accuracy of the Gaussian approximation. Such accuracy is much 
worse than that of the quasi-homogeneous assumption in Eq. (398), that is 



1/ \l NyDy. This fact has interesting consequences. In fact, Eq. (399) can first 
be used to calculate the ideal intensity on the image plane and, then, it may be 
convolved with the line spread function of the lens to give a characterization 
of the intensity distribution with a reduced accuracy. If, for instance, a <^ 



Dzi we can have situation when zi/{dy Ny) l/y NyDy but zi / (ay A^^ ) -C 
l/^fWy and, as a result, accounting for diffraction effects would not modify 
the intensity with accuracy 1 / In spite of this, going back to Eq. (398) to 



calculate the ideal intensity with better accuracy 1/ J NyDy, diffraction effects 



will appreciably modify the intensity within the accuracy 1/J NyD 



When Ny 1 and Dy 



1 tlie accuracy of the quasi-homogeneous ap- 



proximation becomes 1/^ Ny max[l , Dy\ . When TVj^ ~ 1 and Dy 1 it be- 



comes, instead, 1/ ymax.[l, Ny\Dy. In these cases, the accuracy of the quasi- 
homogeneous approximation is comparable to the accuracy of the Gaussian 
approximation. To be specific, when A^j, ^ 1 and Dy ~ 1 the accuracy of the 

quasi-homogeneous approximation is of order \/ J Ny (note that the coherence 



length at the pupil is 
substituted with 



1, that is the diffraction size) and Eq. (398) can be 



G'(0,y,Ay) =exp 



2N, 



exp 



-2DyAf 7(Ay) 



(400) 



without loss of accuracy, because the relative accuracy of the convolution is 
of order 1/ \f^y as the accuracy of the quasi-homogenous approximation. Eq. 
(400) is just Eq. (162). A similar reasoning can be done when Dy ^ 1 and 
Ny ~ 1. In this case the accuracy of the quasi- homogeneous approximation is 

of order 1/^J^y (note that the coherence length at the pupil is "C 1, that 
is much smaller than the diffraction size), and Eq. (398) can be substituted 
with 



(7(0, y. Ay) = exp -2byAy 



>exp 



(y + 

2N, 



(401) 



without loss of accuracy. In fact, neglecting the 7 function in Eq. (398) is 
equivalent to approximate the convolution in Eq. (390) with a Gaussian dis- 
tribution, which can be done with an accuracy of order Dy, the same of 
the quasi-homogenous approximation. Eq. (401) is just Eq. (188). 
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In closing this Section we should stress that, in the case of ideal lenses, the 
intensity on the image plane, that is a scaled version of Eq. (391), is exactly 
given by a convolution of a known universal function and the electron beam 
profile. In the case the electron beam profile is unknown, one may measure 
the intensity and deconvolve Eq. (391) in order to find back the electron beam 
profile. In literature (see, for example, [27]) it is usually accepted that the 
resolution of any electron beam size a inferred from the measurement of the 
radiation intensity distribution on the image plane is limited (in the case of an 
ideal lens) by the diffraction size of the single particle undulator radiation, i.e. 
the resolution is of order \f\L^I {27ra). However, we have seen that Eq. (391) 
is an exact result. Therefore, any measurement of a obtained by deconvolution 
of Eq. (391) is only limited by the finite accuracy of the detector. 



16 Conclusions 

As has been remarked in [28]: " [...] it is very desirable to have a way to model 
the performance of undulator beamlines with significant partial coherent ef- 
fects, and such modeUing would, naturally, start with the source. The calcu- 
lation would involve the knowledge of the partial coherence properties of the 
source itself and of how to propagate partially coherent fields through space 
and through the optical components used in the beamline. [...] it is impor- 
tant to recognize that, although most of these calculations are, in principle, 
straightforward applications of conventional coherence theory (Born and Wolf, 
1980; Goodman, 1985), there is not much current interest in the visible optics 
community. [...] For example, even for the rather simple problem of diffrac- 
tion by an open aperture with partially coherent illumination, we have found 
published solutions only for circular and slit-shaped apertures and only for 
sources consisting of an incoherently illuminated aperture of similar shape to 
the diffracting aperture. Thus, there is no counterpart in these types of Fourier 
Optics problem to the highly developed art of ray tracing in geometrical op- 
tics, not is there anything as simple as a ray to which an exact system response 
can be calculated.". This program of development of Synchrotron Radiation 
theory was formulated more than ten years ago. Operation of third generation 
fight sources also started in this period. This demonstrates that when third 
generation light sources were born, it was immediately recognized that the 
usual theory of Synchrotron Radiation was not adequate to describe them. 
Yet, up to now, no theoretical progress has been made in that direction. The 
present paper, as well as our previous work [2] are devoted to the realization 
of the before mentioned program of development of Synchrotron Radiation 
theory. 

In [2] we described spatial coherence propertied of undulator radiation from 
third-generation light sources in free space. In this paper we aim at an exten- 
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sion of [2]. PrcvioTis scientific works and textbooks postulate that the cross- 
spectral density at the virtual undulator source can be described in terms of 
a Gaussian-Shell model or, even more restrictively, that undulator sources are 
perfectly incoherent. Such assumption is not adequate when treating third 
generation light sources, because the vertical emittance is comparable or even 
much smaller than the radiation wavelength (i.e. <C A/(27r)) in a wide spec- 
tral interval extending from the Amstrong wavelength range up to the soft 
X-rays. 

In this work we combined Statistical Optics methods with Fourier Optics 
techniques in order to describe in an analytical way the propagation of the 
cross-spectral density of Synchrotron Radiation through a lens. In particular, 
we focused our attention on the problem of finding both the intensity and the 
spectral degree of coherence of undulator radiation at the focal and at the 
image plane of the lens. Although our paper is not limited to this situation 
alone, our main result deals with the quite generic case of a large normalized 
horizontal emittance ^ X/{2t:) and an arbitrary vertical emittance Cy. 

Our paper provides physical understanding of a setup of general interest and 
we expect it to be useful for practical estimations in almost all range of the pa- 
rameter space for third generation light sources. We expect that, in the future, 
numerical codes fully capable of dealing with transverse coherence properties 
of Synchrotron Radiation will also be developed, and will be capable of pro- 
viding detailed analysis of particular experimental setups. Our theory will 
be of help to developers of these codes because it provides both benchmarks 
and partially manipulated equations for the field correlation, simpler to treat 
numerically than first principle calculations and still reasonably generic. 

Two basic non-restrictive assumptions made in our theory are the paraxial 
approximation and the resonance approximation. The first is justified by the 
fact that we are treating an ultra-relativistic system. The second means that 
we are working with an undulator composed of a large number of period 
Nyj ^ 1 and that we are interested in frequencies near the fundamental. This 
allows to neglect the vertical polarization component of the field and to treat 
the field within a scalar theory. 

Analytical studies also required the introduction of some restrictive assump- 
tions introduced in our theory, to be relaxed in the future. First, we assumed 
that the beta functions in both directions have their minima in the center of 
the undulator. However, beta functions are (or may be) often tuned around 
the center of the undulator in cases of practical interest. Second, we assumed a 
single converging thin lens with none or very specific pupil functions. However, 
other shapes of the lens can be accounted for by means of numerical convolu- 
tions between a more complicated complex pupil function and our fundamental 
results referring to the case when influence of the pupil is neghgible. Note that 
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the situation of a thin lens is often met in practice in the case of X-ray radi- 
ation, as grazing incidence reflective optics is quite frequently used for image 
formation in this spectral range. Third, we assumed that monochromatization 
is good enough to neglect finite bandwidth of the radiation around the funda- 
mental and also electron beam energy spread. Both these effects may be taken 
into account by an extension of this theory. However, in practical cases of in- 
terest these restrictions are often met. For instance, a monochromator relative 
bandwidth of 10~^ is sufficient to guarantee a reasonably narrow bandwidth 
for undulators up to about 40 periods. Also, with this number of undulator 
periods, such a monochromator guarantees small corrections to our results for 
electron beams with a relative energy spread of 10~^ or better. 

The three before mentioned assumptions are somewhat restrictive. However, 
they are often met in practice and our theory can be extended to the case when 
they are not met. A more restrictive assumption is that radiation frequency is 
tuned at perfect resonance with the fundamental frequency of the undulator. 
Although this last assumption is also often met in practice, our theory cannot 
be easily extended to account for the situation when it is not satisfied because 
the starting point for the calculation of the cross-spectral density is given by 
an expression for the electric field around the fundamental frequency. 

As a final remark, it should be said that we chose not to deal with bending 
magnet sources in this paper. This is left for future investigations. However, 
some estimates based on dimensional analysis suggest that, as discussed in 
Section 7, Geometrical Optics treatments may be sufficient to describe in a 
satisfactory way bending magnet radiation from third generation light sources. 

To conclude, our paper constitutes, to our knowledge, the first satisfactory 
theory describing imaging of undulator sources by a single non-ideal lens. We 
restrict ourselves to the analysis of a single lens for simplicity, the results 
for more complicated optical systems involving a larger quantity of optical 
elements being a straightforward extension of the present work. 
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